## Abstract

We analyze the far field and near field diffraction pattern produced by an amplitude grating whose strips present rough edges. Due to the stochastic nature of the edges a statistical approach is performed. The grating with rough edges is not purely periodic, although it still divides the incident beam in diffracted orders. The intensity of each diffraction order is modified by the statistical properties of the irregular edges and it strongly decreases when roughness increases except for the zero-th diffraction order. This decreasing firstly affects to the higher orders. Then, it is possible to obtain an amplitude binary grating with only diffraction orders -1, 0 and +1. On the other hand, numerical simulations based on Rayleigh-Sommerfeld approach have been used for the case of near field. They show that the edges of the self-images are smoother than the edges of the grating. Finally, we fabricate gratings with rough edges and an experimental verification of the results is performed.

© 2008 Optical Society of America

## 1. Introduction

Diffraction gratings are one of the most important optical components. It can be defined as an element which produces a periodical modulation in the properties of the incident light beam. Its behaviour in the near and far field is well known [1–7]. Amplitude or phase gratings are used in most applications. Also, other kind of gratings is possible, such as polarization gratings [8–10] or gratings with random microscopic irregularities in the topography [11–13]. In the far field the beam is divided into diffraction orders whose directions are given by the well-known grating equation [1]. The intensity of the diffraction orders is obtained as the square of the Fourier coefficients of the grating [2]. In the near field, Talbot effect is produced when the grating is illuminated with a plane wave [14, 15]. Self-images of the grating are formed at Talbot distances given by ${z}_{T}=\lambda /\{1-{[1-{(\lambda /p)}^{2}]}^{1/2}\}$ where λ is the wavelength of the incident wave and *p* is the period of the grating. When the period of the grating is much larger than the wavelength then the Talbot distance simplifies to *z _{T}*=2

*p*

^{2}/λ[16].

Theoretical approaches normally assume that the diffraction gratings present an ideal optical behaviour. However, this assumption is not always right. Flaws or defects can be produced, such as lost strips, during the fabrication process [17–21]. Other non-ideal behaviour is owing to roughness on the surface of the gratings, as it happens for steel tape gratings. Roughness produces a decreasing of the contrast of the self-images [11, 12]. Also, stochastical irregularities in the shape of the edges can be produced. This effect is not normally present in chrome on glass gratings or phase glass gratings, but strips with rough edges can be detected in some other manufacturing processes, such as laser ablation or chemical attack [13, 18].

In this work we theoretically, numerically and experimentally analyze the behaviour of amplitude gratings with rough edges. In particular, the far field intensity distribution and the Talbot effect are studied in detail. Talbot effect is a very important point which must be taken into account in devices that include a diffraction grating which works in near field approach. Rough edges effect can be very important in self-imaging phenomenon, then we will analyse it in detail. In addition, due to the stochastic properties of the edges, a statistical approach needs to be used. We theoretically show that the intensity of the diffraction orders depends on roughness parameters of the edges. Selecting the roughness level for the edges, we can suppress ±3 and higher orders, maintain orders 0 and reduce only slightly orders ±1.

For the near field, a theoretical analysis is not possible since the integrals cannot be solved analytically. Then, a numerical analysis based on Rayleigh-Sommerfeld formalism is performed for determining the properties of the self-images produced by the grating with rough edges [22]. The properties of the self-images are obtained for single realizations and also for statistical averages. In both cases, the edges of the self-images are smoother than the edges of the grating. Finally, we fabricate gratings with rough edges and we analyze the near field behaviour, which is in accordance to the theoretical results.

## 2. Far field approach

Let us consider an amplitude grating with period *p*. When the grating presents an ideal optical behaviour, the transmittance can be defined as a Fourier series *tp*(*ξ*)=∑_{n}
*a _{n}*exp(

*iqnξ*), where

*q*=2

*π*/

*p*and

*a*are the Fourier coefficients of the grating with

