an exercise in distribution theory
DOCX · 42.3 KB
Open DOCX file
A short worked note by Phil dated 5.24.10 in distribution theory. He tests the claim that the integral of cos(k'x)cos(kx) over x equals πδ(k-k') using a Gaussian test function and integrals from Gradshteyn-Ryzhik. He finds the symmetric result πδ(k-k')+πδ(k+k'), gets the sine analog with a minus sign, checks consistency with the Fourier delta, and confirms it via exponentials.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
An Exercise in Distribution Theory PhL 5.24.10
How would you verify (or disprove) this claimed fact?
!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k')
In distribution theory, we could come up with a test function φ(k) being some kind of smooth thing with finite support, and we would then say
!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k') =>
φ(k')!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k') φ(k')
∫dk' φ(k')!Syntax Error, Idx cos(k'x) cos(kx) = ∫dk' π δ(k-k') φ(k')
so this last line is the "meaning" of the distributional statement of the first line. At this point it would be good to ask about the integration range of k. Whatever it is, we will assume that we have chosen k' to be within the range. Then we have
∫dk' φ(k')!Syntax Error, Idx cos(k'x) cos(kx) = π φ(k)
We might then continue along and try to say
LHS = !Syntax Error, Idx cos(kx) ∫dk' φ(k') cos(k'x)
Now the range of the dk' integration seems perhaps more important.
We might follow page 33 Stak and choose the following test function
φa(k) = exp(a2/ (k2- a2) |k| < a and 0 outside
This function is maximal when the exponent is maximal which is when k2 = 0 and we have e-1 = 1/e. So this function peaks at k = 0 and has half width a. Let's now move this so the center is somewhere else
φa;k0(k) = exp( a2/ ([k-k0]2- a2) |k-k0| < a and 0 outside
Now assume for the moment that the range of dk' is (-∞,∞). Then we can consider:
∫dk' φ(k') cos(k'x) = !Syntax Error, Idk exp(a2/ ([k-k0]2- a2) cos(kx)
Is this a doable integral? I don't recall Stak every trying an integral like this. And I don't see anything in GR that involves something like exp[ 1/(ax2+ bx + c)] but I'm sure you could do fractions and relate to simpler integrals, but it sounds unpleasant. I would rather use a Gaussian test function even though it does not have finite support, hope to get away with this. So let's try this Gaussian form
φa;k0(k) = exp( - [k-k0]2/a2)
I am then trying to show this is true
!Syntax Error, Idx cos(kx) !Syntax Error, Idk' φa;k0(k') cos(k'x) = π φa;k0(k)
which is to say, I want to show this is true:
!Syntax Error, Idx cos(kx) !Syntax Error, Idk' exp( - [k'-k0]2/a2) cos(k'x) = π exp( - [k-k0]2/a2)
So let's go compute this integral
!Syntax Error, Idk' exp( - [k'-k0]2/a2) cos(k'x)
Shift integration variable to k" = k'-k0 to get
= !Syntax Error, Idk" exp( -k"2/a2) cos[k"x + k0x]
I find this on GR7 page 488
so set q = 1/a and x = k" and p = x and pλ = xλ = k0x so λ = k0 so we have
= a exp(-x2a2/4) cos(k0x)
So we now want to show that this is true:
!Syntax Error, Idx cos(kx) a exp(-x2a2/4) cos(k0x) = π exp( - [k-k0]2/a2)
Notice how the gaussian is "controlling" the integrals, making them doable. So now we face
a more complicated integral
!Syntax Error, Idx exp(-x2a2/4) cos(kx) cos(k0x) = 2 !Syntax Error, Idx exp(-x2a2/4) cos(kx) cos(k0x)
Luckily this also appears in GR7 page 489
so we have x = x, β = (a/2)2 = a2/4 a = k, b = k0 so we find // 1/(4β) = 1/a2
!Syntax Error, Idx exp(-x2a2/4) cos(kx) cos(k0x) = (1/4) (2/a) { exp(- [k-k0]2/a2) + exp(- [k+k0]2/a2) }
= (1/2a) { exp(- [k-k0]2/a2) + exp(- [k+k0]2/a2) }
So we now have
!Syntax Error, Idx cos(kx) a exp(-x2a2/4) cos(k0x) =
= a * 2 * (1/2a) { exp(- [k-k0]2/a2) + exp(- [k+k0]2/a2) }
= π { exp(- [k-k0]2/a2) + exp(- [k+k0]2/a2) }
= π { φ(k) + φ(-k) }
which is NOT the right answer, but it is sort of "close". We have an extra second term that we don't want to be there! Is this an artifact because I failed to use a proper test function? I don't think so.
Let's go back to our "claimed fact"
!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k')
Mystery question: the LHS is obviously symmetric under k → - k, but that is not true for the RHS! So something is amiss here! Stak has lots of delta function stuff involving sums of trig functions. On page 284 Stak gets into the singular BC world and on page 284 he considers the sine transform situation which he defines on x in (0,∞) only. Then on page 286 this leads to the π/2 delta function completeness. He shows in 4.58 how the issue of the "measure" works out and why you get just dk = dν for him. He claims that the δ completeness "holds rigorously in distribution theory".
Now, from the work I just did above, I think I have in fact proved this result:
!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k') + π δ(k+k')
where both terms must be present if your "range" for k is (-∞,∞). This then has the proper symmetries under negating one or both k factors. How can I show this? In the above work I had πδ(k-k') on the RHS. This caused the integrated RHS to be π φ(k). If instead I had πδ(k+k') = πδ(-k-k') on the RHS, the integrated RHS would have been π φ(-k). If I then added the two RHS's, I would have gotten the integrated RHS to be π [φ(k) + φ(-k) ]. And this is what I got! So I think I can say that I have used distribution theory to prove the above fact. And there would be a corresponding other fact
!Syntax Error, Idk cos(kx) cos(kx') = π δ(x-x') + π δ(x+x')
which you would verify in the same manner. Then of course if your problem is x in (0,∞), you could delete the second term.
I don't think I have ever seen these integrals with the two terms anywhere. So this throws a new monkey wrench into things.
Would the same thing be true for sines? My guess would be this
!Syntax Error, Idx sin(k'x) sin(kx) = π δ(k-k') - π δ(k+k')
It then has the right symmetries. If I were to redo the above analysis, I would have some different integrals. First, we had this one from page 488
and the one for sines is very similar, same page,
the only difference being the appearance of the letters "sin" in place of "cos" on the right. We would then need to look up an integral with two sines. We had this one from page 489
and the sine one would be this
and there you see the minus sign we want to see on the second term.
So here are my conclusions:
!Syntax Error, Idx cos(k'x) cos(kx) = π δ(k-k') + π δ(k+k')
!Syntax Error, Idx sin(k'x) sin(kx) = π δ(k-k') - π δ(k+k')
and you can then write these in various other ways, such as
!Syntax Error, Idx cos(k'x) cos(kx) = (π/2) δ(k-k') + (π/2) δ(k+k')
!Syntax Error, Idx sin(k'x) sin(kx) = (π/2) δ(k-k') - (π/2) δ(k+k')
and other forms obtained by doing x ↔ k and x'↔k'
Now let's go back to the classic Fourier idea
!Syntax Error, Idk e-ikx = 2π δ(x)
!Syntax Error, Idx e-ikx = 2π δ(k)
Can I relate these to my dual -delta expression above? The first could be written
!Syntax Error, Idk e-ik(x-x') = 2π δ(x-x')
or
!Syntax Error, Idk cos[k(x-x')] = 2π δ(x-x')
or
!Syntax Error, Idk [ cos(kx) cos(kx') + sin(kx) sin(kx')] = 2π δ(x-x')
Now using the rules given above,
!Syntax Error, Idk cos(kx) cos(kx') = π δ(x-x') + π δ(x+x')
!Syntax Error, Idk sin(kx) sin(kx') = π δ(x-x') - π δ(x+x')
we can add and the first term doubles and the second cancels, QED! So maybe I have resolved my mystery finally.
Here us another way to verify my claims:
!Syntax Error, Idx cos(k'x) cos(kx) = (1/4) !Syntax Error, Idx [eik'x + e-ik'x] [eikx + e-ikx]
= (1/4) !Syntax Error, Idx [ ei(k'+k)x + ei(k'-k)x + ei(-k'+k)x + ei(-k'-k)x ]
= (1/4) 2π [δ(k'+k) + δ(k'-k) + δ(-k'+k)+ δ(-k'-k)]
= π [δ(k'+k) + δ(k'-k) ]
and that was a lot easier than what I did above. Live and learn! Same idea for the sine result.