Transverse angular shift in the reflection of light beams

Size: px
Start display at page:

Download "Transverse angular shift in the reflection of light beams"

Transcription

1 1 August 000 Optcs Communcatons Transverse angular shft n the reflecton of lght beams Javer Alda ) Optcs Department. UnÕersty Complutense of Madrd, School of Optcs, AÕda. Arcos de Jalon srn, 8037 Madrd, Span Receved 5 January 000; receved n revsed form 6 Aprl 000; accepted 4 May 000 Abstract The reflecton of a lght beam onto a planar nterface can be out of the ncdence plane due to the msmatchng between the symmetres of the beam and the reflecton geometry. The state of polarzaton of the beam s also nvolved n the descrpton of the transverse angular shft. The moment characterzaton of laser beams s used, along wth the power expanson of the reflecton coeffcents and the geometrc factors, to fnd ths transverse angular shft. q 000 Elsever Scence B.V. All rghts reserved. 1. Introducton The reflecton of laser beams at delectrc and metallc nterfaces s a very common phenomenon n practcal systems handlng and transformng laser beams. The analyss of a delectrcrdelectrc nterw1,x restrcted face has been deeply treated by Porras to the plane of ncdence. He uses the characterstc parameters of a lght beam defned n terms of the moments of the beam. In ths paper we apply a smlar procedure to the calculaton of the transverse angular shft. Ths effect means a reflecton out of the plane of ncdence. The transverse angular shft was outlned and demonstrated n the 70 s for total nternal reflecton w3 5 x. Hugonn and Pett also appled a moment calculaton to obtan the dsplacement of the center of gravty of the beam. A few years ago, Nasalsk presented a characterzaton of the non-specular effects n the reflecton of lght beams wx 6. Hs treatment s based on a detaled ) Tel.: q ; fax: q ; e-mal: j.alda@fs.ucm.es analyss of the changes produced n the reflecton. Some of the results of our contrbuton wll be compared wth Nasalsk s calculaton. To characterze the beam we calculate the moments of the ntensty dstrbuton of ts plane wave decomposton. The transformaton of the moments by the reflecton s calculated by usng a Taylor expanson of the reflecton coeffcents and by takng nto account the geometry of the ncdence and the state of polarzaton of the ncomng beam. The reflected ampltude for a delectrcrmetallc plane nterface s obtaned usng the complex nature of the reflecton coeffcents. Once the moments of the reflected beam are obtaned, the angular shfts along the plane of ncdence and along a perpendcular drecton are calculated. The results show that an out-of-plane reflecton s obtaned when the symmetry of the beam does not match wth the symmetry of the reflecton geometry. It also depends on the state of polarzaton of the beam. Along the paper we wll assume that the beam s paraxal. Ths means that t extends wthn a small regon around the prncpal drecton of propagaton r00r$ - see front matter q 000 Elsever Scence B.V. All rghts reserved. PII: S

