Di raction: Boundary Value Problems of Electromagnetic Waves

Size: px
Start display at page:

Download "Di raction: Boundary Value Problems of Electromagnetic Waves"

Transcription

1 Chater 7 Di raction: Boundary Value Problems of Electromagnetic Waves 7. ntroduction n electrostatics, a rescribed otential distribution on a closed surface uniquely determines the otential elsewhere according to the Dirichlet s formulation, s(r )d ; G = on ; rovided G(r; r ), the Green s function, is so chosen that it vanishes on the closed surface : The normal derivative of the Green s can be regarded as a rojection or maing oerator to convert a given distribution of the surface otential s (r ) into the otential at an arbitrary observing osition r: For scalar otentials, the rojection oerator is a scalar function, or more recisely, matrix function. For any vector elds, too, rescribing a vector on a closed surface uniquely determines the vector led elsewhere. However, a vector eld changes its direction as well as its magnitude over a closed surface. A vector eld rescribed on a closed surface is to be converted into a vector eld with de nitive direction and magnitude at the observing osition. A rojection oerator that transforms one vector eld on a closed surface into another elsewhere must be a tensor (or dyadic) comosed of aroriate eigenvectors for the vector eld of concern. For electromagnetic elds, the TE and TM eigenvectors identi ed in Chater 5 can be conveniently used for this urose. n this Chater, a basic formulation will be develoed for vector boundary value roblems of electromagnetic elds, E and B: t will then be alied to some di raction roblems in which boundary electromagnetic elds are known. However, in most ractical roblems, it is di cult to know recisely boundary elds because of the feedback from the di racted elds themselves. A rigorous analysis of di raction roblems usually requires solving integral equations based on the boundary conditions for electromagnetic waves. 7. Vector Green s Theorem For arbitrary vector elds E and F, by exloiting the exansion r (E r F) = r E r F E r r F; (7.)

2 r (F r E) = r F r E F r r E; (7.) the following identity can easily be roven, (E r r F F r r E)dV = (F r E E r F) d: (7.3) ince V E r r F = E [r(r F) r F] (7.4) = r [E(r F)] (r E)(r F) E r F; (7.5) Eq. (7.3) can be modi ed as V (F r E E r F)dV = [(Fr E Er F) n (n r E) F n E r F] d; (7.6) where n is the unit normal vector on the closed surface directed away from the volume V; d = nd: (ee Fig. 7..) This is a mathematical identity and holds for arbitrary vector elds E and F. Figure 7-: For E s ; the electric eld seci ed on the closed surface ; the electric eld o the surface E is uniquely determined. n is the normal unit vector directed away from the region of interest. We now let the vector E be the electric eld and F be one comonent of Green s dyadic. n source free region, the electric eld E satis es the Helmholtz equation, (r + k )E = ; (7.7) and the Green s dyadic satis es the singular inhomogeneous Helmholtz equation (r + k )G = (r r ); (7.8) where is the unit dyadic de ned in terms of aroriate eigenvectors for the electromagnetic elds. Then,

3 Eq. (7.4) yields an exression for the electric eld E in terms of the electric eld seci ed on the closed surface ; E(r) = [(n E)r G + (n r E) G + (n E) (r G)] d; (7.9) where r E = (no charges on ) has been assumed. Note that r G is a vector and r G is a dyadic. For a diagonal Green s dyadic G = G; (7.) it follows that r G = rg; (7.) (n r E) G = (n r E)G; (7.) (n E) (r G) = (n E) rg: (7.3) Therefore, Eq. (7.9) reduces to E(r) = [(n E)rG + (n r E)G + (n E) rg] d: (7.4) The scalar Green s function G satis es and its exlicit exression in terms of sherical harmonic exansion is G(r; r ) = 4 jr r j eikjr r j ( X lx j l (kr ) h () l (kr) = ik j l (kr) h () l (kr ) l= m= l (r + k )G = (r r ); (7.5) ) Y lm (; ) Y lm ; ; ( r > r r < r ) : (7.6) As shown in Chater 6, the solution for the vector otential that satis es the Helmholtz equation r + k A = J; (7.7) can be decomosed into TE and TM modes. n the outer region r > r ; we have the following decomosition, A(r) = ik + i k X lx l= m= l X lx l= m= l l(l + ) h() l (kr)r ry lm (; ) l(l + ) r [h() l (kr)r ry lm (; )] j l (kr )r r Y lm ( ; ) J(r )dv r [j l (kr )r r Y lm ( ; )] J(r (7.8) )dv : This can be written in a more comact form using the Green s dyadic G; A(r) = e ikjr r j 4 jr r j J(r )dv = GJ(r )dv ; (7.9) 3

4 where the Green s dyadic is G(r; r ) = e ikjr r j 4 jr r j = ik X lm + i X k lm Note that the TE eigenvector and TM eigenvector l(l + ) h() l (kr)j l (kr )r ry lm (; )r r Y lm ( ; ) (7.) l(l + ) r [h() l (kr)r ry lm (; )]r [j l (kr )r r Y lm ( ; )]: (7.) h () l (kr)r ry lm (; ); (7.) r [h () l (kr)r ry lm (; )]; (7.3) are orthogonal to each other for given mode numbers l and m and the Green s dyadic is evidently diagonal. For di raction by an aerture in a large, at conducting screen, Eq. (7.4) simli es to E(r) = (n E) rgd: (7.4) aerture This formula was originally found by mythe based on an equivalent layer of magnetic diole (equivalent to a sheet of magnetic current), dm = i! n Ed; (7.5) which roduces a vector otential A (r) = i! (Recall that an oscillating magnetic diole creates a vector otential, rovided kr :) Corresonding electric eld is rg (n E)d: (7.6) A (r) = rg m ' e ikr ik m; E(r) A = rg (n E)d aerture = (n E) rgd: aerture (Exlain why we can use A (r) ; here instead of going through the usual c = r B = r (r in terms of longitudinal and transverse comonents of A:) The tangential comonent of the di racted magnetic eld in the aerture vanishes (that is, in the aerture, the tangential magnetic eld is equal to that 4

5 of the incident wave which is una ected by the aerture) and only the tangential comonent of the electric eld contributes to di raction. Likewise, the normal comonent of the di racted electric eld vanishes in the aerture. These follow from the symmetry roerty of the di racted electric and magnetic eld in front of and behind the conducting late. The di racted electric eld normal to the late is odd with resect to the normal coordinate z, E z (z) = E z ( z); (7.7) and so is the tangential magnetic eld, e z B (z) = e z B ( z) ; (7.8) but the tangential electric eld and normal magnetic eld are even, e z E (z) = e z E ( z) ; B z (z) = B z ( z): (7.9) Therefore in the aerture, the tangential comonent of the di racted magnetic eld should vanish. On the surface of the conducting late away from the aerture, tangential comonent of di racted magnetic eld can exist with a corresonding surface current n H = J s ; (A m ): ince the normal unit vector n changes sign from one side (z < ) to other (z > ) ; J s = n H is even with resect to z; J s (z = ) = J s (z = +) : The surface charges induced on both sides of the conducting late are also equal since = " n E is even with resect to z: We have encountered such symmetry roerties in Chaters and 3 when we analyzed leakage of electric and magnetic elds through a hole in a suerconducting late. n both cases, e ective dioles for elds in regions z < and z > are oosite to each other and erturbed electric and magnetic elds due to the hole have such symmetry roerties which hold in general even for radiation elds as long as the late can be regarded as ideally conducting. The mythe s formula, Eq. (7.4), may be derived urely mathematically as follows. For a at boundary, the normal unit vector n is constant, say n = e z. First, we note (n E) rgd = r (Gn E) d + Gr (n E) d = Gr (n E) d; (7.3) since for a closed surface, r [G (n E)]d = : (7.3) Recalling r E = and exloiting the fact that for a unidirectional vector n; r (n E) = (Er) n (n r) E + n (r E) E (r n) = (n r) E; (7.3) r (n E) = (Er) n+ (n r) E + n (r E) + E (r n) = (n r) E + n (r E) ; (7.33) 5