_{n}*n*integer. We assume that the edges of the strips are not straight, but present a certain random shape. Therefore, the grating cannot be described as its Fourier series expansion since it is not purely periodic, but as a sum of strips. To mathematically characterize the grating, let us assume that the left and right edges of a single strip are described as

*η*=

*fn(ξ)*and

*η*=

*gn(ξ)*, where (

*ξ,n*) are the transversal coordinates at the grating plane and

*n*=…-2,-1,0,1,2…. These functions are not analytical functions, but stochastic with the same statistical parameters. The exact shape of the edges

*f*and

_{n}(ξ)*g*is not known but it is possible to describe these functions by means of some stochastic parameters, such as their correlation length

_{n}(ξ)*T*and their standard deviation σ [23]. Let us assume that

*f*and

_{n}(ξ)*g*(ξ) present the same statistical distribution, that

_{n}*f*(ξ) and

_{n}*g*(ξ) are totally uncorrelated, and

_{n}*f*(ξ) and

_{n}*f*

_{n}^{′}(ξ

^{′}) are also totally uncorrelated except for

*n*=

*n*

^{′}. To describe the edges, a normal distribution with standard deviation

*σ*is assumed. We will also assume a Gaussian autocorrelation coefficient with correlation length

*T*. The two-dimensional transmittance of the grating

*t(ξ,η)*results

