Solving Difficult Integral: Exploiting Orthogonality of Plane Waves

  • Thread starter Morberticus
  • Start date
  • Tags
    Integral
In summary: The integral has shown up due to the Fourier transform:|r-r'|^{-1} = (1/2\pi^2)\int x^{-2} e^{ix(r-r')}dxMy intuition is telling me that, as (r-r') gets larger, the integral should get smaller, as |r-r'|^-1 will get smaller. Or, to use the variable I have been using up to now, |k|^-1 goes to 0 as k goes to infinity. Would this also imply that f(k) goes to infinity as k goes to 0? If
  • #1
Morberticus
85
0
I'm trying to sort out the integral:

[tex]\int x^{-2}e^{ikx}dx[/tex].

At first, I thought I would have to solve it numerically. But I am hoping to exploit the orthogonality of the plane waves:

[tex]\int e^{ikx}dx = 2\pi\delta(k)[/tex]

It has been a while since I have dealt with integral manipulation, so any help would be appreciated.
 
Physics news on Phys.org
  • #2
what are the limits?
 
  • #3
Dickfore said:
what are the limits?

The limits are at minus infinity and infinity.

[edit]-actually good point. If the integral is infinite then I will have to introduce limits.
 
  • #4
The integral is not convergent in the Riemann sense. It has a pole of order 2 at x = 0. Nevertheless, you can formally try this trick. Consider the integral as a function of the parameter k:

[tex]
f(k) \equiv \int_{-\infty}^{\infty}{x^{-2} \, e^{i \, k \, x} \, dx}
[/tex]

Take two derivatives with respect to k. What do you get?
 
  • #5
Dickfore said:
The integral is not convergent in the Riemann sense. It has a pole of order 2 at x = 0. Nevertheless, you can formally try this trick. Consider the integral as a function of the parameter k:

[tex]
f(k) \equiv \int_{-\infty}^{\infty}{x^{-2} \, e^{i \, k \, x} \, dx}
[/tex]

Take two derivatives with respect to k. What do you get?

Hmmm... My guess would be

[tex]f''(k) = \int^{\infty}_{-\infty} e^{ikx}dx[/tex]

which is infinite, right?

It's funny that you mention f(k) because that is precisely what I would need, as the integral will be nested in a larger, numerical integral over k, and I was hoping to save on a lot of computational time by getting f(k) analytically.
 
  • #6
You made a sign error and the integral you got is of the same form as one of the integrals you had already posted.
 
  • #7
Dickfore said:
You made a sign error and the integral you got is of the same form as one of the integrals you had already posted.

Ah yes, I forgot to bring down the 'i' twice (I'm horribly rusty).

So I would get:

[tex]f''(k) = -\int^{\infty}_{-\infty}e^{ikx}dx = -2\pi\delta(k)[/tex]

Though is there a way I can implement this? Would I be allowed to, say, do the following:

[tex]f(k) = \frac{2\pi k^2}{2}\delta(k)[/tex]

I.e. re-integrate f''(k)
 
  • #8
Yes, you are allowed to re-integrate, but your integration is incorrect.
 
  • #9
Dickfore said:
Yes, you are allowed to re-integrate, but your integration is incorrect.

Forgot completely that delta was a function of k.

So would these steps be correct?

[tex]f'(k) = 2\pi[/tex]

[tex]f(k) = 2\pi k[/tex]

[edit]-Oops, sign error. But is the form correct?
 
  • #10
Hint:

Consider the Heaviside step function:

[tex]
\theta (k) = \left\{\begin{array}{rcl}
1 &,& k \ge 0 \\

0 &,& k < 0
\end{array}\right.
[/tex]

What is the derivative of this function?
 
  • #11
Dickfore said:
Hint:

Consider the Heaviside step function:

[tex]
\theta (k) = \left\{\begin{array}{rcl}
1 &,& k \ge 0 \\

0 &,& k < 0
\end{array}\right.
[/tex]

What is the derivative of this function?

Hmmm... The derivative is a delta function. But would that mean I would have to manipulate something like:

[tex]\theta(\infty) = \int^{\infty}_{-\infty} \delta(k)dk [/tex]

as

[tex]f''(k) = -2\pi\delta(k) [/tex]

[tex]f'(k) = -2\pi \int^{\infty}_{-\infty} \delta(k)dk = -2\pi\theta(\infty)[/tex]

[edit]- Have to run but will be back tomorrow. Thanks a million. This should save me weeks of cpu time.
 
  • #12
Morberticus said:
Hmmm... The derivative is a delta function. But would that mean I would have to manipulate something like:

[tex]\theta(\infty) = \int^{\infty}_{-\infty} \delta(k)dk [/tex]

as

[tex]f''(k) = -2\pi\delta(k) [/tex]

[tex]f'(k) = -2\pi \int^{\infty}_{-\infty} \delta(k)dk = -2\pi\theta(\infty)[/tex]

[edit]- Have to run but will be back tomorrow. Thanks a million. This should save me weeks of cpu time.

When you integrate f''(k) you want a function of k again, which you don't get if you integrate over the entire k range. If you integrate over the entire range of k there's no k left.

[tex]f'(k) = f'(-\infty) - 2\pi \int_{-\infty}^k dk'~\delta(k') = -2\pi \theta(k)[/tex]
assuming f and f' vanish at the boundaries. You then do something similar for the next integration.
 
  • #13
Mute said:
When you integrate f''(k) you want a function of k again, which you don't get if you integrate over the entire k range. If you integrate over the entire range of k there's no k left.

[tex]f'(k) = f'(-\infty) - 2\pi \int_{-\infty}^k dk'~\delta(k') = -2\pi \theta(k)[/tex]
assuming f and f' vanish at the boundaries. You then do something similar for the next integration.

I have read up a little on the method you have used (dummy variable?) and I think I have the hang of it. But would it be valid to assume that f(k) and f'(k) are bounded? Would this need to be tested?

The integral has shown up due to the Fourier transform:

[tex]|r-r'|^{-1} = (1/2\pi^2)\int x^{-2} e^{ix(r-r')}dx[/tex]

