Physics 504, Lecture 10 Feb. 24, Geometrical Fiber Optics. Last Latexed: February 21, 2011 at 15:11 1

Size: px
Start display at page:

Download "Physics 504, Lecture 10 Feb. 24, Geometrical Fiber Optics. Last Latexed: February 21, 2011 at 15:11 1"

Transcription

1 Last Latexe: Feruary 1 11 at 15:11 1 Physics 54 Lecture 1 Fe Geometrical Fier Optics The wave guies consiere so far containe their fiels within conucting walls ut we know from stuying total internal reflection in elementary optics that it is also possile to contain fiels y changes in the inex of refraction. n extremely important application is fier optics. fier optic cale is a silica fier with an inex of refraction which varies as a function of raius. The simplest form is a core of raius a an inex of refraction n 1 surroune y a shell of outer raius anhavinganinexof refraction n < n 1. If the angle of incience α>α c =sin 1 n /n 1 there will e total internal reflection an the light will e confine. It is more convenient in iscussing optical fiers to use the angle θ θ α α a which the light makes with the axis of the fier so the conition for total internal reflection is θ<θ max =cos 1 n /n 1. Of course this iscussion was in terms of geometrical optics suitale if the wavelength of the light is negligile compare to the geometrical istances a. Optical fiers come in two categories multimoe an single moe. For multimoe fiers a 5 µm 75 µm an the light use is generally near infrare.85 µm so geometrical optics is a reasonale approach though we shall see that interference effects are still relevant. Single moe fier is smaller a µm an we nee to treat these as waveguies. Consier the simple multimoe fier an efine = n 1 n n 1 1 n n 1 which is often aout.1. s it is small cos θ max 1 1 θ max =1 so θ max. There is an uncertainty principle etween the localization of a wave an its wave numer which limits the numer of moes that can e transmitte. We may think of this quantum mechanically where the ensity ofmoesisgiveny pq/π h D for each moe in D imensions. The 54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 coorinate integral is q = πa anas k k z tan θ max = k z p = h k =π h kz. There are polarizations so the numer of moes that can propagate is roughly N = kza = 1 V wherev := ka is calle the fier parameter. For a multimoe fier N is aout 1 ut for a single moe it is one for each polarization. There is a prolem with this simple arrangement as the istance light travels to get a istance z own the fier is z sec θ so light with ifferent θ values will travel ifferent optical istances to get to the same θ point an will interfere. This can e ameliorate y having more than one transition or even a continuum. In fact for next week you will fin how it can e fixe with a continuous istriution of nr. To analyze such a situation where ɛ x varies smoothly an assuming µ = µ as the fier is not magnetic we may write Maxwell s equations having fourier transforme time to iscuss a particular frequency as ɛe = = ɛ E + ɛ E E H = µ t = iµ ωh H = ɛ E = iωɛe t H =. So E = E + E = iµ ω H = µ ω ɛe = E 1 ɛ E ɛ H = H + H = iω ɛe or H = iω ɛ E iωɛ E = iω ɛ E + µ ω ɛ H. E + µ ω ɛe + 1 ɛ E = ɛ H + µ ω ɛh iω ɛ E =

