Berechnung des Erwartungswerts für kinetische Energie ⟨Ek⟩⟨Ek⟩\langle E_k \rangle für eine bekannte Wellenfunktion

Ich habe eine Wellenfunktion ( A = 1 N M ):

ψ = A X exp [ X 2 2 A ]

für die ich bereits den Normalisierungsfaktor berechnet habe (in meinem anderen Thema ):

A = 2 A π A = 1.06 1 N M N M

Was ich wissen möchte, ist, wie man den Erwartungswert für eine kinetische Energie berechnet. Ich habe versucht, es analytisch zu berechnen, aber ich verliere mich in der Integration:

E k = ψ ¯ T ^ ψ D X = A X exp [ X 2 2 A ] ( 2 2 M D 2 D X 2 A X exp [ X 2 2 A ] ) D X =

An diesem Punkt gehe ich und löse die zweite Ableitung und werde danach fortfahren:

= D 2 D X 2 A X exp [ X 2 2 A ] = A D 2 D X 2 X exp [ X 2 2 A ] = A D D X ( exp [ X 2 2 A ] 2 X 2 2 A exp [ X 2 2 A ] ) = = A ( 2 X 2 A exp [ X 2 2 A ] 1 A D D X X 2 exp [ X 2 2 A ] ) = = A ( X A exp [ X 2 2 A ] 2 X A exp [ X 2 2 A ] + X 3 A 2 exp [ X 2 2 A ] ) = = A ( 3 X A exp [ X 2 2 A ] + X 3 A 2 exp [ X 2 2 A ] )

Ok, jetzt kann ich mit der Integration fortfahren:

= A X exp [ X 2 2 A ] ( 2 2 M A ( 3 X A exp [ X 2 2 A ] + X 3 A 2 exp [ X 2 2 A ] ) ) D X = = A 2 2 2 M X exp [ X 2 2 A ] ( 3 X A exp [ X 2 2 A ] + X 3 A 2 exp [ X 2 2 A ] ) D X = A 2 2 2 M ( 3 X 2 A exp [ X 2 A ] X 4 A 2 exp [ X 2 A ] ) D X = A 2 2 2 M ( 3 X 2 A exp [ X 2 A ] X 4 A 2 exp [ X 2 A ] ) D X Wie löse ich das? =

Das ist der Punkt, an dem ich mir eingestand, dass ich in einem Integral verloren war und das WolframAlpha benutzte, um mir selbst zu helfen. Nun, ich habe ein seltsames Ergebnis . Mein Professor hat das irgendwie verstanden ( M ist eine Masse eines Elektrons), aber ich weiß nicht wie:

= 2 2 M 3 2 A = 3 2 4 M A = 0,058 e v

Kann mir jemand helfen, das letzte Integral zu verstehen? Wie kann ich es lösen? Ist es möglich, analytisch zu sein (es sieht so aus, als hätte es der Professor getan, aber ich bin mir nicht sicher)?

Wolfram/Mathematica kann diese Integrale definitiv ausführen, aber Sie haben bei der Eingabe einige Zeichen übersehen ( e X 2 / A anstatt e X 2 / A ). Wenn Sie Mathematica verwenden, ist es ein guter Rat, Mathematica das beizubringen A > 0 , dh verwenden Sie zum Beispiel "$Assumptions = {a>0}".
Nachdem ich die Eingabe in Wolfram Alpha korrigiert habe, erhalte ich ein Ergebnis 9 / 4 π A während A 2 = 2 / A π A und so bekomme ich immer noch ein anderes Ergebnis als mein Prof: A 2 2 2 M 9 / 4 π A = 2 2 A π A 2 M 9 / 4 π A = 9 2 4 A M

Antworten (2)

In Ihrem Problem benötigen Sie Integrale der Art:

ICH 2 N = X 2 N e X 2 A   D X

Beachte das erstmal ICH 0 = ( π ) 1 2 ( 1 A ) 1 2

Nun ist es leicht zu sehen, dass es eine Wiederholungsbeziehung zwischen den Integralen gibt:

ICH 2 N + 2 = ICH 2 N ( 1 A )

Zum Beispiel,

ICH 2 = ICH 0 ( 1 A ) = 1 2 ( π ) 1 2 ( 1 A ) 3 2 = 1 2 ( π ) 1 2   A 3 2

ICH 4 = ICH 2 ( 1 A ) = 3 2 1 2 ( π ) 1 2 ( 1 A ) 5 2 = 3 2 1 2 ( π ) 1 2   A 5 2

Eine allgemeine Formel lautet:

ICH 2 N = ICH 0   ( 2 N 1 ) ! !   ( A 2 ) N = ( π ) 1 2 2 N   ( 2 N 1 ) ! !   A N + 1 2
Wo ( 2 N 1 ) ! ! = ( 2 N 1 ) ( 2 N 3 ) . . . . . .5 .3 .1

Mein Professor für statistische Physik nennt diese Integrale Laplace-Integrale ICH ( H ) .

ICH ( H ) = 0 X H e A 2 X 2 D X

Beachten Sie, dass

X H e A 2 X 2 D X = 2 ICH ( H )

einige Werte

ICH ( 0 ) = π 2 A , ICH ( 1 ) = 1 2 A 2 , ICH ( 2 ) = π 4 A 3 , ICH ( 3 ) = 1 2 A 4 , ICH ( 4 ) = 3 π 8 A 5

Sie können Brute-Force durch Integrieren von Teilen loswerden X H und verwenden ICH ( 0 ) ein klassisches Ergebnis, oder Sie können Induktion verwenden H oder eine andere Methode.

Du denkst also, ich kann das Integral per partes lösen?
Ja, man muss es nur wissen ICH ( 0 ) und die Laufzeit verkürzen X H bis du bekommst ICH ( 0 ) , und natürlich Symmetrien wie das Integral einer ungeraden Funktion in einem symmetrischen Intervall ist Null und so weiter.