Fall 2004/05 Solutions to Assignment 5: The Stationary Phase Method Provided by Mustafa Sabri Kilic. I(x) = e ixt e it5 /5 dt (1) Z J(λ) =

Similar documents
ASTR 3740 Relativity & Cosmology Spring Answers to Problem Set 4.

Chapter 19 Webassign Help Problems

Then the number of elements of S of weight n is exactly the number of compositions of n into k parts.

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Gravity. David Barwacz 7778 Thornapple Bayou SE, Grand Rapids, MI David Barwacz 12/03/2003

Solutions Practice Test PHYS 211 Exam 2

Quantum Mechanics II

3. Perturbation of Kerr BH

Magnetic Dipoles Challenge Problem Solutions

On a quantity that is analogous to potential and a theorem that relates to it

Estimation and Confidence Intervals: Additional Topics

MATH 417 Homework 3 Instructor: D. Cabrera Due June 30. sin θ v x = v r cos θ v θ r. (b) Then use the Cauchy-Riemann equations in polar coordinates

Chapter 2: Basic Physics and Math Supplements

On the integration of the equations of hydrodynamics

Math 2263 Solutions for Spring 2003 Final Exam

In the previous section we considered problems where the

Magnetic Field. Conference 6. Physics 102 General Physics II

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

CHAPTER 25 ELECTRIC POTENTIAL

Math 124B February 02, 2012

A New Approach to General Relativity

How can you find the dimensions of a square or a circle when you are given its area? When you multiply a number by itself, you square the number.

Riemann s Functional Equation is Not Valid and its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.

3.6 Applied Optimization

Numerical Integration

(a) Unde zeo-bias conditions, thee ae no lled states on one side of the junction which ae at the same enegy as the empty allowed states on the othe si

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 18

Introduction to Laplace Transform Techniques in Circuit Analysis

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13

k. s k=1 Part of the significance of the Riemann zeta-function stems from Theorem 9.2. If s > 1 then 1 p s

Tutorial on Strehl ratio, wavefront power series expansion, Zernike polynomials expansion in small aberrated optical systems By Sheng Yuan

radians). Figure 2.1 Figure 2.2 (a) quadrant I angle (b) quadrant II angle is in standard position Terminal side Terminal side Terminal side

Physics 505 Homework No. 9 Solutions S9-1

Impulse and Momentum

So, if we are finding the amount of work done over a non-conservative vector field F r, we do that long ur r b ur =

MODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...

to point uphill and to be equal to its maximum value, in which case f s, max = μsfn

Physics 4A Chapter 8: Dynamics II Motion in a Plane

Inference for A One Way Factorial Experiment. By Ed Stanek and Elaine Puleo

Chapters 5-8. Dynamics: Applying Newton s Laws

10.2 Parametric Calculus

Force & Motion: Newton s Laws

Section 26 The Laws of Rotational Motion

KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS

Basic propositional and. The fundamentals of deduction

Question 1: The dipole

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

Mechanics and Special Relativity (MAPH10030) Assignment 3

Histogram Processing

Practice Integration Math 120 Calculus I Fall 2015

Image Enhancement: Histogram-based methods

Practice Integration Math 120 Calculus I D Joyce, Fall 2013

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form

( ) F α. a. Sketch! r as a function of r for fixed θ. For the sketch, assume that θ is roughly the same ( )

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

Riemann s Functional Equation is Not a Valid Function and Its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.

LECTURE 14. m 1 m 2 b) Based on the second law of Newton Figure 1 similarly F21 m2 c) Based on the third law of Newton F 12

EN40: Dynamics and Vibrations. Midterm Examination Thursday March

Non-Linear Dynamics Homework Solutions Week 2

C/CS/Phys C191 Shor s order (period) finding algorithm and factoring 11/12/14 Fall 2014 Lecture 22

Homework 7 Solutions

6 Matrix Concentration Bounds

of the contestants play as Falco, and 1 6

13. Adiabatic Invariants and Action-Angle Variables Michael Fowler

Lecture 8 - Gauss s Law