2 ( ) J. AldarOptcs Communcatons In Secton we frst descrbe how the three dmensonal beam s transformed by a plane nterface. We calculate the transformaton by decomposng the beam nto plane waves and by usng the reflecton coeffcents. A detaled analyss of the geometry of the ncdence wll be necessary to correctly apply the reflecton coeffcents. Also the state of polarzaton wll play an mportant role n the calculaton. Secton 3 s devoted to the change of the moments characterzng the beam. A Taylor expanson untl second order s appled around the axs defned by the center of the ncomng beam. In Secton 4 we calculate the transverse angular shft wth the general equaton obtaned to transform the moments. Some cases are treated numercally n ths secton to show up the nfluence of the symmetry of the beam, ts orentaton, ts modal character, and the state of polarzaton. Fnally, Secton 5 summarzes the man conclusons of the paper.. Transformaton of the ampltude dstrbuton A gven laser beam can be descrbed by ts electrc feld. The vectoral character s relevant to fnd ts components along the ncdence and perpendcular planes. If we assume that the laser propagates along the Z-axs, the electrc feld at a gven z plane s gven by ts ampltude dstrbuton n that plane, C Ž x, y., and by ts polarzaton state. In ths paper we have dropped the tme dependence because we are prmarly nterested n the spatal and angular characterzaton. The ampltude dstrbuton s expanded as a sum of plane waves whose ampltudes are gven by the Fourer transform of the ampltude feld dstrbuton, FŽ j,h.. Ths angular spectral dstrbuton, the angular spectrum, also carres nformaton of the polarzaton state of the beam. If the beam propagates around ts man drecton nsde a paraxal angular regon, the magntude of the plane wave decomposton wthn a homogeneous delectrc s expressed as follows, FŽ j,h, z. ½ 5 p 1 sfž j,h. exp z 1y l Ž j qh., l Ž 1. where l s the wavelength wthn the homogeneous delectrc meda, and j and h are the spatal frequences along the X and Y drectons respectvely. The drecton of propagaton for each plane wave s gven by the untary vector ˆ Is lj,lh,( 1y Ž lj qlh.. Ž. ž / The nterface s characterzed by the ndex of refracton n both sdes of the nterface and by ts untary normal vector N, ˆ ˆNs Ž ysn 0,0,ycos 0., Ž 3. where 0 s the angle of ncdence of the center of the ncomng beam. In the reflecton, the ampltude of each plane wave changes accordng wth the value of the reflecton coeffcent, r. Dependng on the polarzaton state of the ncomng wavefront, t wll be necessary to use both the parallel and the perpendcular reflecton coeffcents, or ether one. These coeffcents are gven by the Fresnel formulae n terms of the ndexes of refracton and the angle of ncdence of the ncomng plane wave wx 7. The ncdence angle for each plane wave s gven as: Ž j,h. y1 scos Ž NPI ˆ ˆ.. To obtan the reflected plane wave for each ncomng component t wll necessary: Ž. to apply the vectoral reflecton law to know the drecton of the reflected plane wave, and Ž. to multply the ncdent ampltude dstrbuton by the reflecton coeffcent to obtan the ampltude of the reflected plane wave. The vectoral form of the reflecton law wx 8 provdes the followng expresson for the reflected vector: ˆ Rs ljcos ysn ( 1ylj ylh,lh, ž 0 0 ( ylj sn 0ycos 0 1ylj ylh. 4 Once the drecton of the reflecton s known by usng the prevous equaton, t s possble to transform the ncomng plane wave nto the reflected one by changng the coordnate system and modfyng the ampltude by means of the reflecton coeffcents. The reference system for the ncdent beam was chosen n such a way that Eq. Ž. 1 descrbes the beam Ž see Fg. 1.. Ths means that the Z axs represents the straght lne where the center of gravty of the /

3 ( ) J. AldarOptcs Communcatons The electrc feld vector for any ncomng component of the plane wave decomposton can be wrtten as, F Ž j,h. sf Ž j,h. uˆ qf Ž j,h. uˆ x x y y sf Ž j,h. uˆ qf Ž j,h. u ˆ. Ž H H Fg. 1. Geometry of the coordnate systems for the descrpton of the ncdent and reflected beams. The coordnate system for the reflecton s the mrror mage of the coordnate system for the ncdent beam. beam les. Then, the frst order moments of the ntensty dstrbuton of the angular spectrum are zero for the ncdent beam. Besdes, the reference system s algned wth the plane of ncdence. The XZ plane defnes the man plane of ncdence,.e., the plane of ncdence for the center of the beam. Ths choce means that the projecton on the XZ plane of the electrc feld of the plane wave that propagates wthn ths plane s the parallel component for the center of the beam. The reflected beam s descrbed wth a reference system that s the mrror reflecton of the one used for the ncdent beam. Therefore, the Z X axs remans close to the man drecton of propagaton of the beam, and X X and Y X represent the coordnates on the transversal plane. The actual drecton of the propagaton of the reflected beam s gven by the frst order moments of ts plane-wave decomposton. The departure of ths actual drecton from the Z X axs defnes the angular shfts. The ncdent and the reflected reference systems have opposte characters: one of them s dextro and the other s levo. However, ths can be taken nto account f necessary. The relaton between the ncdent and reflected reference systems makes the coordnate transformaton very X smple: Ž x, y X. sž x, y., and n the spatal-frequency X doman Ž j,h X. sž j,h.. In ths equaton uˆ5 and uˆh are two untary vectors along the parallel and perpendcular drectons respectvely. On the other hand, uˆx and uˆy are the untary vectors along the axes of reference of the transversal plane XY. In general uˆx and uˆ5 do not concde. Therefore, the decomposton n parallel and perpendcular drecton becomes a key ssue that s analyzed as follow. For every plane wave of the angular spectrum t s possble to defne a local ncdence plane. Ths plane s defned by the ncdent and normal vectors ŽEqs. Ž. and Ž 3... The cross product of these two vectors defnes a vector perpendcular to the local ncdence plane. Ths vector s gven as zsi=n. ˆ ˆ The projecton of ths vector on the transversal plane, XY, defnes the drecton of the perpendcular component of the feld. Ths drecton s gven by the Õx and Õy components of the vector z that can be wrtten as, Õx sylhcos 0, Ž 6. ( y 0 0 Õ sljcos y 1y lj y lh sn. 7 By usng these two components we defne an angle, b, that represents the local rotaton of the parallel and perpendcular drectons wth respect to the man parallel and perpendcular drectons. These man drectons concde wth the X and Y drectons respectvely. The angle wll be gven as bsytan y1 Õ rõ. Ž 8. x y The reflected angular spectrum can be calculated by locally rotatng the ncdent angular spectrum to meet the actual parallel and perpendcular drectons. After ths rotaton we apply the reflecton coeffcents, and fnally the rotaton s undone to retreve

4 4 ( ) J. AldarOptcs Communcatons the correct orentaton. Ths transformaton can be wrtten n matrcal form as follows, Fx r cosb ysnb r5 0 cosb snb Fx r s ž / ž / ž 0 r / ž / y H Fy ž / F snb cosb ysnb cosb r cos b q r sn b Ž r y r. snbcosb 5 H 5 H Fx s. ž ž Ž r y r. snbcosb r sn b q r cos b / Fy / 5 H 5 H Ž 9. The last = matrx of the prevous equaton can be denoted as the reflecton matrx. Ths result concdes wth Eq. Ž 61. of Ref. wx 5. As we stated at the begnnng of ths paper, ths analyss s vald for lght beams whose angular extenson s wthn the paraxal range. Ths condton s also related wth the transversalty condton for lght beams. Ths means that once the paraxal condton can be appled then the component of the electrc feld along the propagaton drecton at a gven pont of a curved wavefront can be neglected. Therefore, wthn the paraxal approach, the electrc vector les on the transversal plane XY. 3. Transformaton of the moments The moments of the ntensty dstrbuton of the angular spectrum, m, j, are defned by the followng ntegral, HH m s < FŽ j,h. < j j h k dj dh. Ž 10. j,k The angular man drectons of propagaton and the components of the dvergence tensor are gven by w9,10x ux slm1,0rm 0,0, Ž 11. uy slm0,1rm 0,0, Ž 1. xx,0 0,0 Ž 1,0 0,0. u s 4l m rm y m rm, 13 uxys4l m1,1rm0,0y m1,0 m0,1 r m 0,0, Ž 14. yy 0, 0,0 Ž 0,1 0,0. u s 4l m rm y m rm, 15 both for the ncdent and reflected beam, as a functon of the moments of the angular spectrum of the beam. A superscrpt or r wll be added to dstngush between the ncdent and the reflected beam respectvely. For the ncdent beam, the electrc feld s gven by Eq. Ž. 5 wth respect to an orthogonal reference system on the transversal plane. At ths pont t s necessary to ntroduce a parameter descrbng the state of polarzaton of the ncomng wave. Ths parameter s the followng rato Fx Ž j,h. xž j,h. s. Ž 16. F Ž j,h. y Nasalsk defnes t n terms of the parallel and perpendcular drectons. However, n our case t s more nterestng to use the rato between the x and y components. Ths parameter s complex n nature < < a and wll be expressed as x s x e. Wthn the paraxal approach, as Nasalsk states, the polarzaton rato can be consdered constant for the whole beam, xž j,h. fx 0. If ths condton s not fulflled then the dependence of x wth respect to Ž j,h. must be taken nto account. By usng the rato of polarzaton the squared modulus of the ncdent plane wave s calculated as < F < s < x < q 1 < F <, Ž 17. y and therefore there s a relaton between the moments of the total ntensty and the moments of the ntensty of the component along the drecton Y. Ths relaton s, m, y j,k mj,k s, 18 < x < 0 q1 where the superndex, y means the moment of the component y of the ncdent beam. We assume that x0 s constant on the transversal plane. After some calculaton usng Eq. Ž. 9 the square modulus of the reflected dstrbuton, < F r <, can be wrtten as, r 4 4 ½ 5 H < F < s < r < cos bq< r < sn b ž / q < r < q< r < sn bcos b < x < 5 H 0 < < 4 < < 4 q r sn bq r cos b 5 H ž 5 H / ž 5 H / q < r < q< r < sn bcos b < r < y< r < =snbcosb< x < cosa < F <, Ž 19. where a s the phase of x y

5 ( ) J. AldarOptcs Communcatons As we can see, the result nvolves not only the reflecton coeffcents Ž r, r. 5 H, but also the geometry of the ncdence Ž b., and the polarzaton state of the beam Ž x. 0. These dependences can be obtaned as a functon of the spatal frequences Ž j,h.. In the paraxal regme the values of the reflecton coeffcents for the whole beam are around the values of the coeffcent for the center of the beam,.e., jsh s0, that ncdes at an angle 0. Then, after expandng n powers, we wrte the reflecton coeffcent as, r Ž j,h. sr qr jqr hq r j qr jh 0 qr h q PPP, Ž 0. where the superscrpts ndcate the order of the dervaton wth respect to the varables Ž j,h..be- sdes, the dervatves must be evaluated at jshs0. In the prevous equaton we dropped the subscrpts 5 and H because the expanson s done for both coeffcents. After some calculus t s found that, wthn our reference system, the frst dervatve of r wth respect to h s zero, both for the parallel and the perpendcular drectons. Ths fact s a consequence of the geometry of the problem. Gven an ncdence angle, 0, f we take another ncdence drecton out of the man plane of ncdence Žplane XZ., the reflecton coeffcent wll vary n a symmetrc way wth respect to the value obtaned nsde the man ncdence plane. Therefore, the dependence of r wth respect to h reaches a statonary pont at hs0. Then, wthn the second order approach gven n the prevous equaton, there are only four coeffcents dfferent from zero. These coeffcents are: r 00, r 10, r 0, and r 0. Ths expanson s also referred n Nasalsk s paper n Eqs. Ž 7.a. and Ž 7.b.. However, he uses an exponental form for the reflecton coeffcents that s more adapted to hs formalsm. The stuaton of the local plane of ncdence and the decomposton onto the parallel and perpendcular local drectons are characterzed by the angle b. Ths angle was defned n the prevous secton ŽEq. Ž.. 8. The trgonometrc functons that appear n Eq. Ž 19. can be also expanded untl second order as follows, cos 4 bs1y L h, Ž 1. sn 4 bsl h, Ž. cosbsnbsylh, Ž 3. cos bsn bsl h, Ž 4. where L s lrtan 0. These equatons also provde another restrcton to the expanson. The angle of ncdence must be large enough to allow the approach. For example, f the ncdence s normal, the angle b vares from yprtopr nterchangng the role of the parallel and perpendcular components along the way. In ths case the approxmaton clearly fals. Therefore, we must assure that the transverse extenson of the beam s small compared wth the actual value of the ncdence angle. The transverse extenson of the ncdent beam s characterzed by the parameter uyy. Then, ths condton can be wrt- ten as uyy < 0. Ž 5. Also the polarzaton rato could be expanded n powers of the spatal frequences f necessary. At ths pont, and by usng these prevous equatons t s possble to calculate the moments of the reflected beam n terms of the moments of the ncdent beam. After applyng the prevous expansons for r 5, r H, and the trgonometrc functons, the moments of the reflected beam are as follows 1 N r mj,ks Ý a m, Ž 6 k,l jql,kqm. < x < 0 q1 lqms0 where the coeffcents are gven by < 00 < < < < 00 a s r x q r <, Ž 7. 0,0 5 0 H Ž H H H H. ž / ) 00) 10 a s r r qr r < x < 1, q r 00 r 10 ) qr 00 ) r 10, Ž 8. < 00 < < 00 a sy r y r < Lx < < cosa, Ž 9. 0,1 5 H 0 a s r r qr r q< r < < x < ) 00) 0 10, q r r qr r q< r <, Ž ) 00) 0 10 H H H H H ) 00) 10 1,1 Ž a s r r qr r Ž H H H H. 0 ž ) 00) 10 y r r qr r Lx < < cosa, Ž 31. a s r r qr r ql < r < ) 00) , H / ž < < < < 5 H / y< r < < x < q r r qr r ) 00) H H H H ql r y r. 3

