Inverse Fourier Transform in 2-d

Click For Summary

Discussion Overview

The discussion revolves around the evaluation of a two-dimensional inverse Fourier transform integral involving a specific kernel. Participants explore various methods to solve the integral, including changing to polar coordinates, using identities related to gamma functions, and applying differentiation techniques.

Discussion Character

  • Exploratory
  • Technical explanation
  • Mathematical reasoning

Main Points Raised

  • Nick presents the integral and seeks advice on finding an analytic solution, expressing uncertainty about the tractability of the integral.
  • A participant corrects an earlier approach by suggesting a transformation to polar coordinates and provides a detailed breakdown of the angular integral.
  • Another participant proposes using a gamma function identity to transform the integral into a triple integral, suggesting that it may simplify the evaluation.
  • Nick shares a rewritten form of the integral, applying differentiation with respect to parameters and introducing Bessel functions, while acknowledging the possibility of algebraic errors in their derivation.

Areas of Agreement / Disagreement

Participants present multiple approaches to the problem, and while they share insights and corrections, there is no consensus on a definitive solution or method. The discussion remains open-ended with various competing views on how to tackle the integral.

Contextual Notes

Some methods proposed depend on specific transformations and identities that may not be universally applicable. The discussion includes unresolved mathematical steps and assumptions that could affect the validity of the approaches.

nickthequick
Messages
39
Reaction score
0
Hi all,

I've been trying to solve the following

I = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \frac{x}{(x^2+y^2+d^2)^{\frac{5}{2}}} e^{-i(kx+\ell y)} \ dx \ dy

where d,k,\ell are constants. I haven't been able to put this into a tractable analytic form and I figured I'd consult all of you experts for advice before I resorted to approximation methods. So does anyone see any obvious way of solving this?

Thanks!

Nick
 
Physics news on Phys.org
EDIT:
The angular integral is wrong!Go to polar coordinates. The integral over the angular coordinate is:
<br /> \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos (\phi - \phi_0)} \, d\phi}<br />
where q = \sqrt{k^2 + l^2}, and \tan \phi_0 = l/k. Due to the periodicity of the two functions, we can perform the translation \phi \rightarrow \phi + \phi_0, but keep the same range of integration. Then, we get:
<br /> \int_{0}^{2 \pi}{\cos (\phi + \phi_0) \, e^{-i q \, \rho \, \cos \phi} \, d\phi}<br />
Applying the addition theorem for the cosine, we get:
<br /> \cos \phi_0 \, \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} - \sin \phi_0 \, \int_{0}^{2 \pi}{\sin \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi}<br />
The two integrals are:
<br /> \int_{0}^{2 \pi}{\cos \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} = -2 \pi \, i \, J_1(q \, \rho)<br />
and
<br /> \int_{0}^{2 \pi}{\sin \phi \, e^{-i q \, \rho \, \cos \phi} \, d\phi} = 0<br />
(the last one may be seen by shifting \phi \rightarrow \phi + \pi. Then, we have an integral of an odd function over a symmetric interval.)

Then, the integral over the radial coordinate is:
<br /> I = -2 \pi \, i \, \cos \phi_0 \, \int_{0}^{\infty}{\frac{\rho^2}{(\rho^2 + d^2)^{5/2}} \, J_1(q \, \rho) \, d\rho}<br />
 
Last edited:
Perhaps you can make use of the identity

\frac{1}{z^\nu} = \frac{1}{\Gamma(\nu)} \int_0^\infty dt~t^{\nu-1} e^{-zt},
which holds for z, \nu &gt; 0.

You would then have a triple integral

I_{k\ell} = \frac{1}{\Gamma(5/2)}\int_{-\infty}^\infty dx \int_{-\infty}^\infty dy \int_0^\infty dt~x t^{3/2} \exp\left(-\left[ tx^2 + ikx + ty^2 + i\ell y + d^2t\right]\right).

This is Gaussian in each of x and y, so it looks like you should be able to at least reduce the integral to just one integral over t. Whether or not that will be easier than the integral Dickfore derived or doable at all, I don't know.
 
Hm. Thanks for the suggestions. Here's what I've come up with.

Rewrite I as (I've changed d to h here for clarity of presentation)

I = \frac{i}{3}\partial_k \partial_h\iint \frac{1}{(x^2+y^2+h^2)^{\frac{3}{2}}}e^{-i(kx+\ell y)} \ dx \ dy

=\frac{i}{3}\partial_k \partial_h\int_0^{2\pi}\int_0^{\infty} \frac{1}{(r^2+h^2)^{\frac{3}{2}}} e^{-i\vec{r}\cdot \vec{k}}r\ dr \ d\theta

WLOG, orient \vec{k}=(k,\ell) with the x-axis so that \vec{r}\cdot \vec{k} = |\vec{r}||\vec{k}| \cos \theta. Therefore, we have

I=-i\frac{2\pi}{3} \partial_k\partial_h\int_0^{\infty} \frac{J_o(|\vec{k}| r)}{(r^2+h^2)^{\frac{3}{2}}} r \ dr

where J_o is the Bessel function of the first kind. The integral on the RHS of the above equation can be solved in closed form, so that we find

I = -i\frac{2\pi}{3}\partial_k\partial_h K_o(h|\vec{k}|)= -i\frac{2\pi}{3} h k\ K_o(h\sqrt{k^2+\ell^2})

where K_o is a modified bessel function of the second kind. This might be as good as I can hope for. Also, I might have made some algebra mistakes along the way.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
12K
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
9K