Homework 09

1 Stirling approximation

Show using steepest descent that the leading asymptotics of Γ(z) for z=reiθ,r,Rez>0, are given by the complex Stirling approximation

Γ(z)ezzz(2πz)1/2.

The Gamma function

Γ(z)=0tz1etdt

after performing a change of variables by setting t=zu, which yields:

Γ(z)=zz0uz1ezudu

The method of steepest descent involves analyzing the integrand’s behavior around its saddle points, where the first derivative of the exponent vanishes. For our integral, the exponent is zu+(z1)log(u). The saddle point occurs where its derivative with respect to u is zero. This gives us:

z+z1u=0

u=z1z1for large z

Near the saddle point u=1, we approximate the exponent by a second-order Taylor expansion:

zu+(z1)log(u)z+(z1)log(z)z2(u1)2

The integral can now be approximated as:

Γ(z)zzeze(z1)log(z)ez2(u1)2du

With the Gaussian integral evaluates to:

ez2(u1)2du=2πz

We get:

Γ(z)ezzz(2πz)12

2 Airy function

Consider the time-independent Schrodinger equation for a linear potential:

d2dx2ψ(x)=xψ(x)

  • Show that the “Airy function” defined by

    Ai(x)=12πeixy+iy3/3dy

    solves the equation above.

  • Determine the saddle points of ixy+iy3/3 in the y plane and make a change of variables to the integral so that these saddle points are not moving. Determine the direction of steepest descent on the constant phase contours through these saddle points for each argument of complex x.

  • Find the steepest descent contour for Ai(x) for xeiα for each α. Determine the leading asymptotics for each α.

2.1 Verifying the Airy Function as a Solution

Taking the Fourier transform on both sides of equation we get

F(d2ydx2)=F(xy),

For the L.H.S., we get

F(d2ydx2)=12π+d2ydx2eikxdx==k2F(y).

For the R.H.S. based on the derivative of the Fourier transform with respect to k we get,

ddkF(y)=ddk(12π+yeikxdx),

Therefore,

ddkF(y)=12π+ddk(yeikx)dx=12π+(ix)yeikxdx

We have

F(xy)=idF(y)dk,

So the equation becomes

k2F(y)=idF(y)dk,

the only solution is

F(y)=eik33,

taking the inverse Fourier transform, we get a solution satisfying the differential equation. Hence we have the final solution as

y=Ai(x)=12π+ei(kx+k33)dk.

2.2 Saddle Points

The exponent in the integral is ixy+i3y3. To find the saddle points, we set the derivative of this expression with respect to y equal to zero:

ix+iy2=0

3 Z(ϵ) function

Compute the all-orders asymptotic series for

Z(ϵ)=dxex2ϵx4ϵ0+

and then use Borel summation to evaluate the integral.

The expansion is based on the Taylor series of eϵx4 around ϵ=0:

eϵx4=n=0(ϵx4)nn!

Substituting this series into the integral and integrating term by term, we get:

Z(ϵ)=n=0(ϵ)nn!dxx4nex2

Each term in this series can be integrated using standard methods, often involving gamma functions for even powers of x.

The Borel sum is then given by an integral of the Borel transform against an exponential:

Z~(ϵ)=0etB(tϵ)dt

where B(t) is the Borel transform of the series.