Asymptotic Expansion of I[Graphics:Images/file_gr_1.gif]

[Graphics:Images/file_gr_2.gif]

This figure shows the original path of integration (C1) from -∞ to ∞ on the real axis. If we try to split the ∫ into 2 parts ( [Graphics:Images/file_gr_3.gif]and 1) we will get ∞-∞ which is useless.

The first idea is to deform the contour of the interation from C1 (Above) to C3. In our case we will take [Graphics:Images/file_gr_4.gif] Since we don't have any poles enclosed we have:

[Graphics:Images/file_gr_5.gif]

Now Let's look at the integral on C2 and C4:

[Graphics:Images/file_gr_6.gif]
[Graphics:Images/file_gr_7.gif]

Which will ->0 as R->∞.

[Graphics:Images/file_gr_8.gif]
[Graphics:Images/file_gr_9.gif]
[Graphics:Images/file_gr_10.gif]

This now can be safely split into 2 parts the second of which is:

[Graphics:Images/file_gr_11.gif]
[Graphics:Images/file_gr_12.gif]
[Graphics:Images/file_gr_13.gif]

Ie. this integral has a finite value and infact is equal to ZERO (According to Mathematica; even if it is not zero it is a finite value and can easily be approximated).

Now we are left with the other part of the integral which is:

[Graphics:Images/file_gr_14.gif]

and

[Graphics:Images/file_gr_15.gif]

Now the Saddle Points are:

[Graphics:Images/file_gr_16.gif]
[Graphics:Images/file_gr_17.gif]

Lets Take the substitution z=τ [Graphics:Images/file_gr_18.gif]

[Graphics:Images/file_gr_19.gif]
[Graphics:Images/file_gr_20.gif]

Now take [Graphics:Images/file_gr_21.gif]

[Graphics:Images/file_gr_22.gif]
[Graphics:Images/file_gr_23.gif]
[Graphics:Images/file_gr_24.gif]

And now if we solve for the stationary points we get

[Graphics:Images/file_gr_25.gif]
[Graphics:Images/file_gr_26.gif]

These points are on the unit circule and are:

[Graphics:Images/file_gr_27.gif]

[Graphics:Images/file_gr_28.gif]

Let's evaluate the function at these 3 points

[Graphics:Images/file_gr_29.gif]
[Graphics:Images/file_gr_30.gif]

This means that there is no contribution from the bottom point (-I) and the contribution comes from [Graphics:Images/file_gr_31.gif]. The next plot shows the path of the original integration as well as the path that we will have to take to go through in order to go though the Saddel point.

[Graphics:Images/file_gr_32.gif]

[Graphics:Images/file_gr_33.gif]

[Graphics:Images/file_gr_34.gif]

We have to ensure that the integration on CM1 and CM2 Vanishes as we go to infinity so that we can deform our contour from the parbola to a strait line that goes to ∞. We will look at CM1 and the argument will be identical for CM2. Take [Graphics:Images/file_gr_35.gif] the integration will become:

[Graphics:Images/file_gr_36.gif]
[Graphics:Images/file_gr_37.gif]
[Graphics:Images/file_gr_38.gif]

Now the only thing left is to expand the orignial integration around the saddle points. We will do the saddle point [Graphics:Images/file_gr_39.gif]and the other one will be it's complex conjogate. To do this expansion we will have to make the substitution τ->R [Graphics:Images/file_gr_40.gif]

[Graphics:Images/file_gr_41.gif]
[Graphics:Images/file_gr_42.gif]
[Graphics:Images/file_gr_43.gif]

If we plot the abs value of integrand vs R with η=1we get

[Graphics:Images/file_gr_44.gif]

[Graphics:Images/file_gr_45.gif]

[Graphics:Images/file_gr_46.gif]

Which shows that we indeed have a maximum @ R=1

A Similar plot will show that the phase of the integrand stops CHANGING at this point also.  
So, Now let us just expand f[R [Graphics:Images/file_gr_47.gif]] about the point R=1. We will just take the firs 3 terms.

[Graphics:Images/file_gr_48.gif]
[Graphics:Images/file_gr_49.gif]

Ofcource we don't see anything of the form [Graphics:Images/file_gr_50.gif]since this is a stationarry point. Also if we take m=0.944941-1.63669 I then we can write the integral as:

[Graphics:Images/file_gr_51.gif]

Now Remember we have:

[Graphics:Images/file_gr_52.gif]
[Graphics:Images/file_gr_53.gif]

Which can be shown by a simple change of variables [Graphics:Images/file_gr_54.gif]gassian one. Note that the limits of the integration were changed from 0 to ∞ to -∞ to +∞ since we are just interested in the behaiviour of the function at R=1 and it is exponentially decaying away from that point. That can be trasfered to a Γ function though an other strait forward transformation. We will eventually get the same result which is:

[Graphics:Images/file_gr_55.gif]
[Graphics:Images/file_gr_56.gif]

And in terms of λ we get

[Graphics:Images/file_gr_57.gif]

And we finally multiply by 2 to take the effect of the other significant stationary point in the second quadrant.

This would be the first term of the expantion. If we wanted to find the higher order terms we would proceed to expand the exponent to get powers of [Graphics:Images/file_gr_58.gif],[Graphics:Images/file_gr_59.gif]. We then would expand and multiply the terms of the form [Graphics:Images/file_gr_60.gif]to obtain a polynomial multiplied by [Graphics:Images/file_gr_61.gif]again similar substitutions would trasfer it to a sum of Γ Integrations that can be calculated to arbitrary powers of η and consequently λ.

END OF PROBLEM


Converted by Mathematica      November 26, 2001