Fourier transform/propagator

  • Thread starter RedX
  • Start date
  • Tags
    Fourier
In summary, the conversation discusses the uniqueness of Fourier transforms and how they can lead to different results in different gauges. The expert explains how the pole structure of the integrals affects the residues and how this can result in different answers. They also provide a mathematical explanation for why the residue in the second integral mentioned in the conversation is zero.
  • #1
RedX
970
3
I've got a question about Fourier transforms.

Say you have the integral:

[tex]D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2+i\epsilon} [/tex]

Can I just set this equal to:


[tex]D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+i\epsilon} [/tex]

since the pole structure of the denominator will force k to be on-shell?

But I thought Fourier transforms were unique? So how can two different transform functions give the same position space function?

You don't really see this type of question when using Lorentz covariant gauges, since here k_0 and vector k are split up in non-covariant way. But if you have a non-covariant gauge, you might see this.

It's just weird to me because I thought Fourier transforms were unique. Does this mean that I can use either momentum space representation of the propagator in Feynman diagrams, and I'll get the same answer?
 
Physics news on Phys.org
  • #2
This can't be right, but I can't figure out why.

Take: [tex]\frac{1}{k^2+i\epsilon}=P\frac{1}{k^2}-i\pi \delta(k^2) [/tex].

Then if:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2+i\epsilon }= \int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+i\epsilon }

[/tex]
this implies:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2}= \int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2}

[/tex]
[tex] \int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(\frac{\vec{k}^2}{k_{0}^2})}{k^2}= \int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2}[/tex]
where the integrals are over the principal part (so they are over the real line but avoid k2=0 ). But clearly this isn't true, or if it were, it would be amazing.

My math is terrible, so I've got to be doing something wrong, but to me the pole forces the photon to have its energy squared equal to the momentum squared, so I can't get around this.
 
  • #3
[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2+i\epsilon}
[/tex]

[tex]
D(x, y) = 0.5 \, \int{\frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+ i\epsilon}} + \frac{0.5}{k^{2}_{0}} \int{\frac{d^{4} k}{(2\pi)^{4}} \, e^{i k (x - y)}}
[/tex]

The second integral is:

[tex]
\delta^{4}(x - y)
[/tex]

and the first is just the propagator you thought your original expression would reduce to.
 
  • #4
Dickfore said:
[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2+i\epsilon}
[/tex]

[tex]
D(x, y) = 0.5 \, \int{\frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+ i\epsilon}} + \frac{0.5}{k^{2}_{0}} \int{\frac{d^{4} k}{(2\pi)^{4}} \, e^{i k (x - y)}}
[/tex]

The second integral is:

[tex]
\delta^{4}(x - y)
[/tex]

and the first is just the propagator you thought your original expression would reduce to.

If k doesn't have an arrow on it, I took it to mean a 4-vector. With an arrow on it, I meant 3-vector. So on the second integral you have, the numerator doesn't cancel with one of the factors in the denominator, and the integral over dk is the entire 4-vector space so you can't bring out the timed component ko2. I should have been clearer about notation.

I should also note that by having a ko in the denominator, that changes the pole structure, but superficially, since there is no residue at ko=0 . Indeed, if you want, you can swap [tex]\vec{k}^2 [/tex] with [tex]k_{o}^2 [/tex], so that the pole structure isn't changed when integrating over [tex]k_{o} [/tex] first.
 
  • #5
oh, yes, my bad. what is the signature of your metric?
 
  • #6
Dickfore said:
oh, yes, my bad. what is the signature of your metric?

I'm using +---.

In general, if your 4-momentum squared is the same sign as the epsilon, then that is +---. If they are opposite, then it is -+++.
 
  • #7
So:

[tex]
0.5 + 0.5 \frac{\vec{k}^{2}}{k^{2}_{0}} = 0.5 \left( 1 + \frac{k^{2} - k^{2}_{0}}{k^{2}_{0}}\right) = 0.5 \frac{k^{2}}{k^{2}_{0}}
[/tex]

The [itex]k^{2}[/itex] in the numerator cancels the one in the denominator and you get [itex]k^{2}_{0}[/itex] in the denominator. However, your initial integral had the pole w.r.t. [itex]k_{0}[/itex] at:

