Saddle point method in complex plane

Summary of the saddle point method [1].

The saddle point method is also known as the method of steepest descent or stationary phase. It estimates the asymptotic behavior of an integral of the following form

Eq. (1)

𝐼(𝜆)=𝐶d𝑧e𝜆𝑆(𝑧),𝜆

as 𝜆+. Here, 𝐶 is a contour in the complex 𝑧 plane and 𝑆 is a holomorphic function.

In order to evaluate this integral, we deform 𝐶 such that it passes through the saddle points of 𝑆 (where 𝑆(𝑧)=0) in the steepest descent directions. The idea of the saddle point method is that, as 𝜆+, the main contribution to the integral comes from the neighborhood of the saddle points. For simplicity, let’s assume that 𝑆 has only one saddle point 𝑧=𝑧0; if 𝑆 has multiple saddle points, we need only to sum the contributions from each. Taylor expanding [2] 𝑆 about 𝑧0 and noting that the Gaussian profile has width |𝑧𝑧0|=𝒪︀(𝜆1/2), we have

Eq. (2)

𝐼(𝜆)=𝐶d𝑧exp(𝜆𝑆(𝑧0)+12𝜆𝑆(𝑧0)(𝑧𝑧0)2+𝒪︀(𝜆(𝑧𝑧0)3)+𝒪︀(𝜆(𝑧𝑧0)4))=e𝜆𝑆(𝑧0)𝐶d𝑧exp(12𝜆𝑆(𝑧0)(𝑧𝑧0)2)×(1+𝒪︀(𝜆(𝑧𝑧0)3)+𝒪︀(𝜆(𝑧𝑧0)4)).

Now we need to determine the direction of the steepest descent. Let 𝑧𝑧0=(𝜌/𝜆)ei𝜃 and 𝑆(𝑧0)=|𝑆(𝑧0)|ei𝛼, then the real part of the exponent

Eq. (3)

Re(12𝜆𝑆(𝑧0)(𝑧𝑧0)2)=12|𝑆(𝑧0)|𝜌2Reei(2𝜃+𝛼)=12|𝑆(𝑧0)|𝜌2cos(2𝜃+𝛼)

decreases most rapidly when cos(2𝜃+𝛼)=1, i.e., when 𝜃=±π/2𝛼/2𝜃0 [3]. Choosing 𝐶 to pass through 𝑧0 along the direction determined by the proper 𝜃=𝜃0, we find

Eq. (4)

𝐼(𝜆)=e𝜆𝑆(𝑧0)ei𝜃0𝜆d𝜌exp(12|𝑆(𝑧0)|𝜌2)×(1+𝒪︀(1𝜆)𝜌3+𝒪︀(1𝜆)𝜌4)=2π𝜆|𝑆(𝑧0)|ei𝜃0+𝜆𝑆(𝑧0)(1+𝒪︀(1𝜆)).

In physical problems, one often encounter integrals of the following form

Eq. (5)

𝐼(𝜆)=𝐶d𝑧𝑓(𝑧)e𝜆𝑆(𝑧)

with 𝑓 being holomorphic. In this case, we may repeat the above derivation to obtain

Eq. (6)

𝐼(𝜆)=𝑓(𝑧0)e𝜆𝑆(𝑧0)ei𝜃0𝜆d𝜌exp(12|𝑆(𝑧0)|𝜌2)×(1+𝒪︀(1𝜆)𝜌+𝒪︀(1𝜆)𝜌2)(1+𝒪︀(1𝜆)𝜌3+𝒪︀(1𝜆)𝜌4)=2π𝜆|𝑆(𝑧0)|𝑓(𝑧0)ei𝜃0+𝜆𝑆(𝑧0)(1+𝒪︀(1𝜆)).

For more physical applications of the saddle point method, see [1] and [4].

Bibliography