being $\prod \left({f}_{n}\left(\xi \right),{g}_{n}\left(\xi \right)\right)=\{\begin{array}{cc}1& {f}_{n}\left(\xi \right)<\eta <{g}_{n}\left(\xi \right)\\ 0& {g}_{n}\left(\xi \right)<\eta <{f}_{n+1}\left(\xi \right)\end{array}.$

An example of the grating proposed is shown in Fig. 1.

Light passes through the grating and then propagates a distance z. To determine the far field diffraction pattern, Fraunhoffer approach is used

being λ the wavelength of the incident light, *k=2π/λ*, *x* and *y* the transversal coordinates at the observation plane, *z* the propagation direction, and *U _{i}(ξ,η)* the incident field that, for simplicity, we will consider that it is a monochromatic plane wave in normal incidence,

*U*. The transmission function only affects to the integration limits into eq. (2). Then, the integral is converted into a summatory of finite integrals. Every one of these integrals corresponds to the contribution to the intensity of every strip of the grating, which has rough edges. Besides, we have truncated the summatory which appears in

_{i}(ξ,η)=U_{0}*t(ξ,η)*to be able to give a clearer equation for the intensity pattern. Thus, the intensity results in

$${\int}_{{g}_{n}\left(\xi \right)+\left(4n+1\right)p\u20444}^{{f}_{n}\left(\xi \right)+\left(4n-1\right)p\u20444}d\eta {\int}_{{g}_{n\prime}\left(\xi \prime \right)+\left(4n\prime +1\right)p\u20444}^{{f}_{n\prime}\left(\xi \prime \right)+\left(4n\prime -1\right)p\u20444}{e}^{-i\frac{k}{z}\left[x\left(\xi \prime -\xi \right)+y\left(\eta \prime -\eta \right)\right]}d\eta \prime $$

where *N* is the number of strips and *L* is their length. Performing the integrals in *η* and *η*
^{′} the intensity results

$$\left\{\begin{array}{c}{e}^{ik{\theta}_{y}p\left(n-n\prime \right)}{e}^{ik{\theta}_{y}\left[{f}_{n}\left(\xi \right)-{f}_{n\prime}\left(\xi \prime \right)\right]}+{e}^{ik{\theta}_{y}p\left(n-n\prime \right)}{e}^{ik{\theta}_{y}\left[{g}_{n}\left(\xi \right)-{g}_{n\prime}\left(\xi \prime \right)\right]}\\ -{e}^{ik{\theta}_{y}\frac{p}{2}\left[2\left(n-n\prime \right)+1\right]}{e}^{ik{\theta}_{y}\left[{f}_{n}\left(\xi \right)-{g}_{n\prime}\left(\xi \prime \right)\right]}-{e}^{ik{\theta}_{y}\frac{p}{2}\left[2\left(n-n\prime \right)-1\right]}{e}^{ik{\theta}_{y}\left[{g}_{n}\left(\xi \right)-{f}_{n\prime}\left(\xi \prime \right)\right]}\end{array}\right\},$$

where *θ _{x}=x/z* and

*θ*=

_{y}*y/z*. Assuming the stochastical description of the functions

*fn(ξ)*and

*g*given before, the characteristic functions for this distribution result in [23]

_{n}(ξ)$$\u3008{e}^{i\alpha \left[{f}_{n}\left(\xi \right)-{g}_{n\prime}(\xi \prime )\right]}\u3009={e}^{-{\left(\alpha \sigma \right)}^{2}},$$

$$\u3008{e}^{i\alpha \left[{f}_{n}\left(\xi \right)-{f}_{n}(\zeta \prime )\right]}\u3009=\u3008{e}^{i\alpha \left[{g}_{n}\left(\xi \right)-{g}_{n}(\xi \prime )\right]}\u3009={e}^{-{\left(\alpha \sigma \right)}^{2}}\sum _{m=0}^{\infty}\frac{{\left(\alpha \sigma \right)}^{2m}}{m!}{e}^{-\frac{m{(\xi -\xi \prime )}^{2}}{{T}^{2}}},$$

$$\u3008{e}^{i\alpha \left[{f}_{n}\left(\xi \right)-{f}_{n\prime}(\xi \prime )\right]}\u3009=\u3008{e}^{i\alpha \left[{g}_{n}\left(\xi \right)-{g}_{n\prime}(\xi \prime )\right]}\u3009={e}^{-{\left(\alpha \sigma \right)}^{2}},\left(n\ne n\prime \right).$$

being in our case *α=kθ _{y}*. Performing an averaging process in (4), using the relationships given in (5), and reorganizing the terms, then the average intensity is

$$-N{e}^{-{\left(k{\theta}_{y}\sigma \right)}^{2}}{\int}_{-L\u20442}^{L\u20442}{\int}_{-L\u20442}^{L\u20442}{e}^{ik{\theta}_{x}\left(\xi -\xi \prime \right)}d\xi d\xi \prime $$

$$+N{e}^{-{\left(k{\theta}_{y}\sigma \right)}^{2}}\sum _{m=0}^{\infty}\frac{{\left(k{\theta}_{y}\sigma \right)}^{2}}{m!}{\int}_{-L\u20442}^{L\u20442}{\int}_{-L\u20442}^{L\u20442}{e}^{ik{\theta}_{x}\left(\xi -\xi \prime \right)}{e}^{-\frac{m{\left(\xi -\xi \text{'}\right)}^{2}}{{T}^{2}}}d\xi d\xi \prime \},$$

where the averaging in the intensity is performed on a supposed ensemble of realizations which are obtained, for example, when the grating is moved in the direction parallel to the y axis. Let us assume that the length of the strips *L* is large, although not infinite. Then, the integrals have a simple analytical solution and the average intensity results

$$\sum _{j}{\mathrm{sinc}}^{2}\left[\pi \frac{p}{\lambda}\left(N+1\right)\left({\theta}_{y}-j\frac{\lambda}{p}\right)\right]$$

$$+\frac{2{\mid {U}_{0}\mid}^{2}\mathrm{NLT}\sqrt{\pi}{\sigma}^{2}}{{\left(\lambda z\right)}^{2}}{e}^{-{\left(k{\theta}_{y}\sigma \right)}^{2}}\sum _{m=1}^{\infty}\frac{{\left(k{\theta}_{y}\sigma \right)}^{2\left(m-1\right)}}{m!m}{e}^{-\frac{{\left(kT{\theta}_{x}\right)}^{2}}{4m}}.$$

Comparing this average intensity with that of an amplitude grating without roughness, *σ*=0,

we can normalize the Eq. (7) with respect to the maximum intensity of eq. (8), *I _{0}*=[|

*U*

_{0}|

*L*(

_{p}*N*+1)/(2

*λz*)]

^{2}. Also, since is normally very high, the sinc functions are very narrow and have significant values only when their argument,

*θ*, is nearly zero. Then the mean intensity results

_{y}-jλ/p$$+\frac{8\sqrt{\pi}T{\sigma}^{2}}{\mathrm{NL}{p}^{2}}{e}^{-{\left(k{\theta}_{y}\sigma \right)}^{2}}\sum _{m=1}^{\infty}\frac{{\left(k{\theta}_{y}\sigma \right)}^{2\left(m-1\right)}}{m!m}{e}^{-\frac{{\left(kT{\theta}_{x}\right)}^{2}}{4m}}.$$

where $\overline{\u3008I(x,y)\u3009}=\u3008I(x,y)\u3009\u2044{I}_{0}$ is the normalized average intensity and *a _{j}*=sinc(jπ/2). The first term of (9) corresponds to the diffraction pattern of an amplitude grating but multiplied by a factor which diminishes the intensity of the diffraction orders according to

In Table 1, the Fourier coefficients for the first five diffraction orders is shown for several values of *σ/p*.

On the other hand, the second term of Eq. (9) is produced only by the roughness of the edges, which affects in both directions, *x* and *y*. Roughness produces a Gaussian halo centred in the zero-th order. The width of the halo depends on *T* along the *x* direction and *σ* along the *y* direction. In Fig. 2, the far field diffraction pattern along the y-axis for different values of *T* and *σ* is shown. As it can be seen, when roughness increases, high diffraction orders disappear and the halo grows around the zero-th order. The width of the halo in the y-axis depends on *σ*. For higher values of *σ*, the width of the halo diminishes.

Let us analyze two important cases, such as high roughness and low roughness limits. The high roughness limit occurs when *σ*≫*λ*. The characteristic functions in (5) are still valid, except the autocorrelation function which is now [23]

where *T _{F}*=

*kT/σ*. With this substitution and performing the integrals, average intensity results

$$+\frac{8{T}_{F}\sqrt{\pi}{\sigma}^{2}}{\mathrm{NL}{p}^{2}}{e}^{-{k}^{2}\left[{\left({\theta}_{y}\sigma \right)}^{2}+{({\theta}_{x}{T}_{F}\u20442)}^{2}\right]}.$$

It is a more simplified result than Eq. (9), since the sum over j disappears.

On the other hand, the slight roughness limit occurs when *σ*≪*λ*. In this case, Eq. (5) is valid, except the autocorrelation function

Then, the average intensity results

Under this approach the halo disappears and the far field diffraction pattern is equivalent to that of a perfect grating with modified Fourier coefficients.

## 3. Self-imaging process in the near field

The intensity distribution in the near field can be obtained in a simple way replacing the Fraunhofer kernel with the Fresnel kernel in (3). Unfortunately it is not possible to obtain analytical solutions for these integrals. As a consequence, we have performed a numerical analysis to determine the characteristics of the intensity distribution in the near field. For the numerical implementation we have used a fast-Fourier-transform based direct integration method which uses the Rayleigh-Sommerfeld approach [22].

In first place we are interested in how the self-images of this grating with rough edges are formed. In Fig. 3 we show for comparison the self-images produced by a perfect grating, that is, without rough edges and in Fig. 4 and Fig. 5, two examples of gratings with rough edges and the first three self-images for different values of *σ/p*. Although the edges of the grating present a high roughness, the edges of the self-images are quite smooth. The reason is that Talbot effect is a cooperative effect since the intensity at a given point (*x,y*) of the image is obtained as an integration of the amplitude at the diffraction grating. It performs an averaging in the intensity distribution. In addition, an interferential process happens and produces a kind of speckle in the fringes. For the simulation we have considered a grating with size 150*µm*×300*µm*. Since the algorithm does not consider that the grating is periodic, an edge effect is produced. We show the central region of the intensity pattern to avoid this edge effect.

In Fig. 6, a comparison of the average profiles obtained with the perfect grating and the rough gratings is shown. In Fig. 7 (Media 1) a video is included where the transition between *z*=0 and the first self-image is shown. The intensity distribution at fractional Talbot planes is also shown in Fig. 7 (Media 1) for distances *z=z _{T}*/4,

*z*/3,

_{T}*z*/2, and also the average profile for these particular cases.

_{T}#### 3.1. Average intensity distribution

The grating can be placed in a mobile device and then the intensity pattern will be an average over a group of discrete intensity patterns. To characterize this, we calculate the near field intensity pattern for several realizations and then we perform an averaging in the intensity of these realizations. This procedure is repeated for different self images placed at *z=np ^{2}*/λ, with

*n*=1,2,…. The average intensity of these self-images is shown in Fig. 8 for an ensemble of 100 images. For this case, the self-images are very smooth.

## 4. Experimental approach

To confirm the validity of the results and to use a grating with roughness parameters known, we have manufactured a grating with rough edges using a direct laser photoplotter. The grating is an amplitude grating made of chrome on glass and its period is *p*=100*µm*. The grating is illuminated with a collimated laser diode whose wavelength is λ=0.65*µm*. In the near field approximation, some self-images have been acquired. For this, we have used a CMOS camera (ueye, pixel size: 6×6 microns) and a microscope objective in order to get a better resolution. In Fig. 9 (Media 2) we can observe the image of the grating using an optical microscope and the first three self-images taken with the CMOS camera. These images correspond to just one realization. As it can be seen, experimental results are in total accordance with the numerical results. The shape of the self-images is quite smooth compared to the shape of the strips edges. In the self-images we can also see a defect in one of the strips (rectangle) which gradually disappears as the order of the self-image increases. In Fig. 9 (Media 2) a video with the experimental images is also shown.

The mean profile of these self-images is also shown in Fig. 10 and also the intensity distribution of self-image 15. The intensity at *x*=200 *µm* distribution is very smooth except for a dust particle in the optics that we could not eliminate. Comparing this result with that shown in Fig. 6, we can validate the results given by the numerical analysis.

## 5. Conclusions

In this work, we have analyzed the far field and near field diffraction pattern produced by an amplitude grating whose strips present rough edges. Due to the stochastic nature of the grating a statistical approach is performed. In the far field, the intensity of the diffraction orders strongly decreases in terms of the roughness and the index of the diffraction order. Then, a possible application of this kind of gratings is to obtain amplitude binary gratings with only diffraction orders -1, 0 and +1. For the case of near field, an analytical result is not possible and numerical simulations based on a Rayleigh-Sommerfeld approach have been performed. The self-images are smoother than the grating, since Talbot effect is a cooperative effect. Finally, we have fabricated gratings with rough edges and an experimental verification of the theoretical and numerical results is performed.

## Acknowledgments

This work has been supported by the DPI2005-02860 project of the Ministerio de Educación y Ciencia of Spain and a CENIT project “Tecnologías avanzadas para los equipos y procesos de fabricación de 2015: e-eficiente, e-cológica, e-máquina (eEe)” of the Ministerio de Industria, Turismo y Comercio.

## References and links

**1. **M. Born and E. Wolf, *Principles of Optics* (Pergamon Press, Oxford, 1980).

**2. **J. W. Goodman, *Introduction to Fourier Optics* (McGraw-Hill, New York, 1968).

**3. **E. G. Loewen and E. Popov, *Diffraction gratings and applications* (Marcel Dekker, New York, 1997).

**4. **M. J. Lockyear, A. P. Hibbins, K. R. White, and J. R. Sambles, “One-way diffraction grating,” Phys. Rev. E **74**, 056611 (2006). [CrossRef]

**5. **S. Wise, V. Quetschke, A. J. Deshpande, G. Mueller, D. H. Reitze, D. B. Tanner, and B. F. Whiting, “Phase effect in the diffraction of light: beyond the grating equation,” Phys. Rev. Lett. **95**, 013901 (2005). [CrossRef] [PubMed]

**6. **C. Palmer, *Diffraction Grating Handbook* (Richardson Grating Laboratory, New York, 2000).

**7. **R. Petit, *Electromagnetic Theory of Gratings* (Springer-Verlag, Berlin, 1980). [CrossRef]

**8. **F. Gori, “Measuring Stokes parameters by means of a polarization grating,” *Opt. Lett.*24, 584–586 (1999).

**9. **C. G. Someda, “*Far field of polarization gratings*,” Opt. Lett. **24**, 1657–1659 (1999). [CrossRef]

**10. **G. Piquero, R. Borghi, A. Mondello, and M. Santarsiero, “Far field of beams generated by quasi-homogeneous sources passing through polarization gratings,” Opt. Commun. **195**, 339–350 (2001) [CrossRef]

**11. **F. J. Torcal-Milla, L. M. Sanchez-Brea, and E. Bernabeu, “Talbot effect with rough reflection gratings,” Appl. Opt. **46**, 3668–3673 (2007) [CrossRef] [PubMed]

**12. **L.M. Sanchez-Brea, F. J. Torcal-Milla, and E. Bernabeu, “Talbot effect in metallic gratings under Gaussian illumination,” Opt. Commun. **278**, 23–27 (2007). [CrossRef]

**13. **L. M. Sanchez-Brea, F. J. Torcal-Milla, and E. Bernabeu, “Far field of gratings with rough strips,” J. Opt. Soc. Am. A **25**, 828–833 (2008). [CrossRef]

**14. **W. H. F. Talbot, “Facts relating to optical science,” Philos. Mag. **9**, 401–407 (1836).

**15. **K. Patorski, “The self-imaging phenomenon and its applications,” Prog. Opt. **27**, 1–108 (1989). [CrossRef]

**16. **N. Guérineau, B. Harchaoui, and J. Primot, “Talbot effect re-examined: demonstration of an achromatic and continuous self-imaging regime,” Opt. Commun. **180**, 199–203 (2000). [CrossRef]

**17. **Y. Lu, C. Zhou, and H. Luo, “Talbot effect of a grating with different kind of flaws,” J. Opt. Soc. Am. A **22**, 2662–2667 (2005) [CrossRef]

**18. **P. P. Naulleau and G. M. Gallatin, “Line-edge roughness transfer function and its application to determining mask effects in EUV resist characterization,” Appl. Opt. **42**, 3390–3397 (2003). [CrossRef] [PubMed]

**19. **T. R. Michel, “Resonant light scattering from weakly rough random surfaces and imperfect gratings,” J. Opt. Soc. Am. A **11**, 1874–1885 (1994). [CrossRef]

**20. **V. A. Doroshenko, “Singular integral equations in the problem of wave diffraction by a grating of imperfect flat irregular strips,” Telecommunications and Radio Engineering **57**, 65–72 (2002)

**21. **M. V. Glazov. and S. N. Rashkeev, “Light scattering from rough surfaces with superimposed periodic structures,” Appl. Phys. B **66**, 217–223 (1998) [CrossRef]

**22. **F. Shen and A. Wang, “Fast-Fourier-transform based numerical integration method for the Rayleigh-Sommerfeld diffraction formula,” Appl. Opt. **45**, 1102–1110 (2006) [CrossRef] [PubMed]

**23. **P. Beckmann and A. Spizzichino, *The Scattering of Electromagnetic Waves from Rough Surfaces* (Artech House Norwood, 1987).