[tex]
k^{2}_{0} = \vec{k}^{2} - i \epsilon
[/tex]

so you need to correct:

[tex]
k^{2}_{0} \rightarrow k^{2}_{0} + i \epsilon
[/tex]

and your integral becomes:

[tex]
D(x, y) = 0.5 \, \int{\frac{d^{4}k}{(2 \pi)^{4}} \, \frac{e^{i k (x - y)}}{k^{2}_{0} + i \epsilon}}
[/tex]

which has a different pole structure than what you had assumed, but can be integrated by the same methods.
 
  • #8
thanks. you used -+++, but the gist of what you did is correct. changing to +---:

[tex]
0.5 + 0.5 \frac{\vec{k}^{2}}{k^{2}_{0}} = 0.5 \left( 1 + \frac{-k^{2} + k^{2}_{0}}{k^{2}_{0}}\right) = 1-0.5 \frac{k^{2}}{k^{2}_{0}}
[/tex]

[tex]
D(x, y) =\int{\frac{d^{4}k}{(2 \pi)^{4}} \, \frac{e^{i k (x - y)}}{k^2+i\epsilon} - 0.5 \, \int{\frac{d^{4}k}{(2 \pi)^{4}} \, \frac{e^{i k (x - y)}}{k^{2}_{0} + i \epsilon}}
[/tex]

and the second integral is zero since the residue is zero.

So it is true that:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(.5+.5\frac{\vec{k}^2}{k_{0}^2})}{k^2+i\epsilon } = \int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+i\epsilon}
[/tex]

which is weird as I thought Fourier transforms were unique.
 
  • #9
RedX said:
[tex]
D(x, y) =\int{\frac{d^{4}k}{(2 \pi)^{4}} \, \frac{e^{i k (x - y)}}{k^2+i\epsilon} - 0.5 \, \int{\frac{d^{4}k}{(2 \pi)^{4}} \, \frac{e^{i k (x - y)}}{k^{2}_{0} + i \epsilon}}
[/tex]

and the second integral is zero since the residue is zero.

How is the residue zero?
 
  • #10
Dickfore said:
How is the residue zero?

Because of the k02 term in the denominator, that is a pole of order 2 when integrating over k0.

So what you do is multiply by the integrand by k02, and take one derivative, and set k0=0 to get your residue. In this case this will result in zero.

Another way to see it is to note that the Laurent expansion of 1/x^2 is just 1/x^2. In a general Laurent expansion:

[tex]f(z)=m/z^n+...n/z^2+a/z+b+cz+...dz^n [/tex]

the residue would just be 'a', or the coefficient of 1/z.

edit: oops, forgot about the exponential term. so maybe it's not zero.
 
  • #11
The basic idea I'm trying to get at is this:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}g(k)}{k^2-m^2+i\epsilon }
[/tex]

The integral over k0 is done by taking the residues at [tex]k_0=\pm \sqrt{\vec{k}^2+m^2} [/tex] to get: [tex]
D(x,y)=-i\theta(x_o-y_o)\int \frac{d^3k}{2 E_k(2\pi)^3} e^{-ik(x-y)}g(k)
-i\theta(y_o-x_o)\int \frac{d^3k}{2 E_k(2\pi)^3} e^{-ik(y-x)}g(k)

[/tex]

where now [tex]k_0= \sqrt{\vec{k}^2+m^2} [/tex], including the k in g(k).

Therefore from the outset you should be able to set the k in g(k) on-shell in the 4-dimensional integral representation:[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}g(k)}{k^2-m^2+i\epsilon }
[/tex]

But the math ain't really working out (I think it's complicated by the fact that I set m=0).
 
Last edited:
  • #12
Your g(k) adds further poles to the integrand.
 
  • #13
Dickfore said:
Your g(k) adds further poles to the integrand.