6 6 ( ) J. AldarOptcs Communcatons The angular characterstcs of the reflected beam can be obtaned from the results of Eq. Ž 6.. These angular parameters were defned n Eqs. Ž 11. Ž Transverse angular shft One of the most nterestng consequences of the calculaton presented prevously s that t predcts a reflecton outsde of the plane of ncdence. Transwx 6 verse angular shfts have been already descrbed and demonstrated wx 3 by analyzng the three dmensonal propagaton of the electromagnetc feld along the nterface of two delectrc meda w11 13 x. In the present paper ths result s obtaned by calculatng the characterstc parameters of the reflected beam n three dmensons, extendng the Porras analyss who calculated nsde the plane of ncdence w1, x. Based on the transformaton of plane waves by means of the reflecton coeffcents, our approach allows to descrbe both the delectrcrdelectrc nterface partal reflecton, and the delectrcrmetallc nterface reflecton. The transverse angular shft s gven by Eq. Ž 1.. To better calculate t, we wll assume that the ncdent beam s algned wth the Z axs. Ths means that m1,0 sm 0,1 s0, and therefore ux suy s0. After some substtutons usng Eqs. Ž 7. Ž 3. and Eqs. r 13 15, we fnd the followng results for uy wthn the zero, frst, and second order of approxmaton n the calculaton of the moments of the reflected beam Ž Ns0,1, and respectvely.. u r Ž 0. su s0, Ž 33. y y 1 r uy Ž 1. s Ž a1,0uxyqa0,1u yy., Ž 34. 4l uy r Ž. 0 uxy uyy m,1 m1, m 0,3 a1,0 q a0,1 q a,0 q a1,1 q a 0, 4l 4l m0,0 m0,0 m0,0 sl. a q a uxx q a uxy q a uyy 0,0,0 1,1 0, 4l 4l 4l Ž 35. Fg.. Plot of the a coeffcents, normalzed to the factor Ž< x < q1., j 0, as a functon of the ncdence angle 0. The plots are for the case of an ArrAl plane nterface. Each plot contans several cases of lnear or ellptc polarzatons.

