Chapter 36 Diffraction Masatsugu Sei Suzuki Department of Physics, SUNY at Binghamton (Date: August 15, 2020)
Total Page:16
File Type:pdf, Size:1020Kb
Chapter 36 Diffraction Masatsugu Sei Suzuki Department of Physics, SUNY at Binghamton (Date: August 15, 2020) Interference is the more general concept: it refers to the phenomenon of waves interacting. Waves will add constructively or destructively according to their phase difference. Diffraction usually refers to the spreading wave pattern from a finite-width aperture. 1. Diffraction by a single slit Fig . Single slit diffraction. Imagine the slit divided into many narrow zones, width y (= = a/N). Treat each as a secondary source of light contributing electric field amplitude E to the field at P. We consider a linear array of N coherent point oscillators, which are each identical, even to their polarization. For the moment, we consider the oscillators to have no intrinsic phase difference. The rays shown are all almost parallel, meeting at some very distant point P. If the spatial extent of the array is comparatively small, the separate wave amplitudes arriving at P will be essentially equal, having traveled nearly equal distances, that is E E ()r E (r ) ... E (r ) E ()r 0 0 1 0 2 0 N 0 N The sum of the interfering spherical wavelets yields an electric field at P, given by the real part of i( kr1t ) i( kr2 t ) i( krN t ) E Re[ E0 ()r e E0 ()r e ... E0 ()r e ] i( kr1t ) ik( r2 r 1 ) ik( r3 r 1 ) ik( rN r1 )) Re[ E0 ()r e 1[ e e ... e ]] (( Note )) When the distances r1 and r 2 from sources 1 and 2 to the field point P are large compared with the separation , then these two rays from the sources to the point P are nearly parallel . The path difference r2 – r1 is essentially equal to sin . Here we note that the phase difference between adjacent zone is a k( r r ) k sin k( sin ) 2 1 N k( r3 r2 ) k( r4 r3 ) k( rN rN 1) 2 where k is the wavenumber, k . It follows that k( r2 r1) k( r3 r1) 2 k( r4 r1) 3 (rN r1) (N )1 Thus the field at the point P may be written as i( kr1 t) i i2 i( N )1 E Re[E0()r e 1[ e e ...e ]] We now calculate the complex number given by Z 1 ei ei2 ... ei( N )1 1 eiN 1 ei eiN 2/ (eiN 2/ eiN 2/ ) ei 2/ (ei 2/ ei 2/ ) N sin( ) ei( N 1) / 2 2 sin( ) 2 If we define D to be the distance from the center of the line of oscillators to the point P, that is 1 1 kD (N )1 k sin kr (N )1 kr 2 1 2 1 1 k(D r ) (N )1 1 2 Then we have the form for E as N sin( ) ~ E Re[ E ()r ei( kD t ) 2 ] Re[ Eei t ] 0 sin( ) 2 The intensity distribution within the diffraction pattern due to N coherent, identical, distant point sources in a linear array is equal to c ~ 2 I S 0 E 2 N sin2 ( ) sin2 ( ) sin2 ( ) 2 2 2 I I0 I0 Im 2 2 2 sin ( ) sin ( ) 2 2N 2 in the limit of N→∞, where 2 sin 2 ( ) 2N 2N c I 0 [E (r)]2 0 2 0 c c I I N 2 0 [NE (r)]2 0 E 2 m 0 2 0 2 0 is the phase difference N Nk sin ka sin With k sin where a = N. We make a plot of the relative intensity I/Im as a function of . sin 2 ( ) I 2 2 I m 2 Note that 2 sin ( ) 2 d 2 2 2 IêIm 0.06 0.05 0.04 0.03 0.02 0.01 bê2p -4 -2 0 2 4 The numerator undergoes rapid fluctuations, while the denominator varies relatively slowly. The combined expression gives rise to a series of sharp principal peaks separated by small subsidiary maxima. The principal minimum occurs in directions in direction m such that ka 1 sin m , asin 2 m 2 mm 2 2 m k 2 or m m a 2. Shape of the intensity Here we examine the function form of the relative intensity sin 2 ( ) I 2 , I 2 0 ( ) 2 where 2 2 a sin a . 2 in the small limit of . The function I/I0 is equal to 4/ = 0.405 for = , and is ½ for = 2.78311 = 0.8859 . The main feature of the intensity I/I0 is that the intensity is large only in 1 1 , or sin 2 a 2 a or (sin ) a (diffraction limit) a For the diffraction pattern of a circular aperture of the diameter d, we have the condition of 22.1 d Fig. Fraunhofer diffraction pattern for a circular aperture (Mathematica). We use the FFT program of the Mathematica 3. Phasor diagram In order to discuss the intensity I, we need to calculate the sum defined by i i2 (Ni )1 Z E0 r 1)[( e e ... e ] E 0 1[ ei ei2 ... e (Ni )1 ] N in the complex plane, where k sin . For simplicity, here we solve this problem N geometrically using the phasor diagram (in the real x-y plane). There are N isosceles triangular lattices (see the Figs. below). The first one is of length E0/N and it has a phase equal to zero. The next one is of length E0/N and it has a phase equal to φ. The next one is of length E0/N and it has a phase equal to 2 φ, and soon. So we get an equiangular polygon with N sides. Fig. The resultant amplitude of N = 6 equally spaced sources with net successive phase difference φ. = N φ = 6 φ. ______________________________________________________________________ Fig. The resultant amplitude of N = 36 equally spaced sources with net successive phase difference φ. We now consider the system with a very large N. We may imagine dividing the slit into N narrow strips. In the limit of large N, there is an infinite number of infinitesimally narrow strips. Then the curve trail of phasors become an arc of a circle, with arc length equal to the length E0. The center C of this arc is found by constructing perpendiculars at O and T. The radius of arc is given by E0 E0 R R(N) , or R N in the limit of large N, where R is the side of the isosceles triangular lattice with the vertex angle φ, and the phase difference is given by N ka sin with the value being kept constant. Note that is the change of phase for two rays a separated by , N ka k sin , or sin N Then the amplitude Ep of the resultant electric field at P is equal to the chord OT , which is equal to sin E E 2R sin 2 0 sin E 2 . P 2 2 0 2 Then the intensity I for the single slits with finite width a is given by sin I I ( 2 )2 . m 2 where Im is the intensity in the straight-ahead direction where a The phase difference is given by ka sin 2 sin 2p sin . We make a plot of I/Im as a function of , where p = a/ is changed as a parameter. Fig. The relative intensity in single-slit diffraction for various values of the ratio p = a/. The wider the slit is the narrower is the central diffraction maximum. 4. Diffraction patterns 4.1 Young’s double slit experiment (two slits with finite width) We consider the Young’s double slits (the slits are separated by d). Each slit has a finite width a. Fig. Geometric construction for describing the Young’s double-slit experiment (not to scale). The intensity is given by sin 2 2 2 I I0 cos ( ) , 2 2 where 2d 2a kd sin sin ka sin sin We note that the peak due to the double slit diffraction occurs at 2 cos 1 , 2 m m 2 d The intensity becomes zero when 2 sin 0 , 2n n 2 a Fig. Diffraction pattern with the double slit with the distance d and the single slit with the distance a. d a . 1.0 0.8 0.6 0.4 0.2 -20 -10 0 10 20 1.000 0.500 0.100 0.050 0.010 0.005 0.001 -20 -10 0 10 20 Fig. Intensity ratio I/I0 vs angle (degrees) when p = d/ = 48 and q = a/ = 4. Fig. Fraunhofer diffraction pattern for the Young’s double slits. See my article on the Fraunhofer diffraction in http://physics.binghamton.edu/Sei_Suzuki/suzuki.html Fig. The intensity I as a function of y (x = 300) in the above Fraunhofer diffraction pattern. 4.2 DC SQUID Analogy of the diffraction with double slits and single slit See my article on the Josephson junction and DC SQUID for the detail; Fig. Diffraction effect of Josephson junction. A magnetic field B along the z direction, which is penetrated into the junction (in the normal phase). We consider a junction (1) of rectangular cross section with magnetic field B applied in the plane of the junction, normal to an edge of width w. q 2 J Jsin[ A d l ] , 0 1 ℏ c 1 with q = -2e. We use the vector potential A given by 1 By Bx A( B r )( , ,0) , 2 2 2 A' A (By ,0,0) , where Bxy . 2 Then we have q xb qB J J sin[ (By)dx] J sin[ yW], 0 1 ℏc 0 1 ℏc xa qB dI JLdy J Lsin[ yW ]dy , 1 0 1 ℏc Or t 2/ qB I J L sin[ yW ]dy .