6 we nd for a constant vector n; n (r E) = r (n E) + r (n E) : (7.34) Then Gn (re) d = = Gr (n E) d + Gr(n E)d (n E) rgd (n E)rG d; and Eq. (7.4) reduces to or E(r) = [(n E)rG + (n r E)G + (n E) rg] d = (n E) rg d; whole at surface E(r) = [(n E)rG + (n r E)G] d: whole at surface n the case of di raction by an aerture in a conducting late, n E = excet in the aerture. Therefore, E(r) = (n E) rg d: (7.35) aerture This should be equal to E(r) = [(n E)rG + i!(n B)G] d; (7.36) conducting late since in the aerture where no eld discontinuity exists, both n E and nb vanish. (They are odd functions of z:) Rigorous analysis of di raction by an aerture in a conducting screen is extremely di cult mainly because in Eqs. (7.4) and (7.4), the elds on the surface and in aerture are the di racted elds which are to be found! However, in some limiting cases, sensible analytic solutions can be found. n short wavelength limit ka where a is the size of aerture, geometric otics aroximation may be used. n this case, the eld in the aerture may be aroximated by the incident eld. n long wavelength limit ka ; the eld in the aerture is totally unknown. However, di racted ower can be calculated accurately using e ective dioles. Then the eld intensity in the aerture can be estimated to be kae i E i where E i is the incident electric eld. 7.3 ome Examles of Aerture Di raction 7.3. Long lit Let a lane wave be incident on a at oaque screen having a long slit with an oening width a: The roblem is essentially two dimensional and the di racted wavevector k may be assumed to lie in the x y lane. ince the slit extends from z = to ; we are not allowed to assume r r and the di racted electric 6

7 Figure 7-: Di raction by a long slit of width a in a conducting screen. n short wavelength limit ka ; the electric eld in the slit oening may be aroximated by that of the incident wave. mythe s formula in Eq. (7.4) is alicable. eld has to be calculated from E(r) = (n E) r Gdy ; (7.37) where G is the two-dimensional Green s function satisfying (r + k )G = (r r ) = (x x )(y y ): (7.38) olution for G is G = i 4 H() hk x + (y y ) i ; (7.39) where H () (x) is the Hankel function of the rst kind. ts asymtotic form is H () (x) ' r kx eikx i 4 ; x : (7.4) Using n = n (unit normal directed into the di racted region) and noting r G = rg ; we rewrite Eq. (7.37) as E(r) = r G(n E)dy : (7.4) ubstituting H () hk i x + (y y ) ' = r k eik r k eik i 4 iky i 4 iky sin ; (7.4) 7

8 where = x + y ; (7.43) and is the angle between k and x-axis, we obtain E(r) ' a r k (n E ) k eik i sin 4 ; (7.44) with de ned by = ka sin : (7.45) The function sin = is the familiar Fraunhofer di raction form factor to describe the angular () deendence of the eld intensity. The radiation (di racted) ower er unit length of the slit is P l = = c" jej d: (7.46) = n short wavelength limit ka ; the integration limits may be aroximated by from to ; and we obtain an exected result P = c" je j a; ka : (7.47) l n the oosite limit ka ; the electric eld in the slit oening cannot be relaced by the incident eld because the current induced in the conducting late greatly disturbs the incident eld. The roblem is dual of scattering of electromagnetic waves by a thin, long conducting late which will be analyzed as a searate roblem Circular Aerture Figure 7-3: Geometry of circular aerture of radius a in an oaque screen. Radiation from a circular aerture on a large conducting late is a classical roblem and of ractical imortance. n short wavelength limit ka (a the hole radius), the eld in the aerture may be aroximated 8

9 by that of the incident wave E : The di racted electric eld is given by E(r) = i r k (n E )e ikr e ikr d; (7.48) where in the geometry shown, k r = kr sin sin : The integral can be erformed as follows: a = The di racted electric eld thus reduces to = r dr e ikr sin sin d a J (kr sin )r dr a k sin J (ka sin ): (7.49) E(r) = ia r k (n E )e ikr J (ka sin ) : (7.5) ka sin ince the amlitude of the Bessel function J (x) decreases with x in a manner J (x) _ = x; the electric eld is areciable only at small : n order of magnitude, di racted Poynting ux has an angular sread about the z axis, sin ' ' ka : The rst intensity minimum occurs at the rst root of J (x) = ; ka sin ' 3:83; or ' 3:83 ka = : ; D = a (diameter of the hole). (7.5) D This is the familiar di raction limited resolving ower of circular aertures such as mirrors in telescoes and arabolic microwave antennas. The radiation ower associated with the di racted eld is P = c" r jej d = c" je j a = J (ka sin ) sin sin d: (7.5) ince ka ; the uer limit of the integral may be extended to and sin may be aroximated by : Noting J (ka) d = ; (indeendent of ka); we nd P = c" je j a : (7.53) This is an exected result in the short wavelength limit (ka ) and simly corresonds to the incident ower going through the circular aerture. n long wavelength limit ka ; the above analysis breaks down comletely, for the electromagnetic elds in the aerture are entirely di erent from those associated with the incident wave. For normal incidence on a conducting late, the electric eld vanishes at the late but the magnetic eld is doubled because of comlete 9

10 re ection. This does not mean there exists a tangential magnetic eld of H =H in the aerture because the aerture should not a ect the incident magnetic eld and the tangential comonent of the magnetic eld in the aerture is H = H :What we can do is to nd an e ective magnetic diole moment of the hole using H as the unerturbed eld as we did in analyzing the leakage of magnetic eld through a hole in a suerconducting late. The diole moment of the hole is m = 8a3 3 H = 6a3 3 H ; (7.54) which radiates at a ower P =! 4 m 4" 3c 5 = 64 7 (ka)4 a jh j ; (7.55) where the factor is due to radiation into half sace (behind the screen). The transmission cross-section is t ' P jh j = 64 7 (ka)4 a _ k 4 a 6 : (7.56) The deendence _ k 4 (or! 4 ) is the common feature of Rayleigh scattering of electromagnetic waves by small objects (in this case an aerture). _ a 6 indicates extremely sensitive deendence of the cross section on the hole radius. To order k 6 ; the cross section is ince the radiation ower is roortional to ' 64 7 (ka)4 a + 5 (ka) + ; ka : (7.57) (n E a )e ikr d ' E aa 4 ; where E a is the electric eld in the aerture, it is evident that the eld is of the order of E a ' kae E ; being much smaller than the incident eld E : This is exected since the incident wave is essentially shortcircuited at the late and only a small eld can leak through the aerture. The electric eld in the aerture has been worked out by Bouwkam, E x (x; y) = i 4k xy 3 a x y E ; (7.58) E y (x; y) = i 4k x + y a 3 a x y E = i 4k + sin a E : (7.59) 3 a E x is an odd function of x and y and does not contribute (integrates to ). The eld diverges at the rim = a but is integrable. ntegration over the aerture yields a d de y (x; y) = i 8ka3 3 E : (7.6)

