Volumenladungsdichte des H-Atoms

Ich habe ein Problem, bei dem ich die Volumenladungsdichte eines neutralen Wasserstoffatoms berechnen soll. Das Potenzial ist gegeben

Φ = k e A R R ( 1 + A R 2 )
Jetzt habe ich versucht, die Poisson-Gleichung zu verwenden
Δ Φ = ρ ε 0
was mich zu führt
ρ = Δ ( Q 4 π ϵ 0 = k e a R R ( 1 + a R 2 ) ) = k Δ ( e a R R + a e a R 2 ) = k ( Δ ( e a R R ) + a 2 Δ ( e a R ) )
Jetzt definiere ich F = e a R Und G = 1 R . Der Laplace-Operator des Produkts F G ist dann
Δ ( F G ) = G Δ ( F ) + F Δ ( G ) + ( F ) ( G )
und die Derivate sind
( F ) = a e a R R ^ Δ ( F ) = a 2 e a R
( G ) = 1 R 2 R ^ Δ ( G ) = 4 π δ ( R )
( F ) ( G ) = a e a R R 2
Das Wiedereinsetzen in die ursprüngliche Gleichung ergibt
ρ = k e a R ( a 2 R 4 π δ ( R ) + a R 2 + a 3 2 )
Dies scheint mir jedoch etwas falsch zu sein, da ich erwartet hätte, dass der Ausdruck vom Ursprung an zunimmt und dann nach einiger Zeit abnimmt R = R da das Potential der Elektronenhülle übernehmen sollte.

Kann mir jemand bestätigen, dass dies richtig ist, oder mir zeigen, wo ich den Fehler gemacht habe?

Abgesehen davon, dass ich die Ableitungen wie in kartesischen Koordinaten genommen habe, habe ich versucht, den Laplace-Operator durch Berechnung in Kugelkoordinaten sowie unter Verwendung des sphärischen Laplace-Operators zu berechnen

Δ Φ = 1 R 2 R ( R 2 Φ R )
aber immer noch das gleiche Ergebnis.

Es lohnt sich immer daran zu denken, dass die N = 0 Die Elektronenwellenfunktion ist im Zentrum maximal ...
Sie sagen also, dass die Berechnung irgendwie konsistent ist, da meine Elektronenladung die höchste Dichte bei hat R = 0 ? @dmkee

Antworten (1)

Mir ist nicht klar, warum Sie mit der Antwort, die Sie erhalten, unzufrieden sind. Ich würde vorschlagen, Ihre Erwartungen in Bezug auf die Gesamtladung zu formulieren, die in einer Radiuskugel enthalten ist R ,

4 π 0 R ρ ( R ' ) R ' 2 D R ' .
Dies sollte die positive Ladung des Kerns an geben R 0 + (weil das Proton in diesem Modell punktgroß ist!) und fallen dann monoton auf Null ab R steigt durch 1 / a und darüber hinaus zu R , und Ihre Sphäre enthält immer mehr von der Elektronenwolke, die den Kern umgibt (und neutralisiert). Wenn dies fehlschlägt, müssen Sie unbedingt Ihre Berechnungen überprüfen.

Beachten Sie auch, dass Verwirrung zwischen δ ( R ) Und δ ( R ) kann hier ziemlich schädlich sein.
Ok, danke an alle, ich habe die Berechnungen neu gemacht, indem ich zuerst r>0 betrachtete und danach r->0 und eine Taylor-Entwicklung des Potentials bei ließ R = 0 Jetzt habe ich
ρ = ε 0 Φ = Q ( δ ( R ) a 3 8 π e a R )
Was mir ein bisschen besser erscheint, weil es das nicht mehr enthält 1 R Und 1 R 2 Bedingungen
Nein, die Singular-Begriffe sollten da sein (oder zumindest nicht schaden). Das ist egal ρ haben eine Singularität, solange die Singularität integrierbar ist, was bedeutet, dass Singularitäten des Typs R S ganz bis zu mit S < 3 sind erlaubt.