2 54: Lecture 1 Last Latexe: Feruary 1 11 at 15: : Lecture 1 Last Latexe: Feruary 1 11 at 15:11 4 We can simplify these equations if we assume that ɛ varies slowly compare to a wavelength so ɛ ɛ = ɛω c. The other terms in the equations are of orer ω /c times E or H respectively ɛe H/c sothe ɛterms may e roppe as small. Then the components of oth E an H satisfy 1 + ω c n r ψ r =. 1 This equation which escries the rapily oscillating function ψ can e replace y a more graually varying function y using the Eikonal writing ψ r =e iωs r/c so ψ = iω Se c [ iω ω S ] iωs r/c = c S i e iωs r/c. c This is ω n /c e iωs/c from 1 so n r S S = i c ω S. c/ω soas S varies on the same scale as n r which is slowly compare to 1/ the right han sie can e roppe an we have the eikonal approximation S S = n r. This equation oesn t tell us the irection in which S changes ut it oes tell us that the rate is just n r so in following a particular ray s path we can efine a unit vector ˆk r such that S = n rˆk r. Near a point r we may expan ψ r =e iω S r + r r S /c = e iωs r /c e iωˆk r r n r/c 1 Jackson uses x instea of r in most of these equations to emphasize that n epens only on the transverse location an not on the axial coorinate z. But clearly the wave function ψ an the eikonal S o epen on z along a given ray. so it is locally a plane wave with k = ωn r/c suitale for a wave with spee c/n r as one might expect. If s is the istance measure along the path of the ray r/s = ˆk n r r/s = S so s n r r s = S s = S = n r. s 3 Γ In the penultimate expression Γ represents the ray s path. In general a ray will e irecte in a irection with a large axial component a raial component an an aximuthal component. If the latter is zero the ray will pass through the axis an the ray is calle meriional. Rays with a nonzero azimuthal component will never pass irectly through the axis an are calle skew rays. In terms of wave functions such rays correspon to azimuthal quantum numers m an have vanishing intensity at ρ = on the axis. Following Jackson we will avoi complication an iscuss only meriional rays which is also equivalent for geometric optics at least to iscussing a plane sla. So the ρ irection will e rename x an we consier a ray confine to the xz plane where z is the overall irection of motion. We assume n oesn t epen on z. The x an z components of Eq 3 are nxsinθ =nx s x nxcosθ =n r =. s z Thus nxcosθ is a constant along the path an if θ<θ max the path will turn aroun at x max atwhichcosθ = 1 so the constant value of nxcosθ is just n := nx max ut it is also of course n cos θ. Note that z/s =cosθ = n/n so the path is etermine y rewriting the x component of Eq. 3 as n n x nx = n nx z nx z nx z n x = n z x so n x z = nxnx x = 1 x n x. Whether y inspiration from mechanics viewing 1 n x as a potential or otherwise if we multiply this equation y x := x/z we get 1 n z x = 1 z z n x = n x = n x n

3 54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 5 where the constant of integration n is etermine y x =atx max nx max = n. So the axial istance travele from the point the ray crosse the axis to reaching raius x is zx = x z x x x = n x n x n. The path is important ecause if ifferent rays have ifferent path lengths to get to the same large isplacement z own the fier they will interfere. From one crossing of the axis to the next the ray moves in z aistance xmax x Z = n n x n. The istance travelle is not important ut the optical istance nxs is ecause that etermines the change in phase an the optical istance corresponing to the axial isplacement Z is L opt = xmax nx s z xmax z x x = xmax n x = n x n x. nx nx n n n x n x The phase ifference in a single half-perio is not likely to e important ut when the rays travel a large istance z the total optical path will e L opt z/z. Thus it is ieal if L opt /Z is inepenent of n that is it is the same for all rays. In prolem 8.14 you will show how to accomplish that. In aition to geometrical ispersion coming from not having L opt /Z inepenent of n there can also e frequency ispersion ue to the variation of n x with frequency. This will not effect a very narrow anwith signal one with a very narrow range of ω ut the rate at which information can e carrie is proportional to the anwith so we woul like to have little ispersion in frequency as well. We also want minimal asorption. These two issues for silica encourage using 1.4µm. We will skip Jackson : Lecture 1 Last Latexe: Feruary 1 11 at 15:11 6 Sources of Waves in Wave Guies We now turn our attention to the sources of waves in waveguies. With given sources i.e. ignoring ack reactions the equations are still linear inhomogeneous in the fiels an time-inepenent so we can assume everything has a e iωt time epenence. The normal moes of free waves in a waveguie form a complete set of states for the source-free solutions to Maxwell s equations in the guie though we nee to inclue all the moes incluing those whose cutoff frequency is aove ω. Far from the sources however only the real values of k formoeswithω <ω will have appreciale amplitue. We will expan our fiels in the normal moes which will e inexe y a composite inex which inclues information aout whether the moe is TE or TM or of other nature e.g. TEM an also the inices which specify which moe roughly the angular an raial quantum numers For a given moe we have two values of k either positive right moving an negative left moving or i k for cutoff moes. We will lump the +i k moes in with the positive k ones an write the + part of the moe as E + x y z = [ E x y+ẑe z x y ] e ik z H + x y z = [ H x y+ẑh z x y ] e ik z where E an H are purely transverse an were given in terms of E z or H z for TM or TE moes respectively earlier. There are also moes travelling in the negative z irection. The equations are symmetric uner z z uner which the transverse E x y is unchange ut E z x y changes sign. The magnetic fiel is a pseuo-vector so for it the transverse H x y changes sign ut H z x y is unchange. Thus E x y z = [ E x y ẑe z x y ] e ik z H x y z = [ H x y+ẑh z x y ] e ik z The transverse fiels E form a asis which we can choose to e real an normalize E E µ = δ µ.

