twistor59 wrote:A quick google around reveals some papers which purport to be about the 2D hydrogen atom, like
this one. However they treat the Coulomb potential as the 3d function (i.e. 1/r) just restricted to 2 dimensions instead of using the "correct" 2D potential value, which is log(r).
Turns out to have a solution similar to the full-3D solution.
On a related note, I had a go a while back at computing the Klein Gordon propagator in various dimensions.
In effect, we must solve D
x2P - d
2P/dt
2 - m
2P = delta(x)*delta(t)
Our first stop is
P(x,t) = integrate over all p
n,E : (1/(2pi)
n+1)) * (exp(i*p.x-i*E*t)/(E
2 - p
2 - m
2 + i*eps)
I found
a paper on doing it in 4D, and we can generalize from there. Equation 10:
P(x,t) = integrate over all p
n : -(i/2) * (1/(2pi)
n)) * (exp(i*p.x-i*E*t)/E)
where E = sqrt(p
2 + m
2)
Integrating over all time gives a pure-space version:
P(x) = integrate over all p
n : - (1/(2pi)
n)) * (exp(i*p.x)/E
2)
One can do these integrals with the help of
Abramowitz and Stegun: Handbook of Mathematical Functions -- yes, all the formulas are scanned into there, though not much of the tables. I also have the help of Mathematica -- computer-algebra software. If you are going to do a lot of theoretical work, I recommend getting some computer-algebra software.
Integrating over all angles we can do by turning the P(x,t) and P(x) integrals into
integrate p: 0 to +oo, a: 0 to pi : (1/(2pi)
n) * p
n-1 * (sin(a))
n-2 * W
n-1 * exp(i*p*r*cos(a)) *
{ - (i/2)*exp(-i*E/t)/E or - 1/E
2}
(A&S 9.1.20, p.360) where the n-D total solid angle W
n = 2*pi
n/2/Gamma(n/2) -- Gamma(n+1) = n!.
One gets for P(x,t) and P(x):
integrate p: 0 to +oo : (1/(2pi)
n/2) * J
n/2-1(p*r)/(p*r)
n/2-1 p
n-1 * { - (i/2)*exp(-i*E/t)/E or - 1/E
2}
where J is a Bessel function. The P(x) one is relatively easy, at least with "Hankel-Nicholson integrals" (A&S 14.4.44, p.488):
P(x) = - (1/(2pi)
n/2 (m/r)
n/2-1 * K
n/2-1(m*r)
For n >= 3, the m -> 0 limit is
- (1/pi
n/2) * Gamma(n/2-1) * r
-(n-2) = - 1/(W
n*(n-2)*r
n-2)
(A&S 9.6.9, p.375) -- which one can find more directly.
For n = 2, that limit is + (1/(2pi)) * (log(m*r/2) + (Euler-Mascheroni constant))
For n = 1, that limit is - (1/2)*(1/m - r)
We get some simple forms for some of the solutions:
For n = 1, - (1/2) * (1/m) * exp(-m*r)
For n = 3, - (1/(4*pi)) * (1/r) * exp(-m*r)
With the help of (A&S 9.1.31 p.361, 9.6.28 p.376), different dimensions are interrelated by
P
n+2(r) = - (1/(2pi)) * (1/r)(d/dr)P
n(r)
something also true of P(r,t). This means that we only need to solve the n = 1 and n = 2 cases of P(r,t); the rest is by
P
n+2(r,t) = - (1/(2pi)) * (1/r)(d/dr)P
n(r,t)
I'll save those n = 1 and n = 2 cases for my next post.