|
|
| Line 1: |
Line 1: |
| In [[optics]], the '''Fresnel diffraction''' equation for '''near-field diffraction''', is an approximation of [[Kirchhoff's diffraction formula|Kirchhoff-Fresnel diffraction]] that can be applied to the propagation of waves in the [[near and far field|near field]].<ref>[[Max Born|M. Born]] & E. Wolf, Principles of Optics, 1999, Cambridge University Press, Cambridge</ref> It is used to calculate the [[diffraction pattern]] created by waves passing through an aperture or around an object, when viewed from relatively close to the object. In contrast the diffraction pattern in the [[near and far field|far field]] region is given by the [[Fraunhofer diffraction]] equation.
| | I would like to introduce myself to you, I am Jayson Simcox but I don't like when people use my full name. Invoicing is my profession. It's not free tarot readings ([http://myfusionprofits.com/groups/find-out-about-self-improvement-tips-here/ check these guys out]) a common thing but what I like doing is to climb but I don't have the time lately. For a while I've been in Alaska but I will have are psychics real ([http://kjhkkb.net/xe/notice/374835 check these guys out]) to transfer in a yr or two.<br><br>Also visit my web site love psychics ([http://www.january-yjm.com/xe/index.php?mid=video&document_srl=158289 simply click the up coming website page]) |
| | |
| The near field can be specified by the Fresnel number, {{math|''F''}} of the optical arrangement, which is defined, for a wave incident on an aperture, as:
| |
| | |
| : <math>F = \frac{a^{2}}{L \lambda}</math>
| |
| where
| |
| : <math>a\!</math> is the characteristic size of the aperture
| |
| : <math>L\!</math> is the distance of the observation point from the aperture
| |
| : <math>\lambda\!</math> is the [[wavelength]] of the wave.
| |
| When <math>F \gtrsim 1</math> the diffracted wave is considered to be in the near field, and the Fresnel diffraction equation can be used to calculate its form.
| |
| | |
| [[Image:FresnelDiff9 3 PEM.jpg|250px|thumb|Fresnel diffraction showing central [[Arago spot]]]]
| |
| | |
| The multiple Fresnel diffraction at nearly placed periodical ridges
| |
| ([[ridged mirror]]) causes the [[specular reflection]]; this effect can be used for [[atomic mirror]]s.<ref> | |
| http://www.ils.uec.ac.jp/~dima/PhysRevLett_94_013203.pdf H. Oberst, D. Kouznetsov, K. Shimizu, J. Fujita, F. Shimizu. Fresnel diffraction mirror for atomic wave, [[Physical Review Letters]], '''94''', 013203 (2005).</ref> | |
| | |
| ==Early treatments of this phenomenon==
| |
| Some of the earliest work on what would become known as Fresnel diffraction was carried out by [[Francesco Maria Grimaldi]] in Italy in the 17th century.
| |
| In his monograph entitled "Light,"<ref>''Light," by Richard C. MacLaurin, 1909, Columbia University Press</ref> [[Richard C. MacLaurin]] explains Fresnel diffraction by asking what happens when light propagates, and how that process is affected when a barrier with a slit or hole in it is interposed in the beam produced by a distant source of light. He uses the Principle of [[Huygens–Fresnel principle|Huygens]] to investigate, in classical terms, what transpires. The wave front that proceeds from the slit and on to a detection screen some distance away very closely approximates a wave front originating across the area of the gap without regard to any minute interactions with the actual physical edge.
| |
| | |
| The result is that if the gap is very narrow only diffraction patterns with bright centers can occur. If the gap is made progressively wider, then diffraction patterns with dark centers will alternate with diffraction patterns with bright centers. As the gap becomes larger, the differentials between dark and light bands decrease until a diffraction effect can no longer be detected.
| |
| | |
| MacLaurin does not mention the possibility that the center of the series of diffraction rings produced when light is shone through a small hole may be black, but he does point to the inverse situation wherein the shadow produced by a small circular object can paradoxically have a bright center. (p. 219)
| |
| | |
| In his ''Optics,'',<ref>''Optics,'' Francis Weston Sears, p. 248ff, Addison-Wesley, 1948</ref> Francis Weston Sears offers a mathematical approximation suggested by Fresnel that predicts the main features of diffraction patterns and uses only simple mathematics. By considering the perpendicular distance from the hole in a barrier screen to a nearby detection screen along with the wavelength of the incident light, it is possible to compute a number of regions called half-period elements or [[Fresnel zone]]s. The inner zone will be a circle and each succeeding zone will be a concentric annular ring. If the diameter of the circular hole in the screen is sufficient to expose the first or central Fresnel zone, the amplitude of light at the center of the detection screen will be double what it would be if the detection screen were not obstructed. If the diameter of the circular hole in the screen is sufficient to expose two Fresnel zones, then the amplitude at the center is almost zero. That means that a Fresnel diffraction pattern can have a dark center. These patterns can be seen and measured, and correspond well to the values calculated for them. Figure 9-5, following p. 222, in Sears shows four Fraunhofer patterns in the top row followed by sixteen Fresnel diffraction patterns. Three of them have dark centers. (See the photograph above, and check the book by Sears for very much nicer photographs.
| |
| | |
| ==The Fresnel diffraction integral==
| |
| | |
| [[Image:Diffraction geometry.svg|thumb|350px|Diffraction geometry, showing aperture (or diffracting object) plane and image plane, with coordinate system.]]
| |
| The electric field [[diffraction]] pattern at a point ''(x,y,z)'' is given by:
| |
| | |
| :<math> E(x,y,z)={z \over {i \lambda}} \iint_{-\infty}^{+\infty}{ E(x',y',0) \frac{e^{ikr}}{r^2}}dx'dy' </math>
| |
| | |
| where
| |
| | |
| :<math> E(x',y',0) </math> is the aperture,
| |
| :<math> r=\sqrt{(x-x')^2+(y-y')^2+z^2} </math> , and
| |
| :<math> i \,</math> is the [[imaginary unit]].
| |
| | |
| Analytical solution of this integral is impossible for all but the simplest diffraction geometries. Therefore, it is usually calculated numerically.
| |
| | |
| === The Fresnel approximation ===
| |
| The main problem for solving the integral is the expression of ''r''. First, we can simplify the algebra by introducing the substitution:
| |
| | |
| :<math>\rho^2 = (x-x')^2+(y-y')^2 \,</math>
| |
| | |
| Substituting into the expression for ''r'', we find:
| |
| | |
| :<math> r= \sqrt{\rho^2+z^2} = z \sqrt{ 1 + \frac{\rho^2}{z^2} } </math>
| |
| | |
| Next, using the [[Taylor series]] expansion
| |
| | |
| :<math>\sqrt{1+u} = (1+u)^{1/2} = 1 + \frac{u}{2} - \frac{u^2}{8} + \cdots</math>
| |
| | |
| we can express ''r'' as
| |
| | |
| :<math>r = z \sqrt{ 1 + \frac{\rho^2}{z^2} } </math>
| |
| ::<math> = z \left[ 1 + \frac{\rho^2}{2 z^2} - \frac{1}{8} \left( \frac{\rho^2}{z^2} \right)^2 + \cdots \right] </math>
| |
| ::<math> = z + \frac{\rho^2}{2 z} - \frac{\rho^4}{8z^3} + \cdots</math>
| |
| | |
| If we consider all the terms of Taylor series, then there is no approximation.<ref>There was actually an approximation in a prior step, when assuming <math>e^{i k r}/r</math> is a real wave. In fact this is not a real solution to the vector [[Helmholtz equation]], but to the scalar one. See [[scalar wave approximation]]</ref> Let us substitute this expression in the argument of the exponential within the integral; the key to the Fresnel approximation is to assume that the third term is very small and can be ignored. In order to make this possible, it has to contribute to the variation of the exponential for an almost null term. In other words, it has to be much smaller than the period of the complex exponential, i.e. <math>2 \pi</math>:
| |
| | |
| :<math> k \frac{\rho^4}{8 z^3} \ll 2 \pi</math>
| |
| | |
| expressing ''k'' in terms of the wavelength,
| |
| | |
| :<math>k = { 2 \pi \over \lambda } \,</math>
| |
| | |
| we get the following relationship:
| |
| | |
| :<math> \frac{\rho^4}{z^3 \lambda} \ll 8</math>
| |
| | |
| Multiplying both sides by <math>z^3/\lambda^3</math>, we have
| |
| | |
| :<math> \frac{\rho^4}{\lambda^4} \ll 8 {z^3 \over \lambda^3}</math>
| |
| | |
| or, substituting the earlier expression for ρ<sup>2</sup>,
| |
| | |
| :<math> \frac{[(x-x')^2+(y-y')^2]^2}{ \lambda^4} \ll 8 {z^3 \over \lambda^3}</math>
| |
| | |
| If this condition holds true for all values of ''x'', ''x' '', ''y'' and ''y' '', then we can ignore the third term in the Taylor expression. Furthermore, if the third term is negligible, then all terms of higher order will be even smaller, so we can ignore them as well.
| |
| | |
| For applications involving optical wavelengths, the wavelength λ is typically many orders of magnitude smaller than the relevant physical dimensions. In particular:
| |
| | |
| :<math>\lambda \ll z</math>
| |
| | |
| and
| |
| | |
| :<math>\lambda \ll \rho</math>
| |
| | |
| Thus, as a practical matter, the required inequality will always hold true as long as
| |
| | |
| :<math>\rho \ll z </math>
| |
| | |
| We can then approximate the expression with only the first two terms:
| |
| | |
| :<math> r \approx z + \frac{\rho^2}{2 z} = z + \frac{(x-x')^2 +(y-y')^2}{2 z} </math>
| |
| | |
| This equation, then, is the '''Fresnel approximation''', and the inequality stated above is a condition for the approximation's validity.
| |
| | |
| ===Fresnel diffraction===
| |
| {| align="right" class="infobox" style="font-size:12px; text-align:center;"
| |
| |'''Fresnel diffraction''' occurs when:
| |
| <math>F = \frac{a^2}{L\lambda} \ge 1</math>
| |
| |-
| |
| | '''[[Fraunhofer diffraction]]''' occurs when:
| |
| <math>F = \frac{a^2}{L\lambda} \ll 1</math>
| |
| |-
| |
| |<math>a</math> – aperture or slit size,
| |
| <math>\lambda</math> – wavelength,
| |
| <math>L</math> – distance from the aperture
| |
| |}
| |
| The condition for validity is fairly weak, and it allows all length parameters to take comparable values, provided the aperture is small compared to the path length. For the ''r'' in the denominator we go one step further, and approximate it with only the first term, <math> r \approx z</math>. This is valid in particular if we are interested in the behaviour of the field only in a small area close to the origin, where the values of ''x'' and ''y'' are much smaller than ''z''. In addition, it is always valid if as well as the Fresnel condition, we have <math>L \gg \lambda</math>, where L is the distance between the aperture and the field point.
| |
| | |
| For Fresnel diffraction the electric field at point ''(x,y,z)'' is then given by:
| |
| | |
| :<math> E(x,y,z)=\frac{e^{ikz}}{i \lambda z} \iint_{-\infty}^{+\infty} E(x',y',0)e^{{ik \over 2z}[(x-x')^2+(y-y')^2]}dx'dy' </math>
| |
| | |
| This is the Fresnel diffraction integral; it means that, if the Fresnel approximation is valid, the propagating field is a spherical wave, originating at the aperture and moving along ''z''. The integral modulates the amplitude and phase of the spherical wave.
| |
| Analytical solution of this expression is still only possible in rare cases. For a further simplified case, valid only for much larger distances from the diffraction source see [[Fraunhofer diffraction]]. Unlike Fraunhofer diffraction, Fresnel diffraction accounts for the curvature of the [[wavefront]], in order to correctly calculate the relative [[Phase (waves)|phase]] of interfering waves.
| |
| | |
| ==Alternative forms==
| |
| ===Convolution===
| |
| The integral can be expressed in other ways in order to calculate it using some mathematical properties. If we define the following function:
| |
| | |
| :<math>h(x,y,z) = \frac{e^{i k z}}{i \lambda z} e^{i \frac{k}{2 z} (x^2 + y^2)}</math>
| |
| | |
| then the integral can be expressed in terms of a [[convolution]]:
| |
| | |
| :<math>E(x,y,z) = E(x,y,0) * h (x,y,z) \, </math>
| |
| | |
| in other words we are representing the propagation using a linear-filter modeling. That is why we might call the function ''h(x,y,z)'' the impulse response of free space propagation.
| |
| | |
| ===Fourier transform===
| |
| | |
| Another possible way is through the [[Fourier transform]]. If in the integral we express ''k'' in terms of the wavelength,
| |
| | |
| :<math>k = { 2\pi \over \lambda } </math> | |
| | |
| and we expand each component of the transverse displacement,
| |
| :<math>(x-x')^2 = x^2 + x'^2 -2 x x' \,</math>
| |
| :<math>(y-y')^2 = y^2 + y'^2 -2 y y' \,</math>
| |
| then we can express the integral in terms of the two dimensional Fourier transform. Let us use the following definition:
| |
|
| |
| :<math>G(p,q) = \mathcal{F} \left\{ g(x,y) \right\} \equiv \iint_{-\infty}^{\infty} g(x,y) e^{-i 2 \pi (p x + q y)} dx dy </math>
| |
| | |
| where ''p'' and ''q'' are spatial frequencies (in units of lines/meter). The Fresnel integral can be expressed as:
| |
| | |
| :<math>E(x,y,z) = \frac{e^{i k z}}{i \lambda z} e^{i \frac{\pi}{\lambda z}(x^2 + y^2)} \mathcal{F}
| |
| \left. \left\{ E(x',y',0) e^{i \frac{\pi}{\lambda z} (x'^2 + y'^2)} \right\} \right|_{p = \frac{x}{\lambda z}; q = \frac{y}{\lambda z}}
| |
| </math>
| |
| | |
| :::::<math>= H(p,q) \cdot G(p,q) |_{p = {x \over \lambda z} , q = { y \over \lambda z } }</math>
| |
| | |
| where
| |
| :<math>H(p,q) = \mathcal{F} \left\{ h(x,y) \right\} </math>
| |
| | |
| i.e. first multiply the field to be propagated for a complex exponential, calculate its two dimensional Fourier transform, replace ''(p,q)'' with <math>\left( \frac{x}{\lambda z}, \frac{y}{\lambda z} \right)</math> and multiply it by another factor. This expression is better than the others when the process leads to a known Fourier transform, and the connection with the Fourier transform is tightened in the [[linear canonical transformation]], discussed below.
| |
| | |
| === Linear canonical transformation === | |
| {{main|Linear canonical transformation}}
| |
| From the point of view of the [[linear canonical transformation]], Fresnel diffraction can be seen as a [[Shear mapping|shear]] in the [[time-frequency domain]], corresponding to how the Fourier transform is a rotation in the time-frequency domain.
| |
| | |
| ==See also==
| |
| * [[Fraunhofer diffraction]]
| |
| * [[Fresnel integral]]
| |
| * [[Fresnel zone]]
| |
| * [[Fresnel number]]
| |
| * [[Augustin-Jean Fresnel]]
| |
| * [[Ridged mirror]]
| |
| * [[Fresnel imager]]
| |
| * [[Euler spiral]]
| |
| | |
| ==Notes==
| |
| <references />
| |
| | |
| ==References==
| |
| | |
| * {{cite book |author=Goodman, Joseph W. |title=Introduction to Fourier optics |publisher=McGraw-Hill |location=New York |year=1996 |pages= |isbn=0-07-024254-2 |oclc= |doi=}}
| |
| | |
| {{DEFAULTSORT:Fresnel Diffraction}}
| |
| [[Category:Diffraction]]
| |