7 ( ) J. AldarOptcs Communcatons The zeroth order approach produces a null transverse shft. For the frst order approach we already fnd a non null transverse shft that depends on the elements of the dvergence tensor. The dependence nvolves hgher order moments Ž untl 3rd order. f the approach s extended to the nd order. Eq. Ž 34. can be compared wth Eq. Ž 40. of the Nasalsk s paper. Although Nasalsk splts the calculaton nto the parallel and perpendcular components, the dependence wth the angular dvergence s also found when the Raylegh range dstance used n hs calculatons s wrtten n terms of the dvergence. The second order approach ŽEq. Ž 35.. nvolves the same knd of dependence but n a lttle more complcated fashon. In the present paper t s possble to fnd the transverse angular shft wthout the parallel and perpendcular decomposton. Actually, the geometry of the ncdence and the state of polarzaton are ncluded n the a coeffcents through the Taylor, j expanson used to obtan the coeffcents n terms of L Žsee Eqs. Ž 1. Ž 4.. and the rato of polarzaton, x Numercal examples In ths secton we perform some numercal calculaton of the transverse angular shft n the second order approxmaton for some few typcal examples. These examples nvolve crcular and ellptc Gaussan beams, Hermte Gaussan modes, and several states of polarzaton. The calculaton wavelength s ls633 nm, and the nterfaces are ArrAl Žn Al s 1.44q5.3., ArrAg Ž n s0.q3.44. Ag, and external and nternal reflecton n an ArrGlass Žn s 1.5. nterface. The transformaton of the moments depends on the value of the a, j coeffcents. They are related wth the values of the refracton ndexes nvolved n Fg. 3. Varaton of the longtudnal and transverse angular shft as a functon of the angle of ncdence and the ellptcty of a Gaussan beam for the parallel and perpendcular lnear polarzatons. The axes of the ellptcal ntensty dstrbuton are orented at 458 wth respect to the plane of ncdence. The value of the orthogonal Gaussan dvergences along the axes of the ellpse are: u1 s 0 mrad, and u s 0, 13, and 10 mrad as ndcated n the fgure. The drecton of the axs descrbed by u1 s along the frst-thrd quadrants of the XY plane. The nterface s between two delectrc meda of ndex 1.5 and 1 Ž nternal reflecton..

