How to analyze the asymptotic behaviour of this integral function? Based on the asymptotic analysis of correlation functions at large distence in Physics, now I get a math question. Let the function $$f(x)=\int_{-1}^{1}\sqrt{1-k^2}e^{ikx}dk.$$
Without working out the explicit form of this function, how do we know the asymptotic behaviour of $f(x)$ at large distence $\left |x  \right |\rightarrow+\infty$?

I going to manipulate the integral into a form that I can analyze the method of stationary phase.  Let $k=\cos{\theta}$ and the integral becomes
$$\begin{align}I(x) &= \int_0^{\pi} d\theta \, \sin^2{\theta} \, e^{i x \cos{\theta}}\\ &= \int_0^{\pi} d\theta  \, e^{i x \cos{\theta}} - \int_0^{\pi} d\theta \, \cos^2{\theta} \, e^{i x \cos{\theta}}\\ &= \int_0^{\pi} d\theta  \, e^{i x \cos{\theta}} + \frac{d^2}{dx^2} \int_0^{\pi} d\theta  \, e^{i x \cos{\theta}}\end{align}$$
Now, note that
$$\int_0^{\pi} d\theta  \, e^{i x \cos{\theta}} = \int_0^{\pi/2} d\theta  \, e^{i x \cos{\theta}} + \int_{\pi/2}^{\pi} d\theta  \, e^{i x \cos{\theta}} = 2 \Re{\left [ \int_0^{\pi/2} d\theta  \, e^{i x \cos{\theta}}\right]}$$
Now, we may apply stationary phase.  The stationary point of the integrand is at $\theta=0$; there, we may approximate the argument of the exponential by its Taylor expansion.  Further, because of the oscillatory cancllations, we may simply draw the upper limit of the integral to infinity to first order:
$$\int_0^{\pi/2} d\theta  \, e^{i x \cos{\theta}}  \sim e^{i x} \int_0^{\infty} d\theta \, e^{-(i x/2) \theta^2} = \frac12 e^{i(x-\pi/4)} \sqrt{\frac{2 \pi}{x}} \quad (x\to\infty)$$
Therefore
$$\int_0^{\pi} d\theta  \, e^{i x \cos{\theta}} \sim   \sqrt{\frac{2 \pi}{x}} \cos{\left ( x-\frac{\pi}{4}\right)}\quad (x\to\infty)$$
To get the asymptotic behavior of $I(x)$, it looks like we need to take the second derivative of the above result.  It turns out that
$$\frac{d^2}{dx^2} \left[x^{-1/2} \cos{\left ( x-\frac{\pi}{4}\right)}\right ]= \left (\frac{3}{4} x^{-5/2} - x^{-1/2}\right ) \cos{\left ( x-\frac{\pi}{4}\right)}  + x^{-3/2}\sin{\left ( x-\frac{\pi}{4}\right)}   $$
Note that the $x^{-1/2}$ piece drops out when combined with the original integral, and the $x^{-5/2}$ piece is subdominant.  Thus, the leading asymptotic behavior of the integral $I(x)$ is finally
$$\int_{-1}^1 dk \, \sqrt{1-k^2} \, e^{i k x}  \sim  \sqrt{2 \pi} x^{-3/2} \sin{\left ( x-\frac{\pi}{4}\right)}\quad (x\to\infty)$$
Here is a plot illustrating this behavior against the exact value of the integral (the lower plot is the asymptotic result derived here):

ADDENDUM
One objection to the above derivation might be that I neglected the next term in the asymptotic expansion of the first integral of $I(x)$, which we know goes as $A x^{-3/2} \sin{(x-\pi/4)}$.  How could I ignore that?  It turns out $I$ is the sum of this term and its second derivative, and the only term that remains $O(x^{-3/2})$ vanishes because it is a simple sine term, the sum of itself and its second derivative vanishing.  So the derivation of the asymptotic approximation above remains valid.
