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

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

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

Separation of Variables

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

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

6 Wave equation in spherical polar coordinates

Homework 7 Due 18 November at 6:00 pm

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

Vectors in two dimensions

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

Quantum Mechanics in Three Dimensions

Chapter 4. Electrostatics of Macroscopic Media

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

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

Schrödinger s equation.

Ph12b Solution Set 4

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

1 Heisenberg Representation

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

PHY 114 Summer 2009 Final Exam Solutions

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. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

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

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

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

The Three-dimensional Schödinger Equation

Table of Common Derivatives By David Abraham

Qubit channels that achieve capacity with two states

Moving Charges And Magnetism

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

Physics 2212 GJ Quiz #4 Solutions Fall 2015

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

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

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

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

Physics 115C Homework 4

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

12.11 Laplace s Equation in Cylindrical and

1 Boas, p. 643, problem (b)

Lect. 4 Waveguides (1)

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

PARALLEL-PLATE CAPACITATOR

The Principle of Least Action

Experiment 2, Physics 2BL

5-4 Electrostatic Boundary Value Problems

The Principle of Least Action and Designing Fiber Optics

PDE Notes, Lecture #11

Solution Set #7

Short Intro to Coordinate Transformation

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

Conservation Laws. Chapter Conservation of Energy

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

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

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

Linear First-Order Equations

Determinant and Trace

Answers to Coursebook questions Chapter 5.6

PHYS 414 Problem Set 2: Turtles all the way down

6. Friction and viscosity in gasses

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

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

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

Calculus of Variations

Tutorial Test 5 2D welding robot

4. Important theorems in quantum mechanics

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

Physics 2212 K Quiz #2 Solutions Summer 2016

The total derivative. Chapter Lagrangian and Eulerian approaches

The Press-Schechter mass function

A Model of Electron-Positron Pair Formation

Joel A. Shapiro January 21, 2010

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

De Broglie s Pilot Waves

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

05 The Continuum Limit and the Wave Equation

Partial Differential Equations

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

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

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

The Exact Form and General Integrating Factors

Mathematical Methods wks 5,6: PDEs

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 1B, lecture 8: Integration by parts

Generalization of the persistent random walk to dimensions greater than 1

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

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

ELECTRON DIFFRACTION

Sturm-Liouville Theory

Physics for Scientists & Engineers 2

Appendix: Proof of Spatial Derivative of Clear Raindrop

The Sokhotski-Plemelj Formula

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 θ

SOLUTIONS for Homework #3

Applications of First Order Equations

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

Electromagnetic surface and line sources under coordinate transformations

Exam 2 Review Solutions

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

Introduction to the Vlasov-Poisson system

Euler equations for multiple integrals

Transcription:

Last Latexe: Feruary 1 11 at 15:11 1 Physics 54 Lecture 1 Fe. 4 11 1 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 =

54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 3 54: 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

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 8.11. 54: 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 µ = δ µ.

54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 7 54: 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

54: Lecture 1 Last Latexe: Feruary 1 11 at 15:11 9 54: 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

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.