8 8 ( ) J. AldarOptcs Communcatons the nterface, n, n t, the angle of ncdence, 0, and the state of polarzaton of the ncomng beam defned by x 0. Fg. represents the values of the aj,k coeffcents ŽEqs. Ž 7. Ž 3.. as a functon of the angle of ncdence, 0, for an ArrAl plane nterface. Actually, the plots nclude the normalzaton factor 1rŽ< x < q 1. 0 that allows the comparson between dfferent states of polarzatons. Although the numercal values of these coeffcents are small, they are multpled by the values of the moments of the ncdent beam to obtan the moments of the reflected beam. The moments are calculated n the spatalfrequency doman and they can reach very large values producng a non neglgble angular shft. As we wll see n the cases presented below, the numercal value of the transverse angular shft s, n certan cases, comparable to the longtudnal angular shft obtaned wthn the plane of ncdence. The symmetry of the ncdent beam s a crtcal characterstc to enhance the transverse angular shft. Ths fact s shown n Fg. 3 where we plotted the angular shfts for ellptc Gaussan beams wth dfferent rato between axes. The ellptcal ntensty dstrbutons are orented havng ther axs 458 apart from the plane of ncdence. The case s an nternal reflecton for a delectrcrdelectrc Ž n rn s 1r1.5. nterface. The beams are lnearly polarzed along the parallel and perpendcular man drectons. These drectons concde wth the X and Y axes respectvely. The plots on the left correspond wth the longtudnal angular shft, u r. The plots on the rght x are for the transverse angular shft, uy r. We can see that u r s about twce larger than u r. The transverse x angular shft ncreases wth the ellptcty of the ncomng beam, beng neglgble wth respect to ux r for crcularly symmetrc beams. Our results for ux r y t Fg. 4. Varaton of the longtudnal and transverse angular shfts for several Gaussan Hermte modes. The orentaton of the axes of the ellpse s 458 wth respect to the plane of ncdence Ž as n Fg. 3.. The Gaussan dvergences of the mode Ž 0,0. are u1 s 0 mrad and u s 6.7 mrad. The modes represented are Ž 1,1., Ž,., and Ž 3,3., as shown n the fgure. The states of polarzaton are lnear along the parallel and perpendcular man drectons as ndcated n the fgures. The nterface s ArrAg.