4 54: Lecture 1 Last Latexe: Feruary 1 11 at 15: : Lecture 1 Last Latexe: Feruary 1 11 at 15:11 8 That this can e one an what this normalization requires for the asis of normal moes for E z an H z is elaorate in the slies. With this normalization an from H = Z 1 ẑ E we also have H H µ = 1 Z δ µ an in the time-average power expression P = 1 E H µ ẑ = 1 Z δ µ E H ẑ we have The aove normalization for the E comes from requiring that the z components satisfy orthogonality conitions an normalization conitions ajuste appropriately. For TM waves E = ik γ E z so δ µ = E E µ = k k µ E γ z E zµ = k k µ E γ µ γ z E zµ γ µ = k k µ E γ z E zµ where in the integration y parts thir = the surface term vanishes as E S = an for the fourth E zµ = γµe zµ so E z E zµ = γ δ k µ whilefortewavesthesameargumentforh gives H z H zµ = γ δ Z µ. k For a rectangular waveguie x a y the moes are laelle y integers m an n with TM waves: ψ S = TE waves: E zmn = ψ = iγ mn sin k a E xmn = πm γ mn a a cos E ymn = ψ n S = πn γ mn a sin mπx a mπx a mπx H zmn = ψ = iγ mn cos k Z a E xmn = πn γ mn a cos a E ymn = πm γ mn a a sin a a sin cos sin nπy nπy nπy mπx nπy cos a mπx nπy sin mπx nπy cos where m γmn = π a + n. The functional form of ψ is immeiately apparent from the ounary conitions an for TM moes E = ik ψ/γ anfortemoese = izkẑ ψ/γ. The overall constants are etermine from the normalization E x + Ey = 1 except that for TE moes we nee an extra factor of 1/ foreach n or m which is zero as cos mπx/a =a1 + δ m /..1 Expansion of free waves In our waveguie an aritrary fiel configuration in the asence of sources can e escrie y expaning in normal moes with positive an negative or +i an i wave numers as escrie aove. Thus a total fiel E = E + + E H = H + + H