11 The di racted electric eld can then be found using the mythe s formula which is equivalent to the eld due to an e ective magnetic diole moment in Eq. (7.54). Di raction in the regime ka ' O() is di cult to analyze and requires numerical analysis. The transmission cross-section eaks at ka ' :56 and its value is max ' :8 a. This behavior is similar to the case of scattering by a conducting shere Radiation from an Oen End of a Rectangular Waveguide Figure 7-4: Radiation from an oen end of a rectangular waveguide. The electric eld in the oening is E y (x) = E sin(x=a) and the tangential comonent of the magnetic eld is H x (x) = E(x)= where is the imedance of the TE mode. The electric eld of the TE mode in a rectangular waveguide is given by where a is the width of the waveguide and E(x; z; t) = E sin a x e i(kzz!t) e y ; (7.6) k z =! c (!c =!) ; (7.6) is the axial wavenumber with! c = c a ; (7.63) being the cuto frequency of the TE mode. The accomanying magnetic eld transverse to the axis is H(x; z; t) = E sin a x e i(kzz!t) e x ; (7.64) where = r " (!c =!) ; (7.65) is the imedance of the TE mode. (The axial magnetic eld H z associated with the TE mode does not

12 contribute to radiation from the aerture.) f a waveguide is truncated, the oen radiates electromagnetic waves. The mythe s formula is not alicable because it only ertains to radiation from aertures in a large conducting late. n the general di raction formula where E(r) = we assume r r. Then, and Also, Then, E(r) ' ieikr 4r (n E)r G + G(n r E) + (n E) r G d; (7.66) G(r; r ) = eikjr r j 4 jr r j ; G(r; r ) ' 4r eikr ikr ; r G = rg ' ik 4r eikr r = i!b: ikr : k (n E) +!n B k(n E) e ikr d; (7.67) where n = n is the normal unit directed toward the region where the electric eld is to be evaluated. t is not obvious that the contribution from the second and third terms in the integrand,!n B k(n E) e ikr d; (7.68) is exlicitly transverse (erendicular to k) which should be the case since we are calculating radiation eld. To rove this, we rst show that r (n B) G (r; r ) d; can be exressed in terms of the normal comonent of the electric eld as follows. r G (n B) d = r G (n B) d = r G (B n) d = r G B nd = r (GB) nd = + i! c G (E n) d; Gr B nd where use is made of the identity r A nd = :

13 Therefore, r r G (n B) d = rr = i! c r i! = c G (n B) d + k G (n B) d G (E n) d + k G (n B) d (E n) r Gd + k G (n B) d; and the di racted eelctric eld can alternatively be calculated from E(r) = = (n E)r G + G(n r E) + (n E) r G d rg (n E) + i r r (Gn H) d:!" Then Eq. (7.67) may be rewritten as E(r) = ieikr 4r k n E k (n H) e ikr d; (7.69)!" which is exlicitly transverse to k: As mentioned in Chater 5, this formula was rst derived by chelkuno in term of a ctitious magnetic current to relace the tangential comonent of the electric eld on a boundary. The rice to be aid for introducing a magnetic current is that the Maxwell s equation, r B = ; has to be violated. n the derivation resented here, no magnetic current is assumed. Denoting the surface integral over the oen end of the waveguide by (; ) = = a dx b dy sin a x e ik(x sin cos +y sin sin ) i =a k sin sin (=a) k sin cos ( + e ika sin cos )( e ikb sin sin ); (7.7) we nally nd the electric eld radiated from the oen end of a rectangular waveguide, E(r) = " ikeikr (; ) 4r cos (!c =!) +! sin e + Calculation of the radiation ower is left for an exercise.! # (!c =!) + cos cos e : (7.7) cattering by a Conducting here Revisited This roblem has been solved rigorously in the receding Chater in terms of sherical harmonic exansion. However, the result is not very illuminating hysically articularly in the short wavelength limit ka wherein summation of a large number of functions involving sherical Bessel functions and their derivatives must be erformed. n this limit, geometrical otics aroximation should be able to yield the scattered electric eld rovided the boundary conditions for the electromagnetic elds are aroriately incororated. The incident lane wave sees the cross-section of the shere a which corresonds to the scattering crosssection due to scattering by the illuminated hemisherical surface facing the incident wave. cattering by the illuminated surface is nothing but re ection by a sherical convex mirror. The Poynting ux associated with re ection is uniform (isotroic) and indeendent of the angles and. The total correct cross-section 3

14 Figure 7-5: cattering by a conducting shere in short wavelength limit ka : as revealed from the rigorous analysis was a : The additional a is due to shadow scattering from the hemisherical surface in the shadow of the incident wave where the eld vanish. This vanishing eld can be interreted as cancellation between the eld of the incident wave and that of the scattered wave. hadow scattering is dual of di raction by a circular aerture in a conducting screen which radiates a ower P sh = c" je j a ; ka : (7.7) Although the scattered ower through shadow scattering is identical to that by the illuminated surface, the Poynting ux is sharly eaked behind the shere as in the case of di raction with an angular sread of order ' =a : f the shere is ideally conducting, the boundary conditions for the electric and magnetic elds are E t = ; B n = : (7.73) Let the electric eld of the incident wave be E i and that of scattered eld be E s : The vanishing tangential comonent of the total eld requires that The vanishing normal comonent of the magnetic eld requires that n the sirit of geometrical otics valid in the short wavelength regime, n (E i +E s ) = : (7.74) n (B i +B s ) = : (7.75) E i = c k k ib i ; E s = c k k sb s ; 4

15 where k i is the wavevector of the incident wave and k s is that of the scattered wave. From n (E i +E s ) = n (k i B i +k s B s ) = ; (7.76) and it follows that Likewise, from we nd n k i = n k s ; n k i = n k s ; (7.77) n (B i B s ) = : (7.78) n (k i E i k s E s ) = ; n (E i E s ) = : With these rearations, we now aly Eq. (7.67) to the illuminated and shadow surfaces searately. The contribution from the illuminated surface is E ill (r) = ieikr 4r (k k ) (n E i ) + (k k )(n E i ) e ikr d; where the incident magnetic eld has been eliminated through!b i = k E i : n the geometry shown in Fig. 7-5, the incident wave roagates in the negative z-direction, E i (r) = E e ikz = E e ikr cos : (7.79) Therefore, the electric eld di racted by the illuminated surface is E ill (r) = ieikr 4r (k k ) (n E ) + (k k )(n E ) e ikr cos ikr d: (7.8) imilarly, the contribution from the shadow surface is E sh (r) = ieikr 4r (k + k ) (n E ) (k k )(n E ) e ikr cos ikr d: (7.8) The total di racted eld is given by the sum of E ill and E sh : However, the angular deendence of the eld intensities is entirely di erent as exlained in the introduction, namely, je ill j is insensitive to and ; while je sh j sharly eaked in the direction = (forward scattering). Therefore, the scattered ower can be calculated searately as if the elds were incoherent. Noting r = a on the shere surface, we nd the hase function reduces to kr cos + k r = ka ( + cos ) cos + sin sin cos( ) = ka f(; ; ; ): (7.8) 5

16 By assumtion ka : Therefore, the exonential function e ikaf ; (7.83) raidly oscillates as and are varied. The slow angular deendence in the amlitude function (k k ) (n E ) + (k k )(n E ) can be ignored in the integration over and and can be taken out of the integral, E ill (r) ' i eikr (k k ) (n E i ) + (k k 4r )(n E i ) a ; e ikaf( ; ) sin d d ; where and indicate the angular location on the shere surface which makes the dominant contribution to the hase integral e ikaf( ; ) sin d d : (7.84) ince ka ; major contribution to the integral comes from the angular location where the function f( ; ) = : This determines = ; = : This is a trivial result well exected from otical re ection. At an observing angular location (; ); only the wave re ected at ( = =; = ) can be detected. Let us Taylor exand the hase function f( ; ) about the stationary hase oint = = = and = = ; f(; ; ; ) ' cos! sin ( ) + : (7.86) ntegration over yields " = ex ika cos " ' sin ex ika cos = sin # sin d # d ika cos(=) : (7.87) Note that the integration limits can be extended to because only the region ' = contributes to the 6