Laplace Transformation

Auchmuty High School Mathematics Department Advanced Higher Notes Teacher Version

Simulation of Spatially Correlated Large-Scale Parameters and Obtaining Model Parameters from Measurements

Electromagnetic Theory 1

A Hartree-Fock Example Using Helium

Lecture 3. January 9, 2018

Conditional Convergence of Infinite Products

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Suppose the medium is not homogeneous (gravity waves impinging on a beach,

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)

15 Solving the Laplace equation by Fourier method

Electric Potential and Gauss s Law, Configuration Energy Challenge Problem Solutions

c n ψ n (r)e ient/ h (2) where E n = 1 mc 2 α 2 Z 2 ψ(r) = c n ψ n (r) = c n = ψn(r)ψ(r)d 3 x e 2r/a0 1 πa e 3r/a0 r 2 dr c 1 2 = 2 9 /3 6 = 0.

Precision Spectrophotometry

Σk=1. g r 3/2 z. 2 3-z. g 3 ( 3/2 ) g r 2. = 1 r = 0. () z = ( a ) + Σ. c n () a = ( a) 3-z -a. 3-z. z - + Σ. z 3, 5, 7, z ! = !

, the tangent line is an approximation of the curve (and easier to deal with than the curve).

Lecture 04: HFK Propagation Physical Optics II (Optical Sciences 330) (Updated: Friday, April 29, 2005, 8:05 PM) W.J. Dallas

HRW 7e Chapter 13 Page 1 of 5

Moment-free numerical approximation of highly oscillatory integrals with stationary points

A Tutorial on Multiple Integrals (for Natural Sciences / Computer Sciences Tripos Part IA Maths)

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

DOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS

Green s Identities and Green s Functions

Do not turn over until you are told to do so by the Invigilator.

The second law of thermodynamics - II.

