Paraxial Beam, Gaussian Basics

Size: px
Start display at page:

Download "Paraxial Beam, Gaussian Basics"

Transcription

1 Paraxia Beam, Gaussian Basics ECE 5368/6358 han q e - copyrighted Use soey for students registered for UH ECE 6358/5368 during courses - DO NOT distributed (copyrighted materias). Introduction. Paraxia approximation for Gaussian beam (can skip). Important properties of Gaussian beam. Beam description From the above (we drop the subscript from k for simpicity):

2 LectSet - Gaussian beam basic 5358_p.nb E z, r = p wk  Q z  k z+ r Q z (..a) where: Q z = z - ; = kw = p w (..b) Note: we have E[z,r] in the unit of so that the normaized profie: ength p E z, r r r q= is unitess. Rea E has rea physica unit, of course. The formua is ony for mathematica convenience. We can aso drop the normaization factor for convenience. It doesn t have any specia meaning. We can aso write: where: p E z, r = wk Abs Q z  -k r Im Q z  k z+ r Re - Arg Q z Q z (..a) Q z = z - ; = kw Re = z-â = p w z ; Im = z + z-â and: Arg Q z =-ArcTan z ; z + What are the meaning of a these formuas? What can we infer?. Beam intensity.. Transverse profie For beam intensity: I E = Abs I w k Abs Q z It is a Gaussian profie: with a beam intensity radius: Substitute: And with We obtain: W z = wk Abs Q z  -k r Im Q z  k z+ r Re k Im - k r Im Q z = w k Abs Q z - ==> Im Q z W z = z + = k k = kw W z = w + z w k Abs Q z - r W z Q z = kw z = (..b) (..c) Q z - Arg Q z (..) r W z (..) p W z (..3) + z (..4) = W + z (..5) The transverse profie of a Gaussian beam is a Gaussian profie everywhere, with the radius changing as W diameter). + z. We aso refer to beam spot size as times the beam radius (beam Gaussian beam does not originate from a point: there is no such thing as a point source of ight, but there is a pace where the Gaussian beam has the smaest profie. The smaest beam radius

3 LectSet - Gaussian beam basic 5358_p.nb 3 is caed beam waist: W[]. (see note above aso that shows that the point source of a Gaussian beam is a compex number) HW: pot beam intensity Manipuate k = p ; =.5 k w ;W= w + z ; Show Pot -w + v, w + v, v, -5, 5, AspectRatio Æ.5, PotRange Æ -5, 5, Fiing Æ Æ, Hue.4,,,.7, Graphics Red, Thickness.5, Line z, -5, z, 5, ImageSize Æ 5, 75, PotRange Æ -5, 5, -5, 5, DensityPot w k z + - x +y W, x, -5, 5, y, -5, 5, Mesh Æ Fase, w k CoorFunction Æ GrayLeve, PotRange Æ,, PotPoints Æ 5, 5 z +, z, -5, 5, w,,.5, 4 z w,.. Beam divergence (focusing) Obviousy, at arge z, the beam radius diverges as: W + z Thus, one can define the divergence ange q with: Tan q ~ Sin q = w z z ~ wz (..6) = w (..7a) When we dea with sma ange: Tan q ~ Sin q ~q Aso, with -> p w Tan q ~ Sin q = = = (..7b). p w p w p D where D is the diameter of the beam. As it turns out, this formua with a sight modification is amost generic

4 4 LectSet - Gaussian beam basic 5358_p.nb for any beam: Tan[Divergence ange] = Tan q = m (..7c) D where m is some number that is a characteristic of a specific beam profie, and which is for the case of p Gaussian. It is. for exampe for a top hat beam. (there is no hard definition for the divergence ange either, so the number m is not something fundamenta). Manipuate z w ; Pot w z z,w z, z,,, AspectRatio.5, z PotRange 8, 8, Fiing, Hue.4,,,.7, w,.5, 3 For sma ange: Often, we tak about beam fu width: qº p D D qº 4 p (..8) An exampe for non-gaussian beam: for top hat beam, instead of 4 ~.7, it is.44. p Let s think: sma waist -> arge divergence arge waist -> sma divergence As seen above: waist x divergence = constant. This gives rise to the concept of etendue and reated to the definition of brightness. A simiar formua can be used for focusing (see more in section 4)...3 Longitudina profie Again, for the beam intensity above: I E = Abs I w k Abs Q z wk Abs Q z  -k r Im Q z  k z+ r Re - k r Im Q z = But now, we are interested in variation aong z. At r = the center of the beam w k = w k = w k Abs Q z z + w k Abs Q z - Q z - Arg Q z (..) r W z (..) = I + z z As z->, the intensity drops as as expected for any spherica wave z + z (..9) As the transverse profie of a Gaussian beam gets arger (or smaer), the intensity drops (or increases) as (think Lorentzian) by virtue of power conservation. + z

5 LectSet - Gaussian beam basic 5358_p.nb 5 Pot, z, -,, PotRange Æ A + z Power conservation We now can express the beam intensity I - + z z r W z (..) where I[] is the center of the beam intensity at its waist. (the highest intensity or peak intensity of the beam), Power is obtained by integrating over a area: Int + z - r W z pr r ConditionaExpression p Int W z,re z + W z > FuSimpify p Int W z p Int w z +. W z -> w + z Power is: P = I p w = I A (..a) where: A =pw (..b) is the beam "area". We may ask why the factor? It is because I[] is the peak intensity at the center of the beam waist: I = P p w (..) The power is an integration over a intensity of the whoe beam, not just the peak. Hence, we can think of (..a) as: P = I ave A (..3) with average intensity I ave = I. Think ike the area of a triange. The area is / of the height times the base. Sometime we write the beam intensity as:

6 6 LectSet - Gaussian beam basic 5358_p.nb I r, z = P p w - + z r W z or P p W z - r W z (..4).3 Beam phase.3. Phase front Let's ook at the beam again: E z, r = wk Abs Q z  -k r Im Q z  k z+ r Re - Arg Q z Q z (..a) = p w where: Q z = z - ; = kw z Re = ; Im = z-â z + z-â and: Arg Q z =-ArcTan ; z The phase term is: k z +r Re The phase front is defined as: k z +r Re Q z Q z z + - Arg Q z (.3.) - Arg Q z =constant. (.3.) This surface is a function of z and r. What are the meanings of each term? kz is the pane wave term. k r Re (..b) (..c) Q z the wavefront. Arg Q z =-ArcTan is indepent of r, is just a shift that is dependent on z. z ceary describes the curvature of If we have a pane wave traveing aong a Gaussian wave of the same k, do they traveing in step on the axis?.3. Guoy effect Not quite. Because for pane wave, it is k z ; for Gaussian, it is k z + ArcTan. So as z goes from - to, z The Gaussian is retarded by a p phase shift. This is caed Guoy's effect..3.3 Phase front curvature Phase front term: k z +r Re Q z - Arg Q z (.3.) At arge z, where z>>r, the variation of the term is ~ z +r Re Q z But this is the equation for a parabooida surface: r = constant (.3.3) z + = constant (.3.4) R where R is the radius of curvature. Therefore, we can define: Re = (.3.5a) Q z R z where R z is the radius of curvature of the phase front.

7 LectSet - Gaussian beam basic 5358_p.nb 7 From (..b) Re Q z = Re z-â = R z = z + Now we find out what the meaning of Re Q z = R z ; z : Q z and Im Q z = kw z = Thus, we can write: Q z z z + (.3.5b) (.3.5a) p W z (.3.6) = +Â R z p W z (.3.7) (..b) Remember the comparizon of paraxia spherica wave vs Gaussian: in spherica becomes in Gaussian. Since z Q z in spherica is the radius of curvature, it makes sense that the Re part of is aso the radius of curvature. z Q z Remember: is NOT the Rayeigh range of the beam. We reserve that term ony for, where W[] is p W z p W the beam waist. z =.; k = Pi; w = ; z = kw ; x z phasefrontpot = ContourPot Mod k z + z + z - Arg z - z,*pi == Pi, z,,, x, -5, 5, ImageSize ->,, PotPoints Æ 5, ContourStye Æ Red ; profiepot = Pot -w + z z, w + z, z,,, ImageSize ->,, z PotRange -> -5, 5, PotStye Æ RGBCoor.,., ; Show phasefrontpot, profiepot Ü Graphics Ü The term, where Q z = z -Â z Q z, tes us the beam radius of curvature in the rea part, and the beam transverse profie radius in the imaginary part.

8 8 LectSet - Gaussian beam basic 5358_p.nb 3. Practica techniques to measure Gaussian (or non- Gaussian) beam Discuss in cass ony, and for those who reay dea with aser beam in the ab. 3. Intensity profiing: beam scan method How do we measure the profie of a beam? Why don't we use things ike scanning sit? Scanning of a Gaussian beam W u x y W x - y W y - p W p W p W - erf u W Pot - erf u W. W ->, u, -.5,.5, GridLines -> Automatic, Frame -> True, AxesOrigin ->,, PotRange -> -.5,.5,,

9 LectSet - Gaussian beam basic 5358_p.nb 9 N - erf u W. u -> W N - erf u W. u -> -W N - erf u W. u -> -W N - erf u W. u -> W FindRoot - erf x ==.5, x,.5 x Ø.843 x. x Æ ` 3. Beam divergence (not at thin ens foca pane) How do we determine the divergence of a beam? Not a aser beams are Gaussian (remember that). But first, what is divergence? What is far fied. The distribution of a beam intensity as a function of ange at infinity is defined as FF. We define a function: f[q,f]. Sometime, it is written: f[w], where W={q,f}. If we ook at a Gaussian beam, the resuts we saw above is that the beam profie is - r W + z - r W + z We define. We need to convert this into ange {q,f} as zø. Ø - r W z = - W = Sin q r W z = - z Remember that Sin q = W = W p W Sin q W Sin q (aso use Tan q ). Thus: - Sin q : this is the far fied profie function = p W. So now we have a practica question: how do we determine a beam far-fied? Go out very far, measure the spatia profie, and divide by distance z? how far? cm? m? km? The key condition is z >> : Rayeigh range. Rayeigh range is the range where the beam diverges very itte and the wavefront is neary panar: it is where the pane wave approximation is good. Suppose we have a Gaussian beam with -cm diameter (that is, beam waist=.5 cm), what is the beam divergence? What is its Rayeigh range?

10 LectSet - Gaussian beam basic 5358_p.nb We need to know the waveength aso. Let's =.638 mm; the beam Rayeigh range is:.5 cm ^.638 m. m ^ 4 cm 4. cm Suppose we want to measure at a ocation ~ times the Rayeigh range, it' be: 4 m or. km! How big is the beam there? w + z.5 cm. w ->.5 cm. z -> To measure the far fied of any beam, we can just focus the beam with a ens, then measure the profie at a reasonabe distance (do not confuse this with measuring the foca pane of a thin ens that we wi earn ater on). But then, the divergence we get is just the divergence of the focused beam, not the origina beam! How do we convert? We can convert for a Gaussian beam! But first, et's ook at a concept: etendue. For a beam or an optica system, the product of accepting area and accepting soid ange is caed beam etendue. We can have -D etendue or -D etendue. Earier, we ask a question about Gaussian. Let's see: Beam divergence: Tan q = or q p W º ; p W Beam waist: W ; Thus, -D etendue = W = 4 : constant! p W p (What is its -D etendue?) So now, if a ens transforms a Gaussian beam into another Gaussian beam (just of different beam waist), then if we measure the etendue of one, we knows the etendue of the other! Appy this to the aboratory exercise. 3.3 Laboratory exercise: M and number of times of diffraction-imited Ever wonder when you focus a beam ike aser beam, what is at the focus spot? NOT a beams wi give the same spot size even with the same focusing convergence ange! Perform a Gaussian beam measurement in the ab and report your measurement: ocate and measure beam waist, measure divergence ange. Focus the beam with a ens and measure. 4. Gaussian beam transformation through enses (and parabooida/spherica mirrors) In practice, rarey any beam is used without any transformation through optics.

11 LectSet - Gaussian beam basic 5358_p.nb Except for integrated optics, which use waveguide to conduct a beam, and the beam profie is determined by the waveguide modes, many systems sti use free-space optics. Even at a very macroscopic eve... The beam profie properties are essentia. 4. Review of thin ens and spherica or parabooida mirror

12 LectSet - Gaussian beam basic 5358_p.nb Thin ens Lighty curved mirror Thin ens and ow-curvature spherica or parabooida mirror has an interesting property: they convert input r - k beam into an output beam with an approximatey parabooida wavefront: f for thin ens where f is the ens foca point and  k r R M where R M is the mirror radius of curvature. (derivation comes from paraxia approximation). Convention: f> for a positive ens, and R M < for concave, R M > for convex. Derivation for this comes from ray approximation in which the phase at ocation r is proportiona to r. (Linked to Simpe Optica Eements). Thus, if a ens is ocated at ocation z, and the input wave is U x, y, z, then the output wave at the same r - k ocation is: U x, y, z f. Whatever happens to the wave afterward? at arbitrary z? this is precisey the interesting property of these optics. In particuar, we wi see that Gaussian beam behaves quite nicey with this type of optics and that is the reason we study it here. 4. Gaussian beam after a thin ens. (derivation) A Gaussian beam remains a Gaussian through a thin ens (or equivaent paraboic mirror) but with a different beam waist ocation, vaue, or direction of propagation (the case of mirror). Thin ens What becomes of the beam here? z Let's set z= at the waist of the input beam. Then, et s write the input Gaussian profie at z just before the ens. From: to: E z, r = wk Abs Q z  -k r Im Q z  k z+ r Re - Arg Q z Q z (..a) E z, r -k r Im Q z  k z +r Re - Arg Q z Q z (4..)

13 LectSet - Gaussian beam basic 5358_p.nb 3 w k (We drop the term which is just an overa ampitude factor without any effect on profie and phase Abs Q z  for simpicity). Remember that: Im = = (4..a) Q z kw z p W z and: Re = (4..b) Q z R z Then, by substitute (4..), (4..) becomes: - r W z The output beam just after the ens is: - r W z  k z + r R - Arg Q z  k z + r - Arg Q z R - k r f (4..3) (4..4) = - r W z  kz  k r R - f - Arg Q z If we define: = -, (4..5) R' R f then we can write: - r W z  kz  k r R' - Arg Q z (4..6) Let's ook at the expression. What does it te us? Obviousy, it has a Gaussian ampitude profie: - r W z ; and a parabooida phase term:  kz  k r R'. This is a typica Gaussian profie! (note: The phase term - Arg Q z wk, ike the ampitude term that we drop, is Abs Q z  ONLY a constant phase term without any r-dependency hence, just a factor). The ony difference between this output Gaussian and the input Gaussian is that the phase front has a radius of curvature R' instead of R. They both have the same beam spot radius at z : W z. Thus, this is aso a soution of the Hemhotz equation. The output of a Gaussian beam through a an on-axis thin ens is a Gaussian beam. BUT this Gaussian beam is NOT the same as the input (inspite of the fact that BOTH have the same beam spot radius at z : they have difference phase front at z ). So, to find out the property of this output Gaussian from the ens, we must find a Gaussian beam expression that gives us exacty the beam spot radius: W z and radius of curvature R' at ocation z. Let's write the expression again: - r W z  kz  k r R' - Arg Q z =  k r + R' kw z  kz - Arg Q z Now we reca that in genera, + = R' kw z P z (4..7) where P z is the Q-function (virtua point) for some Gaussian beam. How do we do that? Basicay we must find the beam waist ocation z and beam waist W '. We write the genera expression for a Gaussian beam: - r r  k z+ W z R' z -Arg z-â' (4..8) This is expression for a beam with waist at z=. We shift the coordinate to z : - r W' z-z  k z-z + r R' z-z -Arg z-z - ' (4..9)

14 4 LectSet - Gaussian beam basic 5358_p.nb At z = z : r W' z -z  k z -z + r R' z -z -Arg z -z - ' (4..) Now, this beam must match with the output Gaussian beam from the ens at z, right? So, we compare terms by terms: - Ampitude profie: - r W' z -z = - r W z (4..) which yieds:w ' z - z = W z ; or: W ' + z -z = W ' + z W ' + z -z p W = W ' + z 4 W ' + z -z = W p W ' + z p W p W 4 p W ' + z -z = p W p W ' + z p W ' + z -z ' - Phase profie: r = + z (4..)  k z -z + -Arg z R' z -z -z - '  kz  k r - R z f - Arg z - (4..3) Again, separate the r-dependent term from the others:  k r R' z -z =  k r - R z f f = - (4..4) R' z -z R z f which is precisey the wavefront reation we mention earier. The other term is NOT important here since it is just a constant phase shift (no r or z-dependence). Recaing that = z, we have the equation: R z z + z -z = z - (4..5) z -z + ' z + f Soving the equations Let's rewrite the Eqs. again: ' + z -z =z ' + z (4..) z -z z z -z + ' z + f (4..5) Changing z - z -> -d ' + d =z ' + z (4..6a) d z - d + ' f z + (4..6b) So we have eqs. with unknowns.

15 LectSet - Gaussian beam basic 5358_p.nb 5 FuSimpify Sove zp + d zp == z + z z d, == d + zp f - z, d, zp z + z Genera::spe : Possibe speing error: new symbo name "z" is simiar to existing symbo "zp". d Ø z - f f + f, zp Ø f - z f +z + z f z f - z f +z + z 4.3 Resut discussion 4.3. Reations between the beams The first resut: ' = f z - f + or (4.3.) W ' = W f z - f + = W M (4.3.) where we define: M = f z - f + = f z - f + z - f (4.3.3) Next: d Ø z - f f + f f d - f = z f f - z f +z - f +z z - f +z d = f + z - f M (4.3.4) The quantity M is interesting: it is the ratio between the beam waists. Thus we define M to be the "magnification" factor (as if we use the ens to magnify or focus the beam). We see these reations are somewhat simiar to ray-optics of ens that we are famiiar. In fact we have other reations invoving M: f Rayeigh range: ' =z = M z z - f (4.3.5) + Beam divergence: q' = = = q. (4.3.6) p W ' M p W M Write a code cacuate output beam waist and ocation as a function of input beam parameters 4.3. Meaning of the "magnification" (or reduction) term

16 6 LectSet - Gaussian beam basic 5358_p.nb Let's ook at the term: f M = (4.3.3) z - f + If the beam is point-ike: it has a sma beam waist and short Rayeigh range, then: M = f z - f + = f z - f º + z - f f z - f (4.3.7) z If z - f >> so that º. z - f But this is the typica ray-optics ens reation: z d z-f f Image magnification: W ' = W M In fact, we ook at other reation: d - f = z - f M or: d - f = z - f d - f = f z - f f z - f d - f z - f = f dz - f d + z + f = f d+z dz = f f z + d = f (4.3.8) This is the ens imaging reation in geometrica (ray) optics. Depth of focus: ' = M (4.3.9) A of these reations are compatibe with ray optics. When is this concept not reevant? obviousy, when the beam waist is arge: the beam is more "pane-wave" ike rather than spherica wave emanating from a spot: now we know when we can use simpe geometrica optics and when we shoud not when deaing with Gaussian beam: the Rayeigh ranges << f or the waist is very sma Gaussian beam focusing or mode matching Now we have the other extreme: our beam is penci ike, which means it has a arge Rayeigh range: >> f. Aso, the beam waist shoud aso be cose to the ens, in other words, z is NOT >> f. What wi happen? Let's ook at beam waist: W ' f = W = W z - f + f z - f + = f p W z - f +

17 LectSet - Gaussian beam basic 5358_p.nb 7 Since we assume that >> z - f, we can write: W ' º p W f or the famous Gaussian beam transform formua: p W W ' =f (4.3.) Notice how the beam waists are inverse of each other: input a arge beam, we have a very sma focused spot. Input a sma spot, we have a arge focus spot Foca-foca conjugate Now suppose we want the "symmetric case". Can we make a beam such that the input is simiar to the output except for a dispacement of the waist? What if we just put a Gaussian waist at the foca pane of a ens? That is: z = f or z - f =? From the above: d = f + z - f M = f This means that the output beam waist is at the other foca pane. Sure, this is simpe enough to understand: if we put a beam waist at one foca pane, the ens just bend the phase such that the output is at the other foca pane. But what happens to beam waist size? For that, we need: M = f z - f + = f If we want the two beams aso to be neary equa in waist, we need f ~, in other words, the Rayeigh range of the beam has to be neary equa to the ens f. So we can't arbitrary make beams ike that for a given ens. Either we choose the ens f or we must have beam with proper waist!

18 8 LectSet - Gaussian beam basic 5358_p.nb =. ; f = 5. ; ens = Graphics Hue.45,,,, Disk,, 5, 5, Hue,,,, Thickness.5, Line -,,,, Purpe, PointSize., Point - f,, Point f, ; Manipuate z = p w ; M = f -z - f +z ; d = f - z + f M ;w= M w ; z = p w ; profie = Pot -w z - z +,w + z z - z z, z, -,, PotRange Æ A, PotStye Æ RGBCoor.,.,, Fiing Æ Æ, Hue,.7,,. ; waist = Graphics Red, Thickness.5, Line z, -w, z, w ; profie = If d., Pot -w z - d +,w + z - d z z, z,,, PotRange Æ A, PotStye Æ RGBCoor.,.,, Fiing Æ Æ, Hue.6,,,. z - d, Pot -w +,w + z z - d z, z, d,, PotRange Æ A, PotStye Æ RGBCoor.,.,, Fiing Æ Æ, Hue.6,,,.3 ; waist = Graphics Bue, Thickness.5, Line d, -w, d, w ; Show ens, profie, waist, profie, waist, PotRange Æ -,, -5, 5, z, -,, w, 3,.5, 5 Consider focusing a beam: From the formua: p W W ' =f (4.3.) W ' = f p W If we want a very sma focused spot, we want W as arge as possibe (and f as sma as possibe). But the input beam can t be arger that the ens! otherwise, the power wi not be coected. Hence, we can say that W is imited by the radius of the ens r. Hence: W ' ~ f p r

19 LectSet - Gaussian beam basic 5358_p.nb 9 But the ratio f ~, which is aso caed nunerica aperture (NA) of the ens. We see that it doesn t matter arge r Sin q ens or sma ens, ong f or short f, the ony thing that matter is the NA=Sin[q] The above enses give exacty the same spot size when focusing at optima beam width. W ' ~ (4.3.) p NA That s why NA is the ony rating that matter in enses when it comes to resoution: how fine the detais we can see or how sma a spot we can focus Off-axis beam (advanced ony) Linear decomposition approach (advanced ony) BeamOut beamin_list, ensparm_list, D_, x_, y_, z_ := Modue, W, z, Wz, Rz, R, f,z,k, P, Ta, f,z = ensparm;, W = beamin ; k = *Pi ;Ta=D f ; z = k W ; Wz = W + z z z Rz = ; z + z R = Rz - f ; P = R + ; k Wz ; p P z + P kz kx k y-d-p Ta z+p z+p - k P Ta

20 LectSet - Gaussian beam basic 5358_p.nb BeamIntOut beamin_list, ensparm_list, D_, x_, y_, z_ := Modue, W, z, Wz, Rz, R, f,z,k, P, Ta, f,z = ensparm;, W = beamin ; k = *Pi ;Ta=D f ; z = k W ; Wz = W + z z z Rz = ; z + z R = Rz - f ; R - P = R + k Wz ; 4 k Wz 4 ; R - z + R + k Wz = 4 k Wz 4 z R R - k Wz 4 k Wz R + 4 k Wz 4 Abs p P z + P kx k y-p Ta z+p z+p ens =, *^4 ; * unit micron * beamin =.5,.3*^4 * unit micron * ; Pot3D Re BeamOut beamin, ens,,, y, z, z, 998,, y, -,, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü

21 LectSet - Gaussian beam basic 5358_p.nb ens =, *^4 ; * unit micron * beamin =.5,.3*^4 * unit micron * ; Pot3D Re BeamOut beamin, ens, 3,, y, z, z, 998,, y, -,, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü ens =, *^4 ; * unit micron * beamin =.5,.3*^4 * unit micron * ; Pot3D Re BeamOut beamin, ens,,, y, z, z, 998,, y, -,, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü ens =, *^4 ; * unit micron * beamin =.5,.5 * unit micron * ; Pot3D Re BeamOut beamin, ens,,, y, z, z, 998,, y, -,, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü

22 LectSet - Gaussian beam basic 5358_p.nb ens =, *^4 ; * unit micron * beamin =.5,.5 * unit micron * ; Pot3D Re BeamOut beamin, ens,,, y, z, z, 9 98,, y, -,, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü ens =, *^4 ; * unit micron * beamin =.5,.5 * unit micron * ; Pot3D Re BeamOut beamin, ens, 3,, y, z, z, 9 98,, y, -3, -99, PotPoints ->,, BoxRatios ->,,.5, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü ens =.5, 4.9 *^4 ; * unit micron * beamin = 4.8, * unit micron * ; Pot3D Abs BeamOut beamin, ens,,, y, z ^, z, 5, 5 5, y, -,, PotPoints ->, 5, BoxRatios -> 5, 4, 4, Mesh -> Fase, PotRange -> A, ViewPoint ->,, Ü SurfaceGraphics Ü

23 LectSet - Gaussian beam basic 5358_p.nb 3 When you go into the ab and try to focus a beam to as a sma spot as possibe, what do you do? Project: use whatever enses you can find in the ab, see who can obtain the smaest focused spot. Dr. Bujin Guo is using a pinhoe and try to focus a beam through. Obtain from him necessary parameters to cacuate how much the throughput is (% of power that gets through the pin hoe). Look at the diode beam output. Why do you think the beam diverges in one dimension so much more than the other? Write a theory for x-y asymmetric beam. Diode beam focus: how do you coupe a beam from a diode into a rectanguar waveguide? (which orientation to use?) Where is the beam focused? d - f = z - f d = f + z - f f z very near the foca point. f = z z - f +z - f f º f z z - f z + º z - f f If we want to coimate and expand the beam to shoot far away, where do we put the beam? z f use the reation: d - f = z - f. If the input beam waist is right at the foca point: z z - f +z = f, then: d - f = : the output beam waist is at the other foca point Numerica aperture of a ens If we want a ens of focus or coimate a beam, how do we order it? Let's ook at the externa cavity: Grating: N = g/mm, 5 5 6mm 3 /N [(sin + sin )] = = 75 º, sin =.966, = 5.um sin =.6, = 3.º Lens:.7mm ZnSe f =.7 mm, NA= º 3. º 3 º.5 º w x w y L L ink

24 4 LectSet - Gaussian beam basic 5358_p.nb Does the ens diameter matter? Why do we want "aspheric" ens? what is the reation between distance L and L and foca ength f? Where shoud the ocation of the ens shoud be if we are given a fixed distance L +L. What can the atera sit do for us? Summary Now we know why Gaussian beam is interesting and usefu: it remains Gaussian (athough transformed) under enses or spherica/parabooida mirrors. (A the above resuts for ens are aso appicabe with mirrors). With appropriate condition, we can find a set of Gaussian soutions that remain invariant in a system of optics: aser cavity. We' see that Gaussian turns out to be just one soution among many cavity modes; not surprisingy, it is the simpest soution, the owest order soution of the Hermite-Gaussian function famiy. More importanty: many concepts above are NOT excusive to the Gaussian beam. Even non- Gaussian beams have simiar behaviors: this is because they a are governed by diffraction behavior. What we earn is the diffraction of paraxia ight beam. Gaussian is just a specia case with neat mathematica expression. In the vertica dimension, a diode beam coming out of its facet has a fat wavefront and an ampitude profie ike this: a -a x cos[ k a] e z cos[ k a] e x a cos[ k x] x a Pot the wave profie of the beam as it eaves the diode facet aong the x and z dimension. Use resonabe assumption of everything. (you must figure it out yoursef). Assignment: do a the excersizes in section 4.3 above and this: John measured a beam with a Gaussian intensity profie - r w and obtained w=.5 mm. The waveength was = mm. He then measured the beam divergence which has a fu width haf max (FWHM) of 35 degree. Is this posibe? TanT = p w. w ->.5. ->.73

25 LectSet - Gaussian beam basic 5358_p.nb 5 Pot - Tan q*p 8.7 ^, q, -4, 4, PotRange -> A Advanced demonstration: Large NA beam More accurate cacuation needed for very arge NA beam. Poarization needs be incuded. Iustration:

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Multiple Beam Interference

Multiple Beam Interference MutipeBeamInterference.nb James C. Wyant 1 Mutipe Beam Interference 1. Airy's Formua We wi first derive Airy's formua for the case of no absorption. ü 1.1 Basic refectance and transmittance Refected ight

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

(Refer Slide Time: 2:34) L C V

(Refer Slide Time: 2:34) L C V Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Interference and Diffraction

Interference and Diffraction 0_Interference.nb Chapter 0 Interference and Diffraction 0. Interference 0.. Two Sit Interference Pattern We have aready seen exampes of the interference between two waves. Given two counter propagating

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions Physics 116C Hemhotz s an Lapace s Equations in Spherica Poar Coorinates: Spherica Harmonics an Spherica Besse Functions Peter Young Date: October 28, 2013) I. HELMHOLTZ S EQUATION As iscusse in cass,

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

Keywords: Rayleigh scattering, Mie scattering, Aerosols, Lidar, Lidar equation

Keywords: Rayleigh scattering, Mie scattering, Aerosols, Lidar, Lidar equation CEReS Atmospheric Report, Vo., pp.9- (007 Moecuar and aeroso scattering process in reation to idar observations Hiroaki Kue Center for Environmenta Remote Sensing Chiba University -33 Yayoi-cho, Inage-ku,

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE Progress In Eectromagnetics Research, PIER 30, 73 84, 001 CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE W. Lin and Z. Yu University of

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

LECTURE 10. The world of pendula

LECTURE 10. The world of pendula LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive

More information

Lecture 8 February 18, 2010

Lecture 8 February 18, 2010 Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some

More information

Problem Set 6: Solutions

Problem Set 6: Solutions University of Aabama Department of Physics and Astronomy PH 102 / LeCair Summer II 2010 Probem Set 6: Soutions 1. A conducting rectanguar oop of mass M, resistance R, and dimensions w by fas from rest

More information

Candidate Number. General Certificate of Education Advanced Level Examination June 2010

Candidate Number. General Certificate of Education Advanced Level Examination June 2010 Centre Number Surname Candidate Number For Examiner s Use Other Names Candidate Signature Examiner s Initias Genera Certificate of Education Advanced Leve Examination June 2010 Question 1 2 Mark Physics

More information

b n n=1 a n cos nx (3) n=1

b n n=1 a n cos nx (3) n=1 Fourier Anaysis The Fourier series First some terminoogy: a function f(x) is periodic if f(x ) = f(x) for a x for some, if is the smaest such number, it is caed the period of f(x). It is even if f( x)

More information

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on

More information

arxiv: v1 [physics.flu-dyn] 2 Nov 2007

arxiv: v1 [physics.flu-dyn] 2 Nov 2007 A theoretica anaysis of the resoution due to diffusion and size-dispersion of partices in deterministic atera dispacement devices arxiv:7.347v [physics.fu-dyn] 2 Nov 27 Martin Heer and Henrik Bruus MIC

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

Measurement of acceleration due to gravity (g) by a compound pendulum

Measurement of acceleration due to gravity (g) by a compound pendulum Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams Unit 48: Structura Behaviour and Detaiing for Construction 4.1 Introduction Defection of Beams This topic investigates the deformation of beams as the direct effect of that bending tendency, which affects

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Nuclear Size and Density

Nuclear Size and Density Nucear Size and Density How does the imited range of the nucear force affect the size and density of the nucei? Assume a Vecro ba mode, each having radius r, voume V = 4/3π r 3. Then the voume of the entire

More information

Radiation Fields. Lecture 12

Radiation Fields. Lecture 12 Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

2.1. Cantilever The Hooke's law

2.1. Cantilever The Hooke's law .1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

3.10 Implications of Redundancy

3.10 Implications of Redundancy 118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.

More information

Fast magnetohydrodynamic oscillations in a multifibril Cartesian prominence model. A. J. Díaz, R. Oliver, and J. L. Ballester

Fast magnetohydrodynamic oscillations in a multifibril Cartesian prominence model. A. J. Díaz, R. Oliver, and J. L. Ballester A&A 440, 1167 1175 (2005) DOI: 10.1051/0004-6361:20052759 c ESO 2005 Astronomy & Astrophysics Fast magnetohydrodynamic osciations in a mutifibri Cartesian prominence mode A. J. Díaz, R. Oiver, and J. L.

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 506 Winter 2006 Homework Assignment #6 Solutions Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by

More information

USPAS Course on Recirculated and Energy Recovered Linacs

USPAS Course on Recirculated and Energy Recovered Linacs USPAS Course on Recircuated and Energy Recovered Linacs I. V. Bazarov Corne University D.R. Dougas, G. A. Krafft, and L. Merminga Jefferson Lab Computer Cass: Linear Optics in JLAB IRFEL, Longitudina gymnastics,

More information

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis Laboratory Exercise 1: Penduum Acceeration Measurement and Prediction Laboratory Handout AME 20213: Fundamentas of Measurements and Data Anaysis Prepared by: Danie Van Ness Date exercises to be performed:

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Lecture 17 - The Secrets we have Swept Under the Rug

Lecture 17 - The Secrets we have Swept Under the Rug Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

General Certificate of Education Advanced Level Examination June 2010

General Certificate of Education Advanced Level Examination June 2010 Genera Certificate of Education Advanced Leve Examination June 2010 Human Bioogy HBI6T/P10/task Unit 6T A2 Investigative Skis Assignment Task Sheet The effect of temperature on the rate of photosynthesis

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

Technical Data for Profiles. Groove position, external dimensions and modular dimensions Technica Data for Profies Extruded Profie Symbo A Mg Si 0.5 F 25 Materia number.206.72 Status: artificiay aged Mechanica vaues (appy ony in pressing direction) Tensie strength Rm min. 245 N/mm 2 Yied point

More information

Traffic data collection

Traffic data collection Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road

More information

High-order approximations to the Mie series for electromagnetic scattering in three dimensions

High-order approximations to the Mie series for electromagnetic scattering in three dimensions Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

Theory and implementation behind: Universal surface creation - smallest unitcell

Theory and implementation behind: Universal surface creation - smallest unitcell Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is

More information

General Certificate of Education Advanced Level Examination June 2010

General Certificate of Education Advanced Level Examination June 2010 Genera Certificate of Education Advanced Leve Examination June 2010 Human Bioogy HBI6T/Q10/task Unit 6T A2 Investigative Skis Assignment Task Sheet The effect of using one or two eyes on the perception

More information

Preamble. Flow and Fluid Velocity. In this section of my lectures we will be. To do this we will use as an analogy

Preamble. Flow and Fluid Velocity. In this section of my lectures we will be. To do this we will use as an analogy Preambe Resistance Physics, 8 th Edition Custom Edition Cutne & Johnson Chapter 20.3 Pages 602-605 In this section of my ectures we wi be deveoping the concept of resistance. To do this we wi use as an

More information

18. Atmospheric scattering details

18. Atmospheric scattering details 8. Atmospheric scattering detais See Chandrasekhar for copious detais and aso Goody & Yung Chapters 7 (Mie scattering) and 8. Legendre poynomias are often convenient in scattering probems to expand the

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

should the warm BPMs in LHC be coated with a 100 micron copper layer? (question by Gerhard Schneider)

should the warm BPMs in LHC be coated with a 100 micron copper layer? (question by Gerhard Schneider) shoud the warm BPMs in LHC be coated with a micron copper ayer? (question by Gerhard Schneider) 46 BPMs per beam (6 BPMSW, 8 BPMW, 4 BPMWA, 8 BPMWB) Average beta Injection Top Horizonta beta Vertica beta

More information

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with? Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine

More information

A MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES

A MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES 6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES A MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES Sakae NAGAOKA* *Eectronic

More information

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame Work and energy method EI EI T x-axis Exercise 1 : Beam with a coupe Determine the rotation at the right support of the construction dispayed on the right, caused by the coupe T using Castigiano s nd theorem.

More information

1.2 Partial Wave Analysis

1.2 Partial Wave Analysis February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider

More information

Student Friendly Guide. Cosmic Microwave Background

Student Friendly Guide. Cosmic Microwave Background Student Friendy Guide to the Cosmic Microwave Background Robert D. Kauber www.quantumfiedtheory.info Juy 5, 05 Abstract This document provides what is hopefuy a reativey easy-to-understand introduction,

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

Candidate Number. General Certificate of Education Advanced Level Examination January 2012

Candidate Number. General Certificate of Education Advanced Level Examination January 2012 entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday

More information

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics

Two Kinds of Parabolic Equation algorithms in the Computational Electromagnetics Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa

More information

1. Measurements and error calculus

1. Measurements and error calculus EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

THINKING IN PYRAMIDS

THINKING IN PYRAMIDS ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think

More information

University of California, Berkeley Physics 7A Spring 2009 (Yury Kolomensky) SOLUTIONS TO PRACTICE PROBLEMS FOR THE FINAL EXAM

University of California, Berkeley Physics 7A Spring 2009 (Yury Kolomensky) SOLUTIONS TO PRACTICE PROBLEMS FOR THE FINAL EXAM 1 University of Caifornia, Bereey Physics 7A Spring 009 (Yury Koomensy) SOLUIONS O PRACICE PROBLEMS FOR HE FINAL EXAM Maximum score: 00 points 1. (5 points) Ice in a Gass You are riding in an eevator hoding

More information

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27 Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

Completion. is dense in H. If V is complete, then U(V) = H.

Completion. is dense in H. If V is complete, then U(V) = H. Competion Theorem 1 (Competion) If ( V V ) is any inner product space then there exists a Hibert space ( H H ) and a map U : V H such that (i) U is 1 1 (ii) U is inear (iii) UxUy H xy V for a xy V (iv)

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various

More information

Electromagnetic Waves

Electromagnetic Waves Eectromagnetic Waves Dispacement Current- It is that current that comes into existence (in addition to conduction current) whenever the eectric fied and hence the eectric fux changes with time. It is equa

More information

arxiv:quant-ph/ v3 6 Jan 1995

arxiv:quant-ph/ v3 6 Jan 1995 arxiv:quant-ph/9501001v3 6 Jan 1995 Critique of proposed imit to space time measurement, based on Wigner s cocks and mirrors L. Diósi and B. Lukács KFKI Research Institute for Partice and Nucear Physics

More information

Computer class: Linear Optics in JLAB: Longitudinal Dynamics and BBU

Computer class: Linear Optics in JLAB: Longitudinal Dynamics and BBU USPAS course on Recircuated and Energy Recovered Linacs Ivan Bazarov, Corne University Geoff Krafft and Dave Dougas, JLAB Computer cass: Linear Optics in JLAB: Longitudina Dynamics and BBU JLAB IRFEL Spreadsheet

More information

O9e Fringes of Equal Thickness

O9e Fringes of Equal Thickness Fakutät für Physik und Geowissenschaften Physikaisches Grundpraktikum O9e Fringes of Equa Thickness Tasks 1 Determine the radius of a convex ens y measuring Newton s rings using ight of a given waveength.

More information