# How can stationary phase approx lead to log terms?

In these lecture notes (http://www.math.ucla.edu/~tao/247b.1.07w/notes8.pdf, section 3), Terence Tao approximates the following integral for large $\lambda$ using the stationary phase method:

\begin{align} I_{a,\phi}(\lambda) &:= \int_{\mathbf{R}^d}\ a(x) e^{i \lambda \phi(x) } dx\\ &= \sum_{n=0}^N c_n \lambda^{-n - \frac{d}{2}} e^{i\lambda \phi(x_0)} + O_{N,a,d,\phi}(\lambda^{-N -\frac{d}{2}-1}) \end{align}

Here, $\phi$ has a single non-degenerate critical point. That much is fine. My question is about the following comment (pg 11):

The situation gets significantly more complicated when the det $\nabla^2 \phi$ vanishes; for instance, factors of log $\lambda$ begin to appear in the asymptotic expansion.

-What's a simple example of an integral that has such log $\lambda$ terms?

-Is there a heuristic way of finding the leading log $\lambda$ term (the way the leading term of the usual stationary phase approximation comes from a Gaussian integral)?

Here is an example:

$I(\lambda) = \int_0^1 dx \int_0^x dy\ e^{i \lambda(x^2 - y^2) }$

Note that $\phi(x,y) = x^2 - y^2$ has a single stationary point, $(x,y)= (0,0)$, and that $\text{det}\ \nabla^2 \phi =0$.

Direct evaluation

With some help from Mathematica, we can evaluate this integral in terms of a generalized hypergeometric function and then take $\lambda$ to be large:

\begin{align} I(\lambda) &=\frac{1}{2}\ _2 F_2(1,1;\frac{3}{2}, 2; i \lambda)\\ &\to \frac{i}{4}\frac{\ln \lambda}{\lambda} +O\left( \frac{1}{\lambda} \right) \qquad \left(\lambda \to \infty\right) \end{align}

Stationary phase approximation

The method used below is described in chapter 8 of "Asymptotic Approximations of Integrals", by R. Wong. This will be a heuristic version to leading order, adapted to this example, of what he does rigorously to all orders in a more general setting.

Setting $h(t) = \int_0^1 dx \int_0^x dy\ \delta(x^2-y^2 -t)$, we have:

\begin{align} I(\lambda) = \int_0^1 e^{i \lambda t} h(t) \end{align}

One finds: \begin{align} h(t) = -\frac{1}{4} \ln t + (\text{analytic in } t) \end{align}

Then:

\begin{align} I(\lambda) \sim \int_0^\infty dt\ e^{i \lambda t} \left( - \frac{1}{4} \ln t \right) = \frac{i}{4} \frac{\ln \lambda}{\lambda} + \frac{1}{4}\left( \frac{\pi}{2} + i \gamma \right)\frac{1}{\lambda} \end{align}

We can get the error estimate by noting that $\int_1^\infty dt\ e^{i \lambda t} \ln t$ and $\int_0^1 dt\ e^{i \lambda t}\left( \text{analytic in } t \right)$ are both $O\left(\frac{1}{\lambda}\right)$ as $\lambda\to\infty$.