Boundary Layers and Singular Perturbation Lectures 16 and 17 Boundary Layers and Singular Perturbation. x% 0 Ω0æ + Kx% 1 Ω0æ + ` : 0. (9.

B da = 0. Q E da = ε. E da = E dv

Chem 453/544 Fall /08/03. Exam #1 Solutions

Honors Classical Physics I

Chapter 3 Optical Systems with Annular Pupils

2 Governing Equations

3.1 Random variables

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Tutorial Exercises: Central Forces

Transcription:

8.35 Fall 24/5 Solution to Aignment 5: The Stationay Phae Method Povided by Mutafa Sabi Kilic. Find the leading tem fo each of the integal below fo λ >>. (a) R eiλt3 dt (b) R e iλt2 dt (c) R eiλ co t dt 2. Find the leading tem fo each of the integal below λ >>. I(x) = Conide both cae in which x> o x<. Solution:. (a) I(λ) = R eiλt3 dt Can wite I(λ) = Z + e ixt e it5 /5 dt () Z We obeve that the econd integal i the complex conjugate of the fit, hence I(λ) =2Re Z Moeove, we ee that the point t =i the only tationay phae point, which give the main contibution to the integal, theefoe we can make the appoximation whee I(λ) 2Re J(λ) = =2ReJ(λ) Note that we do not need to calculate the end point contibution, which would come out to be of mallle ode than the contibution of the tationay phae point-we ae only aked the leading tem. The integand of J i ocillatoy. We would like to change the contou into one on which i exponentially deceaing, fo then we we would be able to expe it a a gamma function. Thu we put z 3 = it 3 whee t ieal.thimeanthepathovewhichtheintegalitakenwouldbe z = i /3 t

Thee ae thee oot fo i /3 : e i/6,e 5i/6,e 3i/2. An analyi of thoe path bing that we can only change the domain of integal, which i oiginally the eal axi, into the taightline making an angle /6, which i the only cae whee we can cloe ou contou with a zeo contibution ac. Thu J = e i/6 e λt3 dt Makingonemoevaiablechange we aive at J = ei/6 3λ /3 Theefoe, the anwe i given by. (b) R e iλx2 dx τ = λt 3 e τ τ 2/3 dτ = ei/6 Γ(/3) /3 3λ I(λ) Γ(/3) /3 3λ Thee i no tationay phae point in the domain of the integal, hence thi i an integal in the fom (8.36) of the book, i.e of the fom I(λ) = Z b a e iλu(x) h(x)dx (2) fo which the end point contibution ae impotant. Hence the leading fom can be given by (8.39) in the book, which i e iλu(b) h(b) iλu (b) eiλu(a) h(a) iλu (a) In ou cae, we have only one end point, hence the leading tem i eiλu() h() iλu () = eiλ 2iλ. (c) R eiλ co x dx The tationay phae point ae olution to in x =, which ae x =,. Since tationay phae point contibute moe than the end point, we do not conide end point fo the pupoe of calculating the leading tem. In thi example, the end point ae tationay phae point, o they will aleady be taken cae of when calculating the tationay phae point, o in ome ene, we can ay that thee i no end point contibution. We ue the fomula (8.45) and (8.46) fom the book which ae e i/4 2 λu (x ) eiλu(x ) h(x ) and e i/4 2 λ u (x ) eiλu(x ) h(x ) (3) 2

depending on the ign of u (x ). Since the tationay phae point ae end point, they contibute only half of the quantity they would if they wee inteio point. Thu we finally find 2 2 e i/4 λ eiλ + 2 2 2 ei/4 λ e iλ = λ co(λ 4 ) a the leading tem. 2. The integal () i in the fom (2), with u(t) =t, h(t) =e it5 /5. The integal ha neithe point of tationay phae no finite end point. We fit teat the imple cae x<. We tat with caling the vaiable of integation by the tanfomation t =( x) /4 z, and the integal () become I(x) =( x) /4 e iλf(z) dz with f(z) = z + 5 z5 and Λ =( x) 5/4. Thi lat integand ha two tationay point z =,, which ae found by olving f (z) = +z 4 =. Theefoe the integal can be appoximated by making ue of (3), which give e i/4 2Λ e4iλ/5 + e i/4 2 e 4iΛ/5 theefoe 2 I(x) = ( x) /4 Λ co(4 5 Λ 4 ) = ( x) /4 2 ( x) co(4 5/4 5 ( x)5/4 4 )= 2( x) 3/8 co( 4 5 ( x)5/4 4 ) Fo the cae x>, we cale the oiginal integal with t = x /4 z, and the integal () become I(x) =x /4 e iλf(z) dz (4) with f(z) =z + 5 z5 and Λ = x 5/4. Thi integal till doe not have any tationay phae point, hence we look at the citical point of f by olving f (z) =+z 4 = which give z = e i/4,e 3i/4,e 5i/4,e 7i/4. Putting z = e iθ in f(z) =z + 5 z5 = (co θ + i in θ)+ 5 5 (co 5θ + i in 5θ). So the integand of (4) blow up fo lage in the egion whee in 5θ > and become exponentially mall in the egion whee in 5θ >. Thee ae hown in the figue below along with the citical point of f. The egion in which the integand become exponentially lage ae haded. We obeve that we can defom the domain of ou integal a hown in the figue, to upwad, o that the path of the integal now contain two of the citical point of f. We cannot defom the path of ou integal 3

Figue 4

downwad, a the black egion, ove which the integand become exponentially lagepeventufomdoingo. We calculate the elevant value which mean we have the expanion f(e i/4 ) = 4 5 ei/4,f (e i/4 )=4e 3i/4 f(e 3i/4 ) = 4 5 e3i/4,f (e 3i/4 )=4e i/4 iλf(z) iλ 4 5 ei/4 2Λe i/4 (z e i/4 ) 2 iλf(z) iλ 4 5 e3i/4 2Λe i/4 (z e 3i/4 ) 2 So the contibution fom thoe citical point ae calculated to be x /4 2Λe exp[iλ4 i/4 5 ei/4 ]+x /4 2Λe exp[iλ4 i/4 5 e3i/4 ] which i 2x 3/8 2 exp( 2 2 5 x5/4 )co( 2 2 5 x5/4 8 ) 5