9 ( ) J. AldarOptcs Communcatons concde very well wth those presented n Fg. of Ref. wx 1 takng nto account that n our case the dvergence s 0 mrad. The sgns also change due to the opposte orentaton of the X axs. Beyond the lmt angle and wthn our formalsm the a, j coeffcents untl second order are: a0,0s < x < q1, a sa sa sa s0, and 0 1,0 0,1,0 1,1 a s < x < q 1 L. Ž 36. 0, 0 Therefore, after some substtutons, the transverse angular shft for total nternal reflecton s gven by m 0,3rm r 3 0,0 uy Ž. s4l. Ž 37. 4tan qu 0 yy Ths equaton means that f the beam s symmetrc wth respect to the XZ plane then m 0,3 s0, and there s not transverse angular shft. Therefore the contrbuton of the asymmetry of the beam to the effect would be reported by ths equaton. The beams used n our examples are all symmetrc wth respect to the XZ. Therefore, the calculated transverse angular shft s zero after the lmt angle Ž see Fg. 3.. The multmodal character of the beam also affects the value of the transverse angular shft. In Fg. 4 we plotted the value of uy r calculated for several Her- mte Gaussan modes ncdng onto an ArrAg nterface. The fgure shows up that hgher modes mples hgher angular shfts. Both ux r and uy r are plotted n ths fgure for lnear polarzed Hermte Gaussan beams wth ellptc shape orented at 458 and orthogonal Gaussan dvergences of 0 mrad and 6.7 mrad. The state of polarzaton s lnear pontng along the parallel and the perpendcular drectons. In Fg. 5 we have calculated the longtudnal and transverse angular shft for a Gaussan beam wth crcular shape and dvergence of 0 mrad for an Fg. 5. Varaton of the longtudnal and transverse angular shfts as a functon of the angle of ncdence for a crcularly symmetrc Gaussan beam havng a dvergence of 0 mrad. For the fgures on the top the polarzaton s lnear orented at dfferent angles: 08 Ž parallel., 458, and 90 8 Ž perpendcular.. At the bottom, the fgures corresponds wth ellptcal polarzaton states havng < x < 0 s1. The nterface s ArrGlass Ž ns1.5..

10 10 ( ) J. AldarOptcs Communcatons ArrGlass Ž n s 1.5. ncdence. The plots n the top of the fgure correspond wth lnear polarzaton at 08 Ž parallel man drecton along the X axs., 458, and 908 Ž perpendcular man drecton along the Y axs.. In these cases the transverse angular shft s about two orders of magntude smaller than ux r. The bot- tom of the fgure represents several cases of ellptc polarzed lght. We fx < x < 0 s 1. Now the longtudnal angular shft s about ten tmes smaller than for the lnear polarzaton case, and t does not depend on the phase of x Conclusons In ths paper we have proposed a method for calculatng the angular moments of reflected threedmensonal laser beams wthn the paraxal approach. Ths method can be appled to delectrcrdelectrc and delectrcrmetallc reflectons. Ths study completes prevous contrbutons restrcted to the plane of ncdence. The method uses a power expanson of the reflecton coeffcents and apples t to the three-dmensonal plane-wave decomposton of the ncdent beam. The geometry of the ncdence and the local change of the plane of ncdence s ncorporated to the calculaton. Its dependence wth respect to j and h s also expanded untl second order. The polarzaton state of the beam s assumed to be constant wthn the paraxal approach. As a result of the method, we focus on the transverse angular shft. Ths effect has been descrbed prevously wthn a more pure electromagnetc framework. The same effect s obtaned here wthn the moment characterzaton of laser beams. The actual value of the transverse shft depends on the matchng of the symmetry of the beam and the geometry of the ncdence. It also depends on the multmodal character of the beam and the state of polarzaton of the ncomng beam. Acknowledgements Ths work was prmarly developed durng a stay of the author at the Gnzton Laboratory of the Stanford Unversty. The author s deeply grateful to the program Becas del Amo of the Unversty Complutense of Madrd, for fnancng the stay as a Vstng Scholar. I also want to thank Prof. Anthony Segman for hs support durng the stay. References wx 1 M.A. Porras, Opt. Commun. 131 Ž wx M.A. Porras, Opt. Commun. 135 Ž wx 3 C. Imbert, Phys. Rev. D 5 Ž wx 4 O. Costa de Beauregard, C. Imbert, Y. Levy, Phys. Rev. D 15 Ž wx 5 J.P. Hugonn, R. Pett, J. Optcs 8 Ž Ž n French.. wx 6 W. Nasalsk, J. Opt. Soc. Am. A 13 Ž wx 7 E. Hecht. Optcs, Addson Wesley, wx 8 W.T. Welford, Aberratons of optcal systems, Adam Hlger, wx 9 M.A. Porras, J. Alda, E. Bernabeu, Appl. Opt. 31 Ž w10x J. Alda, J. Alonso, E. Bernabeu, J. Opt. Soc. Am. A 14 Ž w11x F. Falco, T. Tamr, J. Opt. Soc. Am. A 7 Ž w1x J-J Greffet, C. Baylard, Opt. Commun. 93 Ž w13x B.R. Horowtz, T. Tamr, J. Opt. Soc. Am. 61 Ž