My intuition is telling me that, as (r-r') gets larger, the integral should get smaller, as |r-r'|^-1 will get smaller. Or, to use the variable I have been using up to now, |k|^-1 goes to 0 as k goes to infinity. Would this also imply that f(k) goes to infinity as k goes to 0? If it does then I am an idiot, as I am integrating an expression that looks like:

[tex]\frac{1}{2\pi^2}\int^{a}_{b}sin(Cr)sin(Dr')[\int x^{-2} e^{ix(r-r')}dx]sin(Er)sinF(r')drdr'[/tex]and the purpose of the nested integral was ultimately to avoid an infinity when r = r'. I am just copping that the function f(k) might be infinite at k = 0? In that case, I am assuming adding finite limits to the nested integral would give me a finite approximation of f(k) at k = 0 which would suit my purposes fine.

Thanks to both you and Dickfore for your help so far, as I think I better understand how to manipulate the integral. My (hopefully) last question is, is the same line of logic you have shown me valid if the limits are finite?

[edit]- Actually, would I be right in thinking it would be easier, as all I would need is an expression as a function of k and two arbitrary bounds?
 
Last edited:
  • #14
I have tried integrating the intral in mathematica, and it has produced the expression:

[tex]\int x^{-2} e^{ikx}dx = ikEI(ikx) - \frac{e^{ikx}}{x}[/tex]

Ei(ikx) has a simple enough expression that I can approximate with a series expansion. Here is hoping that it will then simply be a case of inserting the two arbitrary bounds in place of x.
 
Last edited:
  • #15
Morberticus said:
I have tried integrating the intral in mathematica, and it has produced the expression:

[tex]\int x^{-2} e^{ikx}dx = ikEI(ikx) - \frac{e^{ikx}}{x}[/tex]

Ei(ikx) has a simple enough expression that I can approximate with a series expansion. Here is hoping that it will then simply be a case of inserting the two arbitrary bounds in place of x.

What limits did you use though?
 
  • #16
Mute said:
When you integrate f''(k) you want a function of k again, which you don't get if you integrate over the entire k range. If you integrate over the entire range of k there's no k left.

[tex]f'(k) = f'(-\infty) - 2\pi \int_{-\infty}^k dk'~\delta(k') = -2\pi \theta(k)[/tex]
assuming f and f' vanish at the boundaries. You then do something similar for the next integration.