17 integral. imilarly, ' ' ex ika cos sin ( ) d ex ika cos sin ( ) d (7.88) ika cos(=) sin(=) : (7.89) Therefore, e ikaf( ; ) d = a i ex ( ika cos(=)) : (7.9) ka cos(=) At the articular angular location (=; ); the magnitude of the vector jk k cos and the normal vector n is in the direction of k k ; n = k j is ; (7.9) ; ; : (7.9) Then, (k k ) (n E ) + (k k )(n E ) = k cos [n (n E ) + n (n E )] = k cos [n (n E ) E ] : (7.93) The vector has a magnitude of E ; that is, n (n E ) E ; (7.94) jn (n E ) E j = E ; and is erendicular to the scattered wavevector k. (Prove this statement.) Therefore, the electric eld scattered (or re ected) by the illuminated surface is E ill (r) = a r eikr ika cos(=) [n (n E ) E ] : (7.95) The amlitude je ill (r)j is indeendent of the angles and as exected from isotroic re ection and scattered ower can readily be found as P ill = c" Er a 4 r = c" E a : (7.96) 7

18 The scattering cross-section due to scattering (essentially re ection) by the illuminated surface is therefore ill = a : (7.97) cattering by the shadow surface is identical to di raction by a circular aerture on a conducting screen analyzed earlier, for the boundary conditions E i + E s = ; B i + B s = ; (7.98) indicate that the elds of the incident wave can be used on the surface of the shadow hemishere. (The sign inversion does not a ect the di racted ower.) The shae of the shadow surface is irrelevant as long as its area rojected normal to the incident wave is circular with an area a : To see this, let us consider the electric eld scattered by the shadow surface in Eq. (7.8), E sh (r) = ieikr 4r (k + k ) (n E ) (k k )(n E ) e i(k k)r d : (7.99) ince there is no stationary hase oint in the shadow (it only occurs in the illuminated surface) and the function e i(k k)r raidly oscillates because of the assumtion ka ; the integral will be nonvanishing only if k k ' which means redominantly forward scattering. The change in the wavevector k = k k is thus erendicular to k and k ; k = k sin e = k sin (cos e x + sin e y ); (7.) where is olar angle of scattered wavevector k now measured form the negative z axis, e is the radial unit vector in the cylindrical coordinates and is the azimuthal angle about the z axis. The hase function becomes (k k) r = ak sin sin cos ; (7.) and the eld amlitude may be aroximated by which is indeendent of and. Then (k + k ) (n E ) (k k )(n E ) ' k cos E ; (7.) (k + k ) (n E ) (k k )(n E ) e i(k k)r d ' ke a sin d d e iak sin sin cos( ) cos a = ke d = k (e z E ) d de ik sin cos( ) a d e ik sin cos( ) ; (7.3) where = a sin ; n is the unit vector normal to the circular at surface, and k E ' k E = is 8

19 noted. The electric eld scattered by the sherical surface in the shadow is thus given by E sh (r) = ieikr r k (e z E ) a d e ik sin cos( ) d: (7.4) This is the negative of the eld di racted by a circular aerture given in Eq. (7.48). ntegrations over and yield E sh (r) = ia r k (e z E ) J (ka sin ) k sin ' ia r k (e z E ) J (ka) ; : (7.5) k The radiation ower associated with shadow di raction is identical to that in the case of circular aerture, Therefore, the total ower reradiated by the shere is and the total scattering cross-section in the short wavelength limit is P sh = c" E a : (7.6) P = P ill + P sh = c" E a ; (7.7) = a ; ka : (7.8) This is consistent with the asymtotic (ka! ) value of the general cross-section derived earlier, where (ka) = k A l = X h (l + ) ja l j + jb l j i ; (7.9) l= j l (ka) h () l (ka) ; B l = d da [aj l(ka)] d (ka)] da [ah() l : (7.) The summation over l in Eq. (7.9) requires u to l ' O(ka) for su cient accuracy and thus becomes cumbersome in otical regime where ka can be huge. Also, it should be noted that the assumtion of ideally conducting shere (" e = ; e = ) entirely breaks down in high frequency (short wavelength) regime. The e ective relative ermittivity of ordinary metals in otical frequency regime is in general comlex and remains of the order of unity. (The relative magnetic ermeability in the regime may be assumed to be unity.) 7.4 cattering by a Knife Edge n scattering by a knife edge, there is no geometrical scale size to seak of excet for k where is the radial distance from the edge ti. n short wavelength otical regime, assuming an incident electric eld along the lane above the edge may rovide the lowest order aroximation. n rigorous aroach, the boundary conditions for the electromagnetic elds must be incororated. We assume an observing oint at (; ): The 9

20 Figure 7-6: Geometry for analyzing scattering by a conducting knife edge. incidence. Bottom: general incidence angle. To: secial case of normal distance between a oint on the y axis and the observing oint is = ( sin y ) + cos ' y sin + y : (7.) Therefore, for small ; the ath di erence between and is y ; and corresonding hase di erence is The electric eld at (x ; y) is thus given by ky : (7.) The function, shown in Fig. (??), y E(y) = E q k F (x) = x is known by Fresnel s integral and its secial values are: i ex t dt: (7.3) i ex t dt; (7.4) F ( ) = ; F () = ( + i); F () = + i:

21 At large ositive y (well above the edge); the incident eld E is recovered as exected. Well below the edge in the shadow region, y! ; the eld vanishes also as exected. At y = ; the eld is E and intensity is 4 of the incident intensity. Fig. (7-8) shows the intensity as a function of the normalized vertical distance q y! y x : F (x) = x i ex t dt Figure 7-7: Fresnel integral F (x) = R x ex i t dt: F ( ) = ; F () = + i: y 4 6 q Figure 7-8: ntensity = as a function of normalized vertical distance y x. (y = ) = 4 ; max ' :37 : The analysis above is rather rimitive because no consideration was given to the boundary conditions of the electric and magnetic elds. For incident electric eld arallel to the conducting late, E = E e z (E-mode), the total eld (incident + scattered) must vanish at the surface of the conductor late. f the

22 incident wave aroaches the late at an angle ; the total wave eld is given by E z (; ) = i k= cos E e ik cos( ) e i t dt i E e ik cos(+ ) k= sin + e i t dt: (7.5) f the incident wave is so olarized that its magnetic eld is arallel to the late, H = H e z (H-mode), the boundary condition is This yields E = at = ; 3 : H z (; ) = i k= cos H e ik cos( ) e i t dt + i H e ik cos(+ ) k= sin + e i t dt: (7.6) The reader should verify that the solutions ( rst formulated by ommerfeld) satisfy the Helmholtz equation and the resective boundary conditions.

23 Problems 7. An electric eld E(; ) = E (; )e + E (; )e is seci ed on the surface of a shere of radius a: Determine the radiation electric eld. 7. A magnetic eld B(; ) = B (; )e + B (; )e is seci ed on the surface of a shere of radius a: Determine the radiation magnetic eld. 7.3 A lane wave is incident normal to a long conducting cylinder having a square cross-section of side a:the electric eld is axial (E z only) and the incident wave falls normal to one of the rectangular surfaces. Derive an integral equation for the surface current density J z (r): (The integral equation can be solved numerically following the rocedure develoed for nding the caacitance of a square conducting late.) 7.4 Reeat the receding roblem for H mode (H z only). 7.5 A oint light source is laced at a distance z on the axis of an oaque disk of radius a ( ): Determine the light intensity along the axis behind the disk as a function of distance z: 7.6 A oint light source is laced at a distance z on the axis of a circular oening of radius a in an oaque screen. Determine the light intensity behind the disk as a function of axial distance z: 7.7 how that in Fresnel di raction of incident wave normal to the late, the maximum intensity occurs at y ' : x = and is given by max ' :37 where is the intensity of the incident wave. 7.8 how that near a knife edge, an E-mode axial electric eld and corresonding magnetic eld H = E z (; ) = A sin ; A cos i! e x + sin e y ; satisfy the boundary conditions. The x comonent of the magnetic eld is discontinuous at = and = : nterret this eculiarity. how that the intensity of a wave scattered by a knife edge is insensitive to the observing angle : (A knife edge aears shiny regardless of observing angle.) 7.9 The transmission cross section of a small (ka ) circular aerture in a conducting late is T = 64 7 k4 a 6 ; where a is the hole radius. how that the magnitude of the electric eld in the aerture should be of the order of E a ' kae ; where E is the incident eld. Note: For incident eld olarized in the x direction, Bouwkam found the electric eld comonent in the aerture resonsible for di raction, E y (x; y) = 8k x + y a 3 a x y E : 3

24 (E x (x; y) is an odd function of both x and y:) ntegration over the aerture yields E y (x; y) d = which is consistent with the e ective magnetic diole, E y (r; ) rdrd = 6 3 ka3 E ; m = 8 3 a3 H : At the edge of a thin conductor, the electric eld arallel to the edge vanishes in a manner E k = ; while the normal comonent diverges as E? = = where is the distance from the edge. For examle, the electric eld at the rim of a charged conducting disk behaves as E? _ q = : a (a ) a uch roerties can be exloited to nd electric eld when the boundary involves shar conducting edge.) 7. n low frequency limit, scatterimg by a conducting shere may be analyzed using diole aroximation. Relevant diole moments are = 4" a 3 E ; m = a 3 H ; where E and H are the elds of the incident lane wave. how that the di erential scattering cross section is " d d = k4 a 6 cos sin + cos cos # ; the total scattering cross section is and that the force exerted on the shere is = 3 k4 a 6 ; F z = 4 3 k4 a 6 " E ; (N). 4