( ) + + REFLECTION FROM A METALLIC SURFACE

( ) + + REFLECTION FROM A METALLIC SURFACE REFLECTION FROM A METALLIC SURFACE For a metallc medum the delectrc functon and the ndex of refracton are complex valued functons. Ths s also the case for semconductors and nsulators n certan frequency

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Boundaries, Near-field Optics

Boundaries, Near-field Optics Boundares, Near-feld Optcs Fve boundary condtons at an nterface Fresnel Equatons : Transmsson and Reflecton Coeffcents Transmttance and Reflectance Brewster s condton a consequence of Impedance matchng

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Problem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that

Problem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that PHY 387 K. Solutons for problem set #7. Problem 1: To prove that under the assumptons at hand, the group velocty of an EM wave s less than c, I am gong to show that (a) v group < v phase, and (b) v group

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations

MAGNUM - A Fortran Library for the Calculation of Magnetic Configurations CRYO/6/34 September, 3, 6 MAGNUM - A Fortran Lbrary for the Calculaton of Magnetc Confguratons L. Bottura Dstrbuton: Keywords: P. Bruzzone, M. Calv, J. Lster, C. Marnucc (EPFL/CRPP), A. Portone (EFDA-

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of

More information

16 Reflection and transmission, TE mode

16 Reflection and transmission, TE mode 16 Reflecton transmsson TE mode Last lecture we learned how to represent plane-tem waves propagatng n a drecton ˆ n terms of feld phasors such that η = Ẽ = E o e j r H = ˆ Ẽ η µ ɛ = ˆ = ω µɛ E o =0. Such

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

Introductory Optomechanical Engineering. 2) First order optics

Introductory Optomechanical Engineering. 2) First order optics Introductory Optomechancal Engneerng 2) Frst order optcs Moton of optcal elements affects the optcal performance? 1. by movng the mage 2. hgher order thngs (aberratons) The frst order effects are most

More information

Lecture 3. Interaction of radiation with surfaces. Upcoming classes

Lecture 3. Interaction of radiation with surfaces. Upcoming classes Radaton transfer n envronmental scences Lecture 3. Interacton of radaton wth surfaces Upcomng classes When a ray of lght nteracts wth a surface several nteractons are possble: 1. It s absorbed. 2. It s

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Differentiating Gaussian Processes

Differentiating Gaussian Processes Dfferentatng Gaussan Processes Andrew McHutchon Aprl 17, 013 1 Frst Order Dervatve of the Posteror Mean The posteror mean of a GP s gven by, f = x, X KX, X 1 y x, X α 1 Only the x, X term depends on the

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

Electron-Impact Double Ionization of the H 2

Electron-Impact Double Ionization of the H 2 I R A P 6(), Dec. 5, pp. 9- Electron-Impact Double Ionzaton of the H olecule Internatonal Scence Press ISSN: 9-59 Electron-Impact Double Ionzaton of the H olecule. S. PINDZOLA AND J. COLGAN Department

More information

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR 5.0, Prncples of Inorganc Chemstry II MIT Department of Chemstry Lecture 3: Vbratonal Spectroscopy and the IR Vbratonal spectroscopy s confned to the 00-5000 cm - spectral regon. The absorpton of a photon

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

CHAPTER II THEORETICAL BACKGROUND

CHAPTER II THEORETICAL BACKGROUND 3 CHAPTER II THEORETICAL BACKGROUND.1. Lght Propagaton nsde the Photonc Crystal The frst person that studes the one dmenson photonc crystal s Lord Raylegh n 1887. He showed that the lght propagaton depend

More information

ECE 107: Electromagnetism

ECE 107: Electromagnetism ECE 107: Electromagnetsm Set 8: Plane waves Instructor: Prof. Vtaly Lomakn Department of Electrcal and Computer Engneerng Unversty of Calforna, San Dego, CA 92093 1 Wave equaton Source-free lossless Maxwell

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula

More information

Supporting Information

Supporting Information Supportng Informaton Water structure at the ar-aqueous nterface of dvalent caton and ntrate solutons Man Xu, Rck Spnney, Heather C. Allen* allen@chemstry.oho-state.edu Fresnel factors and spectra normalzaton

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Pulse Coded Modulation