5 54: Lecture 1 Last Latexe: Feruary 1 11 at 15: : Lecture 1 Last Latexe: Feruary 1 11 at 15:11 1 with E = E H = H can escrie an aritrary fiel configuration in a region that has no sources with the s constant inepenent of z in such a section of the wave guie. The coefficients are etermine y the transverse components E an H along any cross section for E has expansion coefficients + eikz + e ikz while H has expansion coefficients + eikz e ikz. This gives the expansion coefficients taking z = as = 1 E E Z H H.. Localize Sources We have escrie the waves that can propagate in the waveguie ut what actually prouces such waves? We will consier a localize source with current ensity J xe iωt confine to some region z [z z + ] at the ens of which we imagine cross sections enote y S S +. We will assume there are no sources or ostacles to the right of S + or to the left of S soats + there is no amplitue for any moe with negative k or with i k which woul represent left-moving waves or exponential low up as z +. The reverse is true at S so E = + E + H = + H + at S + E = E H = H at S In etween we have the full Maxwell equations with sources E = B t = iωµ H H = J + ɛ E t = J iωɛ E while the normal moes oey Maxwell equations without sources: If we apply the ientity H = iωɛ E E = iωµ H. B = B B we fin E H E H = E H E H E H + E H = iωµ H H + iωɛ E E iωµ H H + E J iωɛ E = J E If we integrate this over the volume etween S an S + using Gauss theorem an the ounary conition that E = at the conuctor s surface E H E H ˆn = J E S where S consists of S + with ˆn =ẑ ans with ˆn = ẑ. Let s take the upper sign. The contriution from S + is can e reuce to an integral over at z =: E + H + E + H + + S + = + E H E H S + z eik +k z = + E Z 1 ẑ E E Z 1 ẑ E = 1 + E Z E 1 E Z E e ik +k z = 1 + δ 1 e ik +k z =. Z Z On the other han the contriution from S is E H S + E + H ẑ = S E H + E H ẑe ik k z so = Z δ = Z = Z V z J E +. V z eik +k z

6 54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 11 The same argument for the lower sign as spelle out in the ook gives the equation with the superscript signs reverse so oth are = Z V J E. In aition to sources ue to currents we may have contriutions ue to ostacles or holes in the conucting ounaries. These can e treate as aitional surface terms in Gauss law y excluing ostacles from the region of integration V ut this requires knowing the full fiels at the surface of the ostacles or the missing parts of the waveguie conuctor. This is treate in 9.5B ut we won t iscuss it here. So finally we are at the en of Chapter 8.

Fiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:

Fiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers: Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total

More information

Physics 504, Lecture 7 Feb 14, Energy Flow, Density and Attenuation. 1.1 Energy flow and density. Last Latexed: February 11, 2011 at 10:54 1

Physics 504, Lecture 7 Feb 14, Energy Flow, Density and Attenuation. 1.1 Energy flow and density. Last Latexed: February 11, 2011 at 10:54 1 Last Latexe: February, at :54 Physics 54, Lecture 7 Feb 4, Energy Flow, Density an ttenuation We have seen that there are iscrete moes for electromagnetic waves with E, B e ikz it corresponing, for real

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2) 1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like)

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like) Stockton Univeristy Chemistry Program, School of Natural Sciences an Mathematics 101 Vera King Farris Dr, Galloway, NJ CHEM 340: Physical Chemistry II Solving the Schröinger Equation for the 1 Electron

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Ph12b Solution Set 4

Ph12b Solution Set 4 Ph Solution Set 4 Feruary 5, 00 The state ψ S un E is transforme to p ψ S un E + p a 0 S γ E + S η E Since γ E an η E aren t mutually orthogonal, we can t use this equation irectly to compute the new ensity

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

PHY 114 Summer 2009 Final Exam Solutions

PHY 114 Summer 2009 Final Exam Solutions PHY 4 Summer 009 Final Exam Solutions Conceptual Question : A spherical rubber balloon has a charge uniformly istribute over its surface As the balloon is inflate, how oes the electric fiel E vary (a)

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses.

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses. Physics 195a Course Notes The K 0 : An Interesting Example of a Two-State System 021029 F. Porter 1 Introuction An example of a two-state system is consiere. involve: 1. Treatment of a ecaying particle.

More information

The Three-dimensional Schödinger Equation

The Three-dimensional Schödinger Equation The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Moving Charges And Magnetism

Moving Charges And Magnetism AIND SINGH ACADEMY Moving Charges An Magnetism Solution of NCET Exercise Q -.: A circular coil of wire consisting of turns, each of raius 8. cm carries a current of. A. What is the magnitue of the magnetic

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Physics 2212 GJ Quiz #4 Solutions Fall 2015

Physics 2212 GJ Quiz #4 Solutions Fall 2015 Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