For the last integration:

[tex]
f(k) = f(-\infty) + \int_{-\infty}^{k}{f'(t) \, dt}
[/tex]

[tex]
f(k) = -2 \, \pi \, \int_{-\infty}^{k}{\theta(t) \, dt}
[/tex]

I would use integration by parts.
 
  • #17
Dickfore said:
What limits did you use though?

I assumed the limits could be inserted the usual way, to get the expression:

[tex]f(k) = ik(Ei(ika)-Ei(ikb)) + (\frac{e^{ikb}}{b}-\frac{e^{ika}}{a})[/tex]

Where a and b are arbitrary limits.
 
  • #18
Hmmm... I was hoping that, as the interval |a-b| approaches infinity, the integral would approach |k|^-1

But when I plot f(k) from my previous post, the real component oscillates, and the imaginary component appears to be quadratic in form. Increasing |a-b| only increases the frequency of oscillation.

I will give the integration by parts a try.
 
  • #19
Dickfore said:
For the last integration:

[tex]
f(k) = f(-\infty) + \int_{-\infty}^{k}{f'(t) \, dt}
[/tex]

[tex]
f(k) = -2 \, \pi \, \int_{-\infty}^{k}{\theta(t) \, dt}
[/tex]

I would use integration by parts.

They way I've approached it:

[tex]\int^{\infty}_{-\infty}\theta(t)dt = \int^{0}_{-\infty}\theta(t)dt + \int^{k}_{0}\theta(t)dt[/tex]

And since [tex]\theta(k) = 1[/tex] for all k > 0:

[tex]\int^{0}_{-\infty}\theta(t)dt + \int^{k}_{0}\theta(t)dt = 0 + \int^{k}_{0}1dt = k[/tex]

There is a fudge: I have divided the integral (-inf,inf) into (-inf,0] (0,inf] so 0 technically isn't a member of the second interval.
 
  • #20
Why do you take the upper limit of the integral to be [itex]\infty[/itex]?
 
  • #21
Dickfore said:
Why do you take the upper limit of the integral to be [itex]\infty[/itex]?

Oops, typo. The limits are -infinity and k, and I divide the integral into -infinity,0 and 0,k.
 
  • #22
What happens if [itex]k < 0[/itex]?
 
  • #23
Dickfore said:
What happens if [itex]k < 0[/itex]?


Hmm... If k < 0 then I am guessing the integral is 0, as [tex]\theta(t) = 0[/tex] across the interval.

So would I end up with:

[tex]f(k) = -2\pi\int^{k}_{-\infty}\theta(t)dt = -2\pi k\theta(k)[/tex]
 
  • #24
doesnt the person only have to do intergration by parts:

UV - ∫vdu

to solve the problem?


or am i thinking of something else?
 
  • #25
Just to say

[tex]\int^{\infty}_{-\infty}e^{ikx}dx[/tex] diverge.
 
  • #26
Morberticus said:
Hmm... If k < 0 then I am guessing the integral is 0, as [tex]\theta(t) = 0[/tex] across the interval.

So would I end up with:

[tex]f(k) = -2\pi\int^{k}_{-\infty}\theta(t)dt = -2\pi k\theta(k)[/tex]

Hmm, if you take the inverse Fourier transform of this:

[tex]
\int_{-\infty}^{\infty}{(-2 \, \pi \, k \, \theta(k) \, e^{i k x} \, \frac{d k}{2 \pi}}
[/tex]

I don't think you will get [itex]x^{-2}[/itex].

The problem with the convergence properties of the integral is, of course, the pole of order 2 on the real axis at [itex]x = 0[/itex]. Let's try something different: Let's consider x as a complex variable and shift the pole either above or below the integration path by an infinitesimal amount and let us close the contour of integration either along a large semicircle in the upper or lower part of the imaginary plane. To see where to close it, consider the absolute value of the factor:

[tex]
|\exp(-i \, k \, z)| = e^{k \, \Im{(z)}}
[/tex]

and require that it tends to zero when [itex]z \rightarrow \infty[/itex]. We see that we must have [itex]k \, \Im{(z)} < 0[/itex], i.e. close it in the upper half-plane if [itex]k < 0[/itex] (notice that your contour is in the counterclockwise direction) or lower half-plane if [itex]k > 0[/itex] (the contour is in the clockwise direction). Also, the residue of the function:

[tex]
f(z) = \frac{\exp{(-i \, k \, z)}}{(z \mp i \, \epsilon)^{2}}
[/tex]

at the residue at the double pole is:

[tex]
\mathrm{Res}\left[f(z), z = \pm i \, \epsilon \right] = \lim_{z \rightarrow \pm i \, \epsilon}{\frac{d}{d z}\left(\exp{(-i \, k \, z)} \right)} = -i \, k \, e^{\pm k \, \epsilon}
[/tex]

Make sure you shift the pole so that it is inside the contour, i.e. take [itex]+i \, \epsilon[/itex] if [itex]k< 0[/itex] or [itex]-i \, \epsilon[/itex] if [itex]k>0[/itex]. Both of these cases can be written as [itex]-i \, \mathrm{sgn}(k) \, \epsilon[/itex].

Noticing that the integral has a zero contribution over an infinite semi-circle, we have:

[tex]
\int_{-\infty}^{\infty}{(x + i \, \mathrm{sgn}(k) \, \epsilon)^{-2} \, e^{-i \, k \, x} \, dx} = \oint_{C_{k}}{\frac{\exp(-i \, k \, z)}{(z + i \, \mathrm{sgn}(k) \, \epsilon)^{2}} \, dz} = -\mathrm{sgn}(k) \, 2 \pi \, i \, \mathrm{Res}\left[f(z), z = -\mathrm{sgn}(k) \, i \, \epsilon \right] = -2 \, \pi \, k \, \mathrm{sgn}(k) \, e^{-\mathrm{sgn}(k) \, k \, \epsilon} = -2 \, \pi \, |k| \, e^{-|k| \, \epsilon}
[/tex]

Finally, taking the limit [itex]\epsilon \rightarrow 0+[/itex], we get:

[tex]
\int_{-\infty}^{\infty}{x^{-2} \, e^{-i \, k \, x} \, dx} = -2 \pi \, |k|
[/tex]

When I perform the inverse Fourier transform, I get [itex]2/x^{2}[/itex], so, maybe you need to take half this value.
 
  • #27
Dickfore said:
Hmm, if you take the inverse Fourier transform of this:

[tex]
\int_{-\infty}^{\infty}{(-2 \, \pi \, k \, \theta(k) \, e^{i k x} \, \frac{d k}{2 \pi}}
[/tex]

I don't think you will get [itex]x^{-2}[/itex].

The problem with the convergence properties of the integral is, of course, the pole of order 2 on the real axis at [itex]x = 0[/itex]. Let's try something different: Let's consider x as a complex variable and shift the pole either above or below the integration path by an infinitesimal amount and let us close the contour of integration either along a large semicircle in the upper or lower part of the imaginary plane. To see where to close it, consider the absolute value of the factor:

[tex]
|\exp(-i \, k \, z)| = e^{k \, \Im{(z)}}
[/tex]

and require that it tends to zero when [itex]z \rightarrow \infty[/itex]. We see that we must have [itex]k \, \Im{(z)} < 0[/itex], i.e. close it in the upper half-plane if [itex]k < 0[/itex] (notice that your contour is in the counterclockwise direction) or lower half-plane if [itex]k > 0[/itex] (the contour is in the clockwise direction). Also, the residue of the function:

[tex]
f(z) = \frac{\exp{(-i \, k \, z)}}{(z \mp i \, \epsilon)^{2}}
[/tex]

at the residue at the double pole is:

[tex]
\mathrm{Res}\left[f(z), z = \pm i \, \epsilon \right] = \lim_{z \rightarrow \pm i \, \epsilon}{\frac{d}{d z}\left(\exp{(-i \, k \, z)} \right)} = -i \, k \, e^{\pm k \, \epsilon}
[/tex]

Make sure you shift the pole so that it is inside the contour, i.e. take [itex]+i \, \epsilon[/itex] if [itex]k< 0[/itex] or [itex]-i \, \epsilon[/itex] if [itex]k>0[/itex]. Both of these cases can be written as [itex]-i \, \mathrm{sgn}(k) \, \epsilon[/itex].

Noticing that the integral has a zero contribution over an infinite semi-circle, we have:

[tex]
\int_{-\infty}^{\infty}{(x + i \, \mathrm{sgn}(k) \, \epsilon)^{-2} \, e^{-i \, k \, x} \, dx} = \oint_{C_{k}}{\frac{\exp(-i \, k \, z)}{(z + i \, \mathrm{sgn}(k) \, \epsilon)^{2}} \, dz} = -\mathrm{sgn}(k) \, 2 \pi \, i \, \mathrm{Res}\left[f(z), z = -\mathrm{sgn}(k) \, i \, \epsilon \right] = -2 \, \pi \, k \, \mathrm{sgn}(k) \, e^{-\mathrm{sgn}(k) \, k \, \epsilon} = -2 \, \pi \, |k| \, e^{-|k| \, \epsilon}
[/tex]

Finally, taking the limit [itex]\epsilon \rightarrow 0+[/itex], we get:

[tex]
\int_{-\infty}^{\infty}{x^{-2} \, e^{-i \, k \, x} \, dx} = -2 \pi \, |k|
[/tex]

When I perform the inverse Fourier transform, I get [itex]2/x^{2}[/itex], so, maybe you need to take half this value.

Strangely enough, I have been messing around with mathematica and am getting something vaguely similar. Thanks.
 

1. What is the concept of orthogonality in solving difficult integrals?

Orthogonality refers to the mathematical property of two functions being perpendicular to each other. In the context of solving difficult integrals, it means that when two plane waves with different frequencies are multiplied together, the integral of their product over a specific range will be equal to zero. This property can be exploited to simplify complex integrals and make them easier to solve.

2. How does the orthogonality of plane waves help in solving difficult integrals?

The orthogonality of plane waves allows us to use the concept of Fourier series, which states that any periodic function can be represented as a sum of sine and cosine functions. By representing a difficult integral as a sum of simpler integrals involving plane waves, we can use the orthogonality property to eliminate certain terms and reduce the complexity of the integral.

3. Are there any limitations to using the orthogonality of plane waves in solving integrals?

While the orthogonality of plane waves can be a powerful tool in solving difficult integrals, it is limited to integrals involving periodic functions. It may not be applicable to integrals involving non-periodic functions or those with discontinuities.

4. Can the orthogonality of plane waves be used in higher dimensions?

Yes, the concept of orthogonality can be extended to higher dimensions, such as in solving integrals over surfaces or volumes. In these cases, plane waves are replaced by spherical or cylindrical waves, and the same principles of orthogonality can be applied to simplify the integrals.

5. How can I apply the concept of orthogonality to solve a specific difficult integral?

To apply the concept of orthogonality to a specific integral, you will need to first decompose the integral into a sum of plane waves using Fourier series. Then, you can use the orthogonality property to simplify the integral by eliminating certain terms. The final step is to solve the simpler integrals and combine the solutions to get the final answer.

Similar threads

Replies
2
Views
830
Replies
8
Views
1K
Replies
3
Views
1K
Replies
3
Views
1K
Replies
10
Views
2K
Replies
2
Views
1K
Replies
2
Views
2K
Replies
8
Views
3K
Replies
1
Views
830
Back
Top