Pulse Coded Modulation Pulse Coded Modulaton PCM (Pulse Coded Modulaton) s a voce codng technque defned by the ITU-T G.711 standard and t s used n dgtal telephony to encode the voce sgnal. The frst step n the analog to dgtal

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Modelli Clamfim Equazione del Calore Lezione ottobre 2014 CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g

More information

8.592J: Solutions for Assignment 7 Spring 2005

8.592J: Solutions for Assignment 7 Spring 2005 8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

New Method for Solving Poisson Equation. on Irregular Domains

New Method for Solving Poisson Equation. on Irregular Domains Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

Programming Project 1: Molecular Geometry and Rotational Constants

Programming Project 1: Molecular Geometry and Rotational Constants Programmng Project 1: Molecular Geometry and Rotatonal Constants Center for Computatonal Chemstry Unversty of Georga Athens, Georga 30602 Summer 2012 1 Introducton Ths programmng project s desgned to provde

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton

More information

Homework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich

Homework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich Homework 4 Contact: frmmerm@ethz.ch Due date: December 04, 015 Nano Optcs, Fall Semester 015 Photoncs Laboratory, ETH Zürch www.photoncs.ethz.ch The goal of ths problem set s to understand how surface

More information

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION

CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR OF SEPARATION PRODUCTS OUTPUT DETERMINATION Górnctwo Geonżynera Rok 0 Zeszyt / 006 Igor Konstantnovch Mladetskj * Petr Ivanovch Plov * Ekaterna Nkolaevna Kobets * Tasya Igorevna Markova * CHARACTERISTICS OF COMPLEX SEPARATION SCHEMES AND AN ERROR

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information

Color Rendering Uncertainty

Color Rendering Uncertainty Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:

More information

Fresnel's Equations for Reflection and Refraction

Fresnel's Equations for Reflection and Refraction Fresnel's Equatons for Reflecton and Refracton Incdent, transmtted, and reflected beams at nterfaces Reflecton and transmsson coeffcents The Fresnel Equatons Brewster's Angle Total nternal reflecton Power

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Note: Please use the actual date you accessed this material in your citation.

Note: Please use the actual date you accessed this material in your citation. MIT OpenCourseWare http://ocw.mt.edu 6.13/ESD.13J Electromagnetcs and Applcatons, Fall 5 Please use the followng ctaton format: Markus Zahn, Erch Ippen, and Davd Staeln, 6.13/ESD.13J Electromagnetcs and

More information

Effect of Losses in a Layered Structure Containing DPS and DNG Media

Effect of Losses in a Layered Structure Containing DPS and DNG Media PIERS ONLINE, VOL. 4, NO. 5, 8 546 Effect of Losses n a Layered Structure Contanng DPS and DNG Meda J. R. Canto, S. A. Matos, C. R. Pava, and A. M. Barbosa Insttuto de Telecomuncações and Department of

More information

Conservation of Angular Momentum = "Spin"

Conservation of Angular Momentum = Spin Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Summer Term 01, Abbe School of Photoncs, FSU Jena, Prof. Thomas Pertsch Computatonal Photoncs Semnar 03, 7 May 01 Implementaton of the Matr Method calculaton of the transfer matr

More information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

Solutions to Problem Set 6

Solutions to Problem Set 6 Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem 4.7.7 Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral:

More information

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

GEO-SLOPE International Ltd, Calgary, Alberta, Canada   Vibrating Beam GEO-SLOPE Internatonal Ltd, Calgary, Alberta, Canada www.geo-slope.com Introducton Vbratng Beam Ths example looks at the dynamc response of a cantlever beam n response to a cyclc force at the free end.

More information

Microscopy with self-reconstructing beams

Microscopy with self-reconstructing beams Mcroscopy wth self-reconstructng beams Floran O. Fahrbach 1,, Phlpp Smon und Alexander Rohrbach 1, 1 Centre for Bologcal Sgnallng Studes (boss), Unversty of Freburg, Germany Laboratory for Bo- and ano-photoncs,

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

On the correction of the h-index for career length

On the correction of the h-index for career length 1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat

More information

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Important Instructions to the Examiners:

Important Instructions to the Examiners: Summer 0 Examnaton Subject & Code: asc Maths (70) Model Answer Page No: / Important Instructons to the Examners: ) The Answers should be examned by key words and not as word-to-word as gven n the model

More information

One Dimensional Axial Deformations

One Dimensional Axial Deformations One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

Microlens quality assessment using the Extended Nijboer-Zernike diffraction theory

Microlens quality assessment using the Extended Nijboer-Zernike diffraction theory Mcrolens qualty assessment usng the Extended Njboer-Zernke dracton theory J.J.M. Braat, S. van Haver, S.F. Perera Delt Unversty o Technology Department o Imagng Scence and Technology Optcs Research Group

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Frequency dependence of the permittivity

Frequency dependence of the permittivity Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information