INTERFERENCE, DIFFRACTION AND SCATTERING OF ELECTROMAGNETIC WAVES

INTERFERENCE, DIFFRACTION AND SCATTERING OF ELECTROMAGNETIC WAVES Chapter 13 INTERFERENCE, DIFFRACTION AND SCATTERING OF ELECTROMAGNETIC WAVES 13.1 Introduction Waves having the same frequency (thus the same wavelength) can interfere with each other. Interference is

More information

CHAPTER 25. Answer to Checkpoint Questions

CHAPTER 25. Answer to Checkpoint Questions CHAPTER 5 ELECTRIC POTENTIAL 68 CHAPTER 5 Answer to Checkoint Questions. (a) negative; (b) increase. (a) ositive; (b) higher 3. (a) rightward; (b),, 3, 5: ositive; 4: negative; (c) 3, then,, and 5 tie,

More information

On the relationship between sound intensity and wave impedance

On the relationship between sound intensity and wave impedance Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi

More information

Focal Waveform of a Prolate-Spheroidal IRA

Focal Waveform of a Prolate-Spheroidal IRA Sensor and Simulation Notes Note 59 February 6 Focal Waveform of a Prolate-Sheroidal IRA Carl E. Baum University of New Mexico Deartment of Electrical and Comuter Engineering Albuquerque New Mexico 873

More information

The individual electric and magnetic waves are in phase. The fields peak at the same position at the same time.

The individual electric and magnetic waves are in phase. The fields peak at the same position at the same time. 1 Part 3: Otics 3.1: Electromagnetic Waves An electromagnetic wave (light wave) consists of oscillating electric and magnetic fields. The directions of the electric and magnetic fields are erendicular.

More information

A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAPE

A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAPE '. ' 220 HILlS TECHNICAL REVIEW VOLUME 20 A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAE 621.317.331:538.632.083 Resistivity and Hall-coefficient measurements

More information

Chapter 6: Sound Wave Equation

Chapter 6: Sound Wave Equation Lecture notes on OPAC0- ntroduction to Acoustics Dr. Eser OLĞAR, 08 Chater 6: Sound Wave Equation. Sound Waves in a medium the wave equation Just like the eriodic motion of the simle harmonic oscillator,

More information

Highly improved convergence of the coupled-wave method for TM polarization

Highly improved convergence of the coupled-wave method for TM polarization . Lalanne and G. M. Morris Vol. 13, No. 4/Aril 1996/J. Ot. Soc. Am. A 779 Highly imroved convergence of the couled-wave method for TM olarization hilie Lalanne Institut d Otique Théorique et Aliquée, Centre

More information

Chater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de

Chater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms

More information

E o e associated with a light field (both the real part and the. ikr t. under the assumptions that J free

E o e associated with a light field (both the real part and the. ikr t. under the assumptions that J free Reiew Problems Chaters 1-5 True and False Questions E1. T or F: The otical index of any material aries with frequency. E2. T or F: The frequency of light can change as it enters a crystal. E3. T or F:

More information

Meshless Methods for Scientific Computing Final Project

Meshless Methods for Scientific Computing Final Project Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking

More information

Introduction to Group Theory Note 1

Introduction to Group Theory Note 1 Introduction to Grou Theory Note July 7, 009 Contents INTRODUCTION. Examles OF Symmetry Grous in Physics................................. ELEMENT OF GROUP THEORY. De nition of Grou................................................

More information

References: 1. Cohen Tannoudji Chapter 5 2. Quantum Chemistry Chapter 3

References: 1. Cohen Tannoudji Chapter 5 2. Quantum Chemistry Chapter 3 Lecture #6 Today s Program:. Harmonic oscillator imortance. Quantum mechanical harmonic oscillations of ethylene molecule 3. Harmonic oscillator quantum mechanical general treatment 4. Angular momentum,

More information

Wave Drift Force in a Two-Layer Fluid of Finite Depth

Wave Drift Force in a Two-Layer Fluid of Finite Depth Wave Drift Force in a Two-Layer Fluid of Finite Deth Masashi Kashiwagi Research Institute for Alied Mechanics, Kyushu University, Jaan Abstract Based on the momentum and energy conservation rinciles, a

More information

Focused azimuthally polarized vector beam and spatial magnetic resolution below the diffraction limit

Focused azimuthally polarized vector beam and spatial magnetic resolution below the diffraction limit Research Article Vol. 33, No. 11 / November 2016 / Journal of the Otical Society of America B 2265 Focused azimuthally olarized vector beam and satial magnetic resolution below the diffraction limit MEHDI

More information

NONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA)

NONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA) NONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA) Note: SFA will automatically be taken to mean Coulomb gauge (relativistic or non-diole) or VG (nonrelativistic, diole-aroximation). If LG is intended (rarely),

More information

pp physics, RWTH, WS 2003/04, T.Hebbeker

pp physics, RWTH, WS 2003/04, T.Hebbeker 1. PP TH 03/04 Accelerators and Detectors 1 hysics, RWTH, WS 2003/04, T.Hebbeker 2003-12-03 1. Accelerators and Detectors In the following, we concentrate on the three machines SPS, Tevatron and LHC with

More information

On a Markov Game with Incomplete Information

