Wasserstoffatom in 2 Raumdimensionen mit logarithmischem Potential

Ich habe versucht, die Wellenfunktion eines Elektrons in einem Wasserstoffatom zu lösen, das auf nur zwei Dimensionen (im Flachland) beschränkt ist. Zuerst habe ich die elektrostatische potentielle Energie im Flachland herausgefunden. Das hat sich herausgestellt

U = Z e 2 2 π ϵ 0 Protokoll r .
Hier habe ich gewählt r = 1 als Oberfläche des Nullpotentials.

Die Schrödinger-Gleichung in zwei Dimensionen ist

( 1 r r ( r r Ψ ) + 1 r 2 θ 2 Ψ ) + U Ψ = E Ψ .
Um übermäßiges Tippen zu vermeiden, habe ich mich für die Einstellung entschieden / 2 μ zu einem. Verwenden Ψ = R ( r ) Y ( θ ) Die Gleichung kann in Radial- und Winkelgleichungen unterteilt werden.

Die Winkelgleichung ist

θ 2 Y = m 2 Y
und die Radialgleichung ist
1 r r ( r r R ) + ( E + Z e 2 2 π ϵ 0 Protokoll r m 2 r 2 ) R = 0.

Ich kann die Winkelgleichung leicht lösen und quantisieren, aber ich stecke seit Tagen in der Radialgleichung fest. Ich weiß nicht, wie ich die radiale Differentialgleichung lösen soll (ich bin nicht sehr gut darin, DE zweiter Ordnung zu lösen). Ich bin auch nicht in der Lage herauszufinden, wie man die Gesamtenergie des Elektrons quantisiert.

Fragen:

  1. Wie löst man die Radialwellenfunktionsgleichung und wie quantisiert sich die Energie in einem solchen Atom?

  2. Als ich nach Online-Quellen suchte, um Hilfe zu erhalten, stellte ich fest, dass jeder die potenzielle Energie in der Form betrachtete 1 / r anstatt Protokoll ( r ) . Nach dem Gaußschen Gesetz ist das elektrostatische Potential in zwei Dimensionen logarithmisch. Warum halten sie das Potenzial für die Form 1 / r ?

Wikipedias Numerische Methoden für ODEs erwähnt mehrere Algorithmen.
Denken Sie nie, dass Sie der Erste sind, der über ein Problem nachgedacht hat: newton.ex.ac.uk/research/qsystems/portnoi/JMP_43_4681.pdf und web.math.ku.dk/~schlicht/DL/2013/HydrogenAtom. pdf
@Gert Ich habe beide Papiere, die Sie hier gepostet haben, bereits gesehen. Beide unterscheiden sich von dem Problem, das ich löse (oder löste). Das in beiden Artikeln verwendete Potential ist das Coulomb-Potential, das im Flachland nicht anwendbar ist, es ist nur in begrenzten zwei Dimensionen im 3D-Raum anwendbar (wie von @Vadium erklärt). Ich habe in meinem Beitrag auch nie gesagt, dass ich der erste bin, der darüber nachdenkt, ich habe nur gesagt, dass ich hier ein Problem habe und Hilfe brauche.

Antworten (2)

Ich stimme zu, dass diese Gleichung wahrscheinlich keine Lösung in Form von Elementarfunktionen zulässt. Mit ein wenig Graben fand ich jedoch dieses Papier: Atabek et. Al. Phys. Rev. A 1974 zu Ihren Fragen (sie verwenden ein logarithmisches 2D-Potential und diskutieren das Spektrum der Energieeigenwerte).

Eine nützliche Sache, die sie tun, ist Ersatz R m = r ½ f m den Begriff proportional zu eliminieren R ' in die Gleichung und um eine Gleichung der Form zu erzeugen: f m ' ' = g ( r , m ) f m womit man besser arbeiten kann.

Zur numerischen Lösung:

Was Sie hier haben, ist ein "2-Punkt-Randwertproblem (BVP)" (da Sie den Wert von angeben R m , oder gleichwertig f m , bei r = 0 und Sie möchten, dass es verschwindet r ), die häufig mit sogenannten „Shooting Codes“ behandelt werden.

Wenn Sie mit einfacheren (einstufigen, expliziten) numerischen Methoden wie RK4 vertraut sind, können Sie einfach über Werte von iterieren E mit einer Vermutung für R ' ( 0 ) und überprüfen Sie den Wert der Lösung ( R ) bei einem großen Endwert von r . Abhängig von Ihren Einheiten sogar so etwas wie r f = 10 scheint zu genügen. Dann können Sie Ihr Suchraster manuell anpassen E bis zu diesem Endwert, R ( r f ) , ist ungefähr null.

Ich habe dieses RK4-Iterationsschema kodiert, nur um die qualitative Form der radialen Wellenfunktionen zu überprüfen (mit m = 0 ) in der verlinkten Veröffentlichung von 1974 berichtet und sie scheinen damit übereinzustimmen (vgl. Abb. 4 der Veröffentlichung).

Eine Einschränkung : Achten Sie nicht auf die numerischen Werte von E oder R m = 0 In der folgenden Abbildung sind die Einheiten wahrscheinlich etwas verwechselt. Außerdem normalisierte ich die Spitzen von R m zu 1 zum Vergleich quer m 's.

Trotzdem sind hier zur qualitativen Überprüfung drei Eigenfunktionen der Gleichung, eine für m = 0 , und zwei für m = 1 . Diese wurden mit der iterativen RK4-Methode gefunden.

Geben Sie hier die Bildbeschreibung ein

Ich bezweifle, dass diese Gleichung lösbar ist, obwohl es ratsam wäre, in einem Buch über spezielle Funktionen wie Abramovitz & Stegun oder Gradshtein & Ryzhik nachzuschlagen.

Ihre Gleichung scheint für ein Flachland richtig zu sein . Viele Probleme betreffen jedoch die reale 3D-Welt, in der die Bewegung auf zwei oder sogar eine Dimension beschränkt ist. In diesem Fall bleibt das Potenzial 3D: 1 / r . Der eindimensionale Fall ist bemerkenswert, da die Bindungsenergie divergiert, was lange Zeit ein Problem für die Analyse von Exzitonen in Kohlenstoffnanoröhren darstellte.