Well, say [tex]g(k)=\frac{k_{0}^2}{\vec{k}^2} [/tex], or [tex]g(k)=1+(k_{0}^2-{\vec{k}^2) [/tex]. These shouldn't affect the pole when integrating just over [tex]k_0[/tex]. So it seems that you should be able to set:

[tex]g(k)=\frac{k_{0}^2}{\vec{k}^2}=1[/tex]

and

[tex]g(k)=1+(k_{0}^2-{\vec{k}^2)=1 [/tex]

from the very beginning.

But if this is possible, then Fourier transforms are not unique, since g(k)=1 gives the same result when integrated over momentum space as the above g(k)'s.

Quantum mechanically Fourier transforms are unique because the momentum basis forms a linearly independent basis.
 
  • #14
So, how exactly is:

RedX said:
[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}g(k)}{k^2-m^2+i\epsilon }
[/tex]

the same as:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2-m^2+i\epsilon }
[/tex]

I don't know what math is not working out for you exactly.
 
  • #15
Dickfore said:
So, how exactly is:

[tex]

D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}g(k)}{k^2-m^2+i\epsilon }

[/tex]

the same as:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2-m^2+i\epsilon }
[/tex]

I don't know what math is not working out for you exactly.

If I choose g(k)=1, they are the same. But what if I choose [tex]g(k)=1+(k_{0}^2-\vec{k}^2) [/tex]. Will the two different g(k)'s give the same result, since eventually [tex](k_{0}^2=\vec{k}^2) [/tex] when integrating over ko?

I chose [tex]g(k)=.5(1+\frac{\vec{k}^2}{k_{0}^2}) [/tex], which when you set [tex](k_{0}^2=\vec{k}^2) [/tex], gives you 1 again. However, as you pointed out, the pole structure changes so this is a bad example, so then I switched it to:

[tex]g(k)=.5(1+\frac{k_{0}^2}{\vec{k}^2}) [/tex]

So basically choose any g(k) such that when you put k on-shell, it forces g(k)=1.
 
  • #16
RedX said:
If I choose g(k)=1, they are the same. But what if I choose [tex]g(k)=1+(k_{0}^2-\vec{k}^2) [/tex]. Will the two different g(k)'s give the same result, since eventually [tex](k_{0}^2=\vec{k}^2) [/tex] when integrating over ko?

I chose [tex]g(k)=.5(1+\frac{\vec{k}^2}{(k_{0}^2}) [/tex], which when you set [tex](k_{0}^2=\vec{k}^2) [/tex], gives you 1 again. However, as you pointed out, the pole structure changes so this is a bad example, so then I switched it to:

[tex]g(k)=.5(1+\frac{k_{0}^2}{\vec{k}^2}) [/tex]

So basically choose any g(k) such that when you put k on-shell, it forces g(k)=1.

Clearly you have some misunderstandings on how to apply the residue theorem in complex integration problems. The g(k) you chose, namely:

[tex]
g(k)=.5(1+\frac{\vec{k}^2}{(k_{0}^2})
[/tex]

has extra poles in the complex [itex]k_{0}[/itex]-plane, so your logic from post #11 does not apply here.
 
  • #17
Dickfore said:
Clearly you have some misunderstandings on how to apply the residue theorem in complex integration problems. The g(k) you chose, namely:

[tex]
g(k)=.5(1+\frac{\vec{k}^2}{(k_{0}^2})
[/tex]

has extra poles in the complex [itex]k_{0}[/itex]-plane, so your logic from post #11 does not apply here.

Right, I messed up on that one. But what if I had this one instead:

[tex]
g(k)=.5(1+\frac{k_{0}^2}{\vec{k}^2})
[/tex]

This does not introduce a pole in the complex [itex]k_{0}[/itex]-plane.

The idea is that two transform functions, one with g(k)=1 and one with [tex]
g(k)=.5(1+\frac{k_{0}^2}{\vec{k}^2})
[/tex], give the same result. This shouldn't be possible if the momentum basis is linearly independent.
 
  • #18
Do you have a non-zero mass or no?
 
  • #19
I just thought of something. All these "g(k)" you had presented are not Lorentz invariant, so I don't think your expressions make physical sense, because they have different value in different reference frames.
 
  • #20
Dickfore said:
Do you have a non-zero mass or no?

Since there are computers that can do simple Fourier transforms, I thought I'd experiment with some simple ones ( http://www.wolframalpha.com/input/?i=fourierx[(1)/(x^2+1)] and http://www.wolframalpha.com/input/?i=fourier[(-x^2)/(x^2+1)] ):

[tex]\frac{1}{\sqrt{2\pi}} \int \frac{e^{-ikx}}{x^2+1} dx=\sqrt{\frac{\pi}{2}}e^{|k|} [/tex][tex]\frac{1}{\sqrt{2\pi}} \int \frac{e^{-ikx}(-x^2)}{x^2+1} dx=\sqrt{\frac{\pi}{2}}e^{|k|}-\sqrt{2 \pi} \delta(k) [/tex]

Since the pole is at [tex]x=\pm i [/tex], one would think that the second transform should be the same as the first.

And they almost are the same except for a funny delta function. First of all the second integral doesn't even satisfy Jordan's lemma as the integrand doesn't go to zero as x goes to infinity. But I guess that's okay.

The only difference between the two Fourier transforms are at a single point, k=0. This would correspond to (x0-y0)=0 in the previous examples (I switched up k and x in this post).

Now when (x0-y0)=0, I think the propagator, defined as the time-ordered product, kind of screws up, because what it the fields are taken at the same time? Actually you can evaluate the propagator at equal times.

As for Lorentz invariance, it can be broken so long as it's part of the gauge condition. If the theory is gauge invariant, you can transform to a non-covariant gauge, thereby breaking Lorentz invariance, but your theory should still give you the same answer since your theory was gauge invariant.

addendum:

I guess you're right:

[tex]
D(x,y)=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(\frac{k_{0}^2}{\vec{k}^2})}{k^2+i\epsilon } =
\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}(\frac{k^2+\vec{k}^2}{\vec{k}^2})}{k^2+i\epsilon}=\int \frac{d^4k}{2 \pi^4} \frac{e^{ik(x-y)}}{k^2+i\epsilon}+
\delta(x_0-y_0)\int \frac{d^3k}{2 \pi^3} \frac{e^{i\vec{k}(\vec{x}-\vec{y})}}{\vec{k}^2+i\epsilon}