v r 1 E β ; v r v r 2 , t t 2 , t t 1 , t 1 1 v 2 v (3) 2 ; v χ αβγδ r 3 dt 3 , t t 3 ) βγδ [ R 3 ] exp +i ω 3 [ ] τ 1 exp i k v [ ] χ αβγ , τ 1 dτ 3

v r 1 E β ; v r v r 2 , t t 2 , t t 1 , t 1 1 v 2 v (3) 2 ; v χ αβγδ r 3 dt 3 , t t 3 ) βγδ [ R 3 ] exp +i ω 3 [ ] τ 1 exp i k v [ ] χ αβγ , τ 1 dτ 3 ON CLASSICAL ELECTROMAGNETIC FIELDS PAGE 58 VII. NONLINEAR OPTICS -- CLASSICAL PICTURE: AN EXTENDED PHENOMENOLOGICAL MODEL OF POLARIZATION : As an introuction to the subject of nonlinear optical phenomena,

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

1 Boas, p. 643, problem (b)

1 Boas, p. 643, problem (b) Physics 6C Solutions to Homework Set #6 Fall Boas, p. 643, problem 3.5-3b Fin the steay-state temperature istribution in a soli cyliner of height H an raius a if the top an curve surfaces are hel at an

More information

Lect. 4 Waveguides (1)

Lect. 4 Waveguides (1) Lect. 4 Waveguies (1) - Waveguie: Confines an guies EM waves Metallic, Dielectric, Plasmonic - We are intereste in ielectric waveguie Total internal reflection b refractive inex ifferences n A B C D n

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Experiment 2, Physics 2BL

Experiment 2, Physics 2BL Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor

More information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

The Principle of Least Action and Designing Fiber Optics

The Principle of Least Action and Designing Fiber Optics University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

Solution Set #7

Solution Set #7 05-455-0073 Solution Set #7. Consier monochromatic light of wavelength λ 0 incient on a single slit of with along the -ais an infinite length along y. The light is oserve in the Fraunhofer iffraction region

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N Physics 15 Solution Set Winter 018 1. At time t = 0, the wave function of a free particle moving in a one-imension is given by, ψ(x,0) = N where N an k 0 are real positive constants. + e k /k 0 e ikx k,

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE

Prof. Dr. Ibraheem Nasser electric_charhe 9/22/2017 ELECTRIC CHARGE ELECTRIC CHARGE Introuction: Orinary matter consists of atoms. Each atom consists of a nucleus, consisting of protons an neutrons, surroune by a number of electrons. In electricity, the electric charge

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Get Solution of These Packages & Learn by Video Tutorials on WAVE OPTICS

Get Solution of These Packages & Learn by Video Tutorials on   WAVE OPTICS FREE Downloa Stuy Package from website: wwwtekoclassescom & wwwmathsbysuhagcom Get Solution of These Packages & Learn by Vieo Tutorials on wwwmathsbysuhagcom PRINCIPLE OF SUPERPOSITION When two or more

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Determinant and Trace

Determinant and Trace Determinant an Trace Area an mappings from the plane to itself: Recall that in the last set of notes we foun a linear mapping to take the unit square S = {, y } to any parallelogram P with one corner at

More information

Answers to Coursebook questions Chapter 5.6

Answers to Coursebook questions Chapter 5.6 Answers to Courseook questions Chapter 56 Questions marke with a star (*) use the formula for the magnetic fiel create y a current μi ( = ) which is not on the syllaus an so is not eaminale See Figure

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

Supplementary Information Parametric down-conversion photon-pair source on a nanophotonic chip

Supplementary Information Parametric down-conversion photon-pair source on a nanophotonic chip Supplementary Information Parametric own-conversion photon-pair source on a nanophotonic chip Xiang Guo, Chang-Ling Zou, Carsten Schuck, Hojoong Jung, Risheng Cheng, an Hong Tang Department of Electrical

More information