On a Markov Game with Incomplete Information On a Markov Game with Incomlete Information Johannes Hörner, Dinah Rosenberg y, Eilon Solan z and Nicolas Vieille x{ January 24, 26 Abstract We consider an examle of a Markov game with lack of information

More information

Principles of Computed Tomography (CT)

Principles of Computed Tomography (CT) Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel

More information

we get our formula for the saddle-point integral (5.294) h(w) e xf(w). (5.301) xf 00 (w) (w j )

we get our formula for the saddle-point integral (5.294) h(w) e xf(w). (5.301) xf 00 (w) (w j ) 5.3 The Abel-Plana Formula and the Casimir E ect 9 and using f (w) = e i and + = to show that e i = e i i = e i = f (w) (5.3) we get our formula for the saddle-oint integral (5.94) I(x) / h(w) e xf(w).

More information

Statics and dynamics: some elementary concepts

Statics and dynamics: some elementary concepts 1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and

More information

Positivity, local smoothing and Harnack inequalities for very fast diffusion equations

Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract

More information

LECTURE 3 BASIC QUANTUM THEORY

LECTURE 3 BASIC QUANTUM THEORY LECTURE 3 BASIC QUANTUM THEORY Matter waves and the wave function In 194 De Broglie roosed that all matter has a wavelength and exhibits wave like behavior. He roosed that the wavelength of a article of

More information

Radiation properties of leaky modes near the spectral gap region for semi-infinite printed-circuit lines

Radiation properties of leaky modes near the spectral gap region for semi-infinite printed-circuit lines RADIO SCIENCE, VOL. 38, NO. 3, 151, doi:1.129/22rs2777, 23 Radiation roerties of leaky modes near the sectral ga region for semi-infinite rinted-circuit lines Frank J. Villegas, 1 David R. Jackson, and

More information

REFLECTION AND TRANSMISSION BAND STRUCTURES OF A ONE-DIMENSIONAL PERIODIC SYSTEM IN THE PRESENCE OF ABSORPTION

REFLECTION AND TRANSMISSION BAND STRUCTURES OF A ONE-DIMENSIONAL PERIODIC SYSTEM IN THE PRESENCE OF ABSORPTION Armenian Journal of Physics, 0, vol. 4, issue,. 90-0 REFLECTIO AD TRASMISSIO BAD STRUCTURES OF A OE-DIMESIOAL PERIODIC SYSTEM I THE PRESECE OF ABSORPTIO A. Zh. Khachatrian State Engineering University

More information

CMSC 425: Lecture 4 Geometry and Geometric Programming

CMSC 425: Lecture 4 Geometry and Geometric Programming CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas

More information

Spin Diffusion and Relaxation in a Nonuniform Magnetic Field.

Spin Diffusion and Relaxation in a Nonuniform Magnetic Field. Sin Diffusion and Relaxation in a Nonuniform Magnetic Field. G.P. Berman, B. M. Chernobrod, V.N. Gorshkov, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 V.I. Tsifrinovich

More information

Chapter 5. Transient Conduction. Islamic Azad University

Chapter 5. Transient Conduction. Islamic Azad University Chater 5 Transient Conduction Islamic Azad University Karaj Branch 1 Transient Conduction Many heat transfer roblems are time deendent Changes in oerating conditions in a system cause temerature variation

More information

Higher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface

Higher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface Higher order theory for analytic saddle oint aroximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface P.F Daley ABSTACT The high frequency solution to the roblem of a Ρ P and Ρ S reflected

More information

ZEEMAN EFFECT: p...(1). Eigenfunction for this Hamiltonian is specified by

ZEEMAN EFFECT: p...(1). Eigenfunction for this Hamiltonian is specified by ZEEMAN EFFECT: Zeeman Effect is a magneto-otical henomenon discovered by Zeeman in 1896. He observed that when an atom (light soce) is laced in an external magnetic field, the sectral lines it emits are

More information

8.7 Associated and Non-associated Flow Rules

8.7 Associated and Non-associated Flow Rules 8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more

More information

Complex Analysis Homework 1

Complex Analysis Homework 1 Comlex Analysis Homework 1 Steve Clanton Sarah Crimi January 27, 2009 Problem Claim. If two integers can be exressed as the sum of two squares, then so can their roduct. Proof. Call the two squares that

More information

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)

More information

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system

More information

Potential Theory JWR. Monday September 17, 2001, 5:00 PM