[/tex]

So the Fourier transforms aren't the same when the times are equal. The last term in fact is a Yukawa term and proportional to:

[tex]\delta(x_0-y_0) \frac{1}{|\vec{x}-\vec{y}|} [/tex]
 
Last edited:

1. What is a Fourier transform/propagator?

A Fourier transform/propagator is a mathematical operation that transforms a function from its original domain (such as time or space) to a new domain (such as frequency or momentum). It is commonly used in signal processing and quantum mechanics to analyze and manipulate data.

2. How does a Fourier transform/propagator work?

A Fourier transform/propagator involves breaking down a function into a series of sine and cosine waves, which are then combined to create a representation of the original function in the new domain. This allows for the analysis and manipulation of the function in a different way.

3. What is the difference between a Fourier transform and a Fourier propagator?

A Fourier transform is a mathematical operation that transforms a function from one domain to another, while a Fourier propagator is a function that propagates a wavefunction in time. In other words, a Fourier transform is a tool used for analysis, while a Fourier propagator is a tool used for prediction and simulation.

4. What are the applications of Fourier transforms/propagators?

Fourier transforms/propagators are used in a wide range of fields, including signal processing, image and sound analysis, communications, quantum mechanics, and more. They are particularly useful for analyzing periodic or oscillating functions, and for understanding how systems evolve over time.

5. Are there any limitations or drawbacks to using Fourier transforms/propagators?

One limitation of Fourier transforms/propagators is that they are best suited for functions that are smooth and well-behaved. They may not work well for functions with sharp discontinuities or singularities. Additionally, Fourier transforms/propagators can be computationally intensive for large datasets, and may require advanced mathematical knowledge to use effectively.

Similar threads

  • Quantum Physics
Replies
4
Views
792
Replies
6
Views
1K
  • Quantum Physics
Replies
2
Views
642
  • Quantum Physics
Replies
1
Views
1K
  • Quantum Physics
Replies
2
Views
768
  • Quantum Physics
Replies
13
Views
1K
  • Quantum Physics
Replies
13
Views
3K
  • Quantum Physics
Replies
0
Views
657
Replies
2
Views
243
  • Calculus and Beyond Homework Help
Replies
5
Views
320
Back
Top