Physics 2212 K Quiz #2 Solutions Summer 2016

Physics 2212 K Quiz #2 Solutions Summer 2016 Physics 1 K Quiz # Solutions Summer 016 I. (18 points) A positron has the same mass as an electron, but has opposite charge. Consier a positron an an electron at rest, separate by a istance = 1.0 nm. What

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

The Press-Schechter mass function

The Press-Schechter mass function The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

Joel A. Shapiro January 21, 2010

Joel A. Shapiro January 21, 2010 Joel A. shapiro@physics.rutgers.edu January 21, 20 rmation Instructor: Joel Serin 325 5-5500 X 3886, shapiro@physics Book: Jackson: Classical Electrodynamics (3rd Ed.) Web home page: www.physics.rutgers.edu/grad/504

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

De Broglie s Pilot Waves

De Broglie s Pilot Waves De Broglie s Pilot Waves Bohr s Moel of the Hyrogen tom: One way to arrive at Bohr s hypothesis is to think of the electron not as a particle but as a staning wave at raius r aroun the proton. Thus, nλ

More information

A General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication

A General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication A General Analytical Approximation to Impulse Response of 3- Microfluiic Channels in Molecular Communication Fatih inç, Stuent Member, IEEE, Bayram Cevet Akeniz, Stuent Member, IEEE, Ali Emre Pusane, Member,

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Mathematical Methods wks 5,6: PDEs

Mathematical Methods wks 5,6: PDEs Mathematical Methos wks 5,6: PDEs John Magorrian, magog@thphys.ox.ac.uk These are work-in-progress notes for the secon-year course on mathematical methos. The most up-to-ate version is available from http://www-thphys.physics.ox.ac.uk/people/johnmagorrian/mm.

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information

Revisiting Fresnel & refractive index. What is the refractive index of a dielectric. Metals and plasmons

Revisiting Fresnel & refractive index. What is the refractive index of a dielectric. Metals and plasmons Revisiting Fresnel & refractive ine What is the refractive ine of a ielectric Metals an plasmons Squeezing plasmons in a nanowire Moe with 150 nm SPP l < 1 mm At l 1.550 mm Snell s law Generic solution

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Physics for Scientists & Engineers 2

Physics for Scientists & Engineers 2 Capacitors Physics for Scientists & Engineers 2 Spring Semester 2005 Lecture 12 Capacitors are evices that can store electrical energy Capacitors are use in many every-ay applications Heart efibrillators

More information

Appendix: Proof of Spatial Derivative of Clear Raindrop

Appendix: Proof of Spatial Derivative of Clear Raindrop Appenix: Proof of Spatial erivative of Clear Rainrop Shaoi You Robby T. Tan The University of Tokyo {yous,rei,ki}@cvl.iis.u-tokyo.ac.jp Rei Kawakami Katsushi Ikeuchi Utrecht University R.T.Tan@uu.nl Layout

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ Department of Physics University of arylan College Park, arylan PHYSICS 4 Prof. S. J. Gates Fall 5 Final Exam Dec. 6, 5 This is a OPEN book examination. Rea the entire examination before you begin to work.

More information

SOLUTIONS for Homework #3

SOLUTIONS for Homework #3 SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

ECE341 Test 2 Your Name: Tue 11/20/2018

ECE341 Test 2 Your Name: Tue 11/20/2018 ECE341 Test Your Name: Tue 11/0/018 Problem 1 (1 The center of a soli ielectric sphere with raius R is at the origin of the coorinate. The ielectric constant of the sphere is. The sphere is homogeneously

More information

Electromagnetic surface and line sources under coordinate transformations

Electromagnetic surface and line sources under coordinate transformations PHYSICAL REVIEW A 8, 3382 29 Electromagnetic surface an line sources uner coorinate transformations Steven A. Cummer,* Nathan Kuntz, an Bogan-Ioan Popa Department of Electrical an Computer Engineering

More information

Exam 2 Review Solutions

Exam 2 Review Solutions Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information