Potential Theory JWR. Monday September 17, 2001, 5:00 PM Potential Theory JWR Monday Setember 17, 2001, 5:00 PM Theorem 1 (Green s identity). Let be a bounded oen region in R n with smooth boundary and u, v : R be smooth functions. Then ( ( ) u v v u dv = u

More information

Averaging sums of powers of integers and Faulhaber polynomials

Averaging sums of powers of integers and Faulhaber polynomials Annales Mathematicae et Informaticae 42 (20. 0 htt://ami.ektf.hu Averaging sums of owers of integers and Faulhaber olynomials José Luis Cereceda a a Distrito Telefónica Madrid Sain jl.cereceda@movistar.es

More information

The oerators a and a obey the commutation relation Proof: [a a ] = (7) aa ; a a = ((q~ i~)(q~ ; i~) ; ( q~ ; i~)(q~ i~)) = i ( ~q~ ; q~~) = (8) As a s

The oerators a and a obey the commutation relation Proof: [a a ] = (7) aa ; a a = ((q~ i~)(q~ ; i~) ; ( q~ ; i~)(q~ i~)) = i ( ~q~ ; q~~) = (8) As a s They can b e used to exress q, and H as follows: 8.54: Many-body henomena in condensed matter and atomic hysics Last modied: Setember 4, 3 Lecture. Coherent States. We start the course with the discussion

More information

Do Gravitational Waves Exist?

Do Gravitational Waves Exist? Universidad Central de Venezuela From the electedworks of Jorge A Franco etember, 8 Do Gravitational Waves Exist? Jorge A Franco, Universidad Central de Venezuela Available at: htts://works.beress.com/jorge_franco/13/

More information

Estimating Time-Series Models

Estimating Time-Series Models Estimating ime-series Models he Box-Jenkins methodology for tting a model to a scalar time series fx t g consists of ve stes:. Decide on the order of di erencing d that is needed to roduce a stationary

More information

Quantitative estimates of propagation of chaos for stochastic systems with W 1, kernels

Quantitative estimates of propagation of chaos for stochastic systems with W 1, kernels oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract

More information

Feedback-error control

Feedback-error control Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller

More information

P043 Anisotropic 2.5D - 3C Finite-difference Modeling

P043 Anisotropic 2.5D - 3C Finite-difference Modeling P04 Anisotroic.5D - C Finite-difference Modeling A. Kostyukevych* (esseral echnologies Inc.), N. Marmalevskyi (Ukrainian State Geological Prosecting Institute), Y. Roganov (Ukrainian State Geological Prosecting

More information

u = 1 (B 2 + E2 E B (16.2) + N = j E (16.3) One might be tempted to put u and N into a 4-vector N and write the equation in the form

u = 1 (B 2 + E2 E B (16.2) + N = j E (16.3) One might be tempted to put u and N into a 4-vector N and write the equation in the form Chater 6 Energy-momentum tensor (Version., 3 November 7 Earlier, we obtained for the energy density and flux u = (B + E µ c (6. We also had a continuity equation N = µ E B (6. u t + N = j E (6.3 One might

More information

Andrea Mantile. Fractional Integral Equations and Applications to Point Interaction Models in Quantum Mechanics TESI DI DOTTORATO DI RICERCA

Andrea Mantile. Fractional Integral Equations and Applications to Point Interaction Models in Quantum Mechanics TESI DI DOTTORATO DI RICERCA DOTTORATO DI RICERCA in MATEMATICA APPLICATA E INFORMATICA Ciclo XVI Consorzio tra Università di Catania, Università di Naoli Federico II, Seconda Università di Naoli, Università di Palermo, Università

More information

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites

More information

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS

PHYSICS 44 {CLASSICAL MECHANICS 2003 SUMMER FINAL EXAM { SOLUTIONS PHYSICS 44 {CLASSICAL MECHANICS 3 SUMMER FINAL EXAM { SOLUTIONS Problem 1 Consider the following Atwood machine comosed of the two masses m 1 and m and attached to a stand with the hel of a sring of strength

More information

Topological-phase effects and path-dependent interference in microwave structures with magnetic-dipolar-mode ferrite particles

Topological-phase effects and path-dependent interference in microwave structures with magnetic-dipolar-mode ferrite particles Toological-hase effects and ath-deendent interference in microwave structures with magnetic-diolar-mode ferrite articles Abstract M. Berezin, E.O. Kamenetskii, and R. Shavit Microwave Magnetic Laboratory

More information

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density Adv. Studies Theor. Phys., Vol. 7, 13, no. 1, 549 554 HIKARI Ltd, www.m-hikari.com A Note on Massless Quantum Free Scalar Fields with Negative Energy Density M. A. Grado-Caffaro and M. Grado-Caffaro Scientific

More information

Numerical Linear Algebra

Numerical Linear Algebra Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and

More information

Dispersion relation of surface plasmon wave propagating along a curved metal-dielectric interface

Dispersion relation of surface plasmon wave propagating along a curved metal-dielectric interface Disersion relation of surface lasmon wave roagating along a curved metal-dielectric interface Jiunn-Woei Liaw * and Po-Tsang Wu Deartment of Mechanical Engineering, Chang Gung University 59 Wen-Hwa 1 st

More information

CHAPTER 2: SMOOTH MAPS. 1. Introduction In this chapter we introduce smooth maps between manifolds, and some important

CHAPTER 2: SMOOTH MAPS. 1. Introduction In this chapter we introduce smooth maps between manifolds, and some important CHAPTER 2: SMOOTH MAPS DAVID GLICKENSTEIN 1. Introduction In this chater we introduce smooth mas between manifolds, and some imortant concets. De nition 1. A function f : M! R k is a smooth function if

More information

Chapter 3. GMM: Selected Topics

Chapter 3. GMM: Selected Topics Chater 3. GMM: Selected oics Contents Otimal Instruments. he issue of interest..............................2 Otimal Instruments under the i:i:d: assumtion..............2. he basic result............................2.2

More information

CR extensions with a classical Several Complex Variables point of view. August Peter Brådalen Sonne Master s Thesis, Spring 2018

CR extensions with a classical Several Complex Variables point of view. August Peter Brådalen Sonne Master s Thesis, Spring 2018 CR extensions with a classical Several Comlex Variables oint of view August Peter Brådalen Sonne Master s Thesis, Sring 2018 This master s thesis is submitted under the master s rogramme Mathematics, with

More information

Equivalence of Wilson actions

Equivalence of Wilson actions Prog. Theor. Ex. Phys. 05, 03B0 7 ages DOI: 0.093/te/tv30 Equivalence of Wilson actions Physics Deartment, Kobe University, Kobe 657-850, Jaan E-mail: hsonoda@kobe-u.ac.j Received June 6, 05; Revised August

More information

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty

More information

Session 5: Review of Classical Astrodynamics

Session 5: Review of Classical Astrodynamics Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve

More information

Radiation Torque Exerted on a Spheroid: Analytical Solution

Radiation Torque Exerted on a Spheroid: Analytical Solution Cleveland State University EngagedScholarshi@CSU Physics Faculty Publications Physics Deartment 7-1-8 Radiation Torque Exerted on a Sheroid: Analytical Solution Feng Xu James A. Lock Cleveland State University,

More information

Elements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley

Elements of Asymptotic Theory. James L. Powell Department of Economics University of California, Berkeley Elements of Asymtotic Theory James L. Powell Deartment of Economics University of California, Berkeley Objectives of Asymtotic Theory While exact results are available for, say, the distribution of the

More information

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.

More information

ε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T

ε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T High High Power Power Laser Laser Programme Programme Theory Theory and Comutation and Asects of electron acoustic wave hysics in laser backscatter N J Sircombe, T D Arber Deartment of Physics, University

More information

The directivity of the forced radiation of sound from panels and openings including the shadow zone

The directivity of the forced radiation of sound from panels and openings including the shadow zone The directivity of the forced radiation of sound from anels and oenings including the shadow zone J. Davy RMIT University, Alied Physics, GPO Box 476V, 3001 Melbourne, Victoria, Australia john.davy@rmit.edu.au

More information

Torsion. Click here to check the Animation. Introduction. Basic Assumptions. Assumptions. Torsion Formula. Stress Formula. Stresses on Inclined Planes

Torsion. Click here to check the Animation. Introduction. Basic Assumptions. Assumptions. Torsion Formula. Stress Formula. Stresses on Inclined Planes Torsion ntroduction Basic Assumtions Torsion Formula Stresses on nclined Planes Angle of Twist in Torsion Torsion of Circular Elastic Bars: Formulae Click here to check the Animation Assumtions Stress

More information

Dimensional perturbation theory for Regge poles

Dimensional perturbation theory for Regge poles Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West

More information

Elementary Analysis in Q p

Elementary Analysis in Q p Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some

More information

Casimir Force Between the Two Moving Conductive Plates.

Casimir Force Between the Two Moving Conductive Plates. Casimir Force Between the Two Moving Conductive Plates. Jaroslav Hynecek 1 Isetex, Inc., 95 Pama Drive, Allen, TX 751 ABSTRACT This article resents the derivation of the Casimir force for the two moving

More information

Central Force Motion Challenge Problems

Central Force Motion Challenge Problems Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has

More information

Chapter 7 Rational and Irrational Numbers

Chapter 7 Rational and Irrational Numbers Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers

More information

Positive decomposition of transfer functions with multiple poles

Positive decomposition of transfer functions with multiple poles Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.

More information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. Score normalization technical details We now discuss the technical details of the score normalization method. SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules

More information

SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY

SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in

More information

Lecture 1.2 Units, Dimensions, Estimations 1. Units To measure a quantity in physics means to compare it with a standard. Since there are many

Lecture 1.2 Units, Dimensions, Estimations 1. Units To measure a quantity in physics means to compare it with a standard. Since there are many Lecture. Units, Dimensions, Estimations. Units To measure a quantity in hysics means to comare it with a standard. Since there are many different quantities in nature, it should be many standards for those

More information

2016-r1 Physics 220: Worksheet 02 Name

2016-r1 Physics 220: Worksheet 02 Name 06-r Physics 0: Worksheet 0 Name Concets: Electric Field, lines of force, charge density, diole moment, electric diole () An equilateral triangle with each side of length 0.0 m has identical charges of

More information

Entropic forces in dilute colloidal systems

Entropic forces in dilute colloidal systems Entroic forces in dilute colloidal systems R. Castañeda-Priego Instituto de Física, Universidad de Guanajuato, Lomas del Bosque 103, Col. Lomas del Camestre, 37150 León, Guanajuato, Mexico A. Rodríguez-Lóez

More information

Partially Quenched Chiral Perturbation Theory and the Replica Method

Partially Quenched Chiral Perturbation Theory and the Replica Method Partially Quenched Chiral Perturbation Theory and the Relica Method P. H. Damgaard and K. Slittorff The Niels Bohr Institute Blegdamsvej 7 DK-200 Coenhagen Ø Denmark March 24, 2000 Abstract We describe

More information

Multiparameter entanglement in quantum interferometry

Multiparameter entanglement in quantum interferometry PHYSICAL REVIEW A, 66, 023822 200 Multiarameter entanglement in quantum interferometry Mete Atatüre, 1 Giovanni Di Giusee, 2 Matthew D. Shaw, 2 Alexander V. Sergienko, 1,2 Bahaa E. A. Saleh, 2 and Malvin

More information

A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands

A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands A Secial Case Solution to the Persective -Point Problem William J. Wolfe California State University Channel Islands william.wolfe@csuci.edu Abstract In this aer we address a secial case of the ersective

More information

Section 4: Electromagnetic Waves 2

Section 4: Electromagnetic Waves 2 Frequency deendence and dielectric constant Section 4: Electromagnetic Waves We now consider frequency deendence of electromagnetic waves roagating in a dielectric medium. As efore we suose that the medium

More information

On split sample and randomized confidence intervals for binomial proportions

On split sample and randomized confidence intervals for binomial proportions On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have

More information

integral invariant relations is not limited to one or two such

integral invariant relations is not limited to one or two such The Astronomical Journal, 126:3138 3142, 2003 December # 2003. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT ORBIT INTEGRATION BY SCALING AND ROTATION FOR CONSISTENCY

More information

Design and Correction of Optical Systems

Design and Correction of Optical Systems Design and Correction of Otical Systems Lecture 5: Wave aberrations 017-05-19 Herbert Gross Summer term 017 www.ia.uni-jena.de Preliminary Schedule - DCS 017 1 07.04. Basics 1.04. Materials and Comonents

More information

Phase transition. Asaf Pe er Background

Phase transition. Asaf Pe er Background Phase transition Asaf Pe er 1 November 18, 2013 1. Background A hase is a region of sace, throughout which all hysical roerties (density, magnetization, etc.) of a material (or thermodynamic system) are

More information

#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS

#A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS #A37 INTEGERS 15 (2015) NOTE ON A RESULT OF CHUNG ON WEIL TYPE SUMS Norbert Hegyvári ELTE TTK, Eötvös University, Institute of Mathematics, Budaest, Hungary hegyvari@elte.hu François Hennecart Université

More information

Convergence performance of the coupled-wave and the differential methods for thin gratings

Convergence performance of the coupled-wave and the differential methods for thin gratings Convergence erformance of the couled-wave and the differential methods for thin gratings Philie Lalanne To cite this version: Philie Lalanne. Convergence erformance of the couled-wave and the differential

More information

A QUATERNIONIC APPROACH to GEOMETRY of CURVES on SPACES of CONSTANT CURVATURE

A QUATERNIONIC APPROACH to GEOMETRY of CURVES on SPACES of CONSTANT CURVATURE INTERNATIONAL JOURNAL OF GEOMETRY Vol. 3 (2014), No. 1, 53-65 A QUATERNIONIC APPROACH to GEOMETRY of CURVES on SPACES of CONSTANT CURVATURE TUNA BAYRAKDAR and A. A. ERG IN Abstract. We construct the Frenet-Serret

More information

CHAPTER 3: TANGENT SPACE

CHAPTER 3: TANGENT SPACE CHAPTER 3: TANGENT SPACE DAVID GLICKENSTEIN 1. Tangent sace We shall de ne the tangent sace in several ways. We rst try gluing them together. We know vectors in a Euclidean sace require a baseoint x 2

More information

The role of current loop in harmonic generation from magnetic metamaterials in two polarizations

The role of current loop in harmonic generation from magnetic metamaterials in two polarizations The role of current loo in harmonic generation from magnetic metamaterials in two olarizations Iman Sajedian 1,2, Inki Kim 2, Abdolnasser Zakery 1 and Junsuk Rho 2,3* 1 Deartment of Physics, College of

More information

THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN

THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN INTERNATIONAL JOURNAL OF GEOMETRY Vol. 5 (2016), No. 1, 31-38 THE ERDÖS - MORDELL THEOREM IN THE EXTERIOR DOMAIN PETER WALKER Abstract. We show that in the Erd½os-Mordell theorem, the art of the region

More information

MATH 2710: NOTES FOR ANALYSIS

MATH 2710: NOTES FOR ANALYSIS MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite

More information

Airy Functions. Chapter Introduction

Airy Functions. Chapter Introduction Chater 4 Airy Functions 4. Introduction Airy functions are named after the English astronomer George Biddell Airy (8 892). Airy s first mathematical work was on the diffraction henomenon, namely, the Airy

More information

Beam-Beam Stability in Electron-Positron Storage Rings

Beam-Beam Stability in Electron-Positron Storage Rings Beam-Beam Stability in Electron-Positron Storage Rings Bjoern S. Schmekel, Joseh T. Rogers Cornell University, Deartment of Physics, Ithaca, New York 4853, USA Abstract At the interaction oint of a storage

More information

Spin as Dynamic Variable or Why Parity is Broken

Spin as Dynamic Variable or Why Parity is Broken Sin as Dynamic Variable or Why Parity is Broken G. N. Golub golubgn@meta.ua There suggested a modification of the Dirac electron theory, eliminating its mathematical incomleteness. The modified Dirac electron,

More information

Minimal Surfaces in R 3

Minimal Surfaces in R 3 Minimal Suaces in R 3 AJ Vargas May 10, 018 1 Introduction The goal of this note is to give a brief introduction to the theory of minimal suaces in R 3, and to show how one would go about generalizing

More information

A Social Welfare Optimal Sequential Allocation Procedure

A Social Welfare Optimal Sequential Allocation Procedure A Social Welfare Otimal Sequential Allocation Procedure Thomas Kalinowsi Universität Rostoc, Germany Nina Narodytsa and Toby Walsh NICTA and UNSW, Australia May 2, 201 Abstract We consider a simle sequential

More information

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS #A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,

More information

The diagonal Born±Oppenheimer correction to molecular dynamical properties

The diagonal Born±Oppenheimer correction to molecular dynamical properties 6 January Chemical Physics Letters 333 () 459±464 www.elsevier.nl/locate/clett The diagonal Born±Oenheimer correction to molecular dynamical roerties Sohya Garashchuk a, *, John C. Light a, Vitaly A. Rassolov

More information

Waveguide Coupler I. Class: Integrated Photonic Devices Time: Fri. 8:00am ~ 11:00am. Classroom: 資電 206 Lecturer: Prof. 李明昌 (Ming-Chang Lee)

Waveguide Coupler I. Class: Integrated Photonic Devices Time: Fri. 8:00am ~ 11:00am. Classroom: 資電 206 Lecturer: Prof. 李明昌 (Ming-Chang Lee) Waveguide Couler I Class: Integrated Photonic Devices Time: Fri. 8:am ~ 11:am. Classroom: 資電 6 Lecturer: Prof. 李明昌 (Ming-Chang Lee) Waveguide Couler n 1 > n n Waveguide 1 n 1 n Waveguide n 1 n How to switch

More information

Phase velocity and group velocity (c) Zhengqing Yun,

Phase velocity and group velocity (c) Zhengqing Yun, Phase velocity and grou velocity (c) Zhengqing Yun, 2011-2012 Objective: Observe the difference between hase and grou velocity; understand that the grou velocity can be less than, equal to, and greater

More information

Very Very Important Questions along with their Solutions for Board Exams. Physics

Very Very Important Questions along with their Solutions for Board Exams. Physics 1 Very Very Imortant Questions along with their Solutions for Board Exams. Physics SECTION A (Each Question carries 1 Mark) 1. The lot of the variation of otential difference across a combination of three

More information

Convex Optimization methods for Computing Channel Capacity

Convex Optimization methods for Computing Channel Capacity Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem

More information

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &

More information