Erhalten des Transmissionskoeffizienten des Strahls bei einem linearen Potential

Geben Sie hier die Bildbeschreibung ein

Ich möchte den Transmissionsgrad bestimmen T für einen Teilchenstrahl

Ψ ( X , T ) = A Ö e ich k X e ich E Ö T
mit Energie
E = 2 k 2 2 M
Einfall auf ein negatives lineares Potential (siehe Diagramm; ich entschuldige mich für seine Grobheit, ich habe es selbst gemacht) definiert durch:

v ( X ) = { 0 für  X < 0 | v Ö | ( X / L ) für  0 < X < L | v Ö | für  X > L

Der Transmissionskoeffizient T wird unter Verwendung des Wahrscheinlichkeitsdichtestroms gefunden J als:

T = | J T R A N S | | J ich N C |

Wo | J T R A N S | Und | J ich N C | sind die übertragenen und einfallenden Wahrscheinlichkeitsströme. Wahrscheinlichkeit aktuell J kann bezogen werden durch:

J = ich 2 M { Ψ ( X , T ) X ( Ψ ( X , T ) ) Ψ ( X , T ) X ( Ψ ( X , T ) ) }

Für den Fall eines einzelnen Step-down-Potentials (grüne Linie im Diagramm) können wir also die TISE in jeder Region lösen und Randbedingungen verwenden, um Koeffizienten und Abhängigkeiten aufzulösen k 1 Und k 2 , die entweder einfallende/reflektierte Strahlen bzw. durchgelassene Strahlen darstellen. Lassen Sie die Position des einzelnen Schritts bei sein X = 0 anstatt der L / 2 im Diagramm dargestellt.

ψ ( X ) = { A Ö e ich k 1 X + B e ich k 1 X für  X < 0 C e ich k 2 X + D e ich k 2 X für  X > 0
Wo
( k 2 ) 2 = 2 M 2 ( E + | v Ö | )

A Ö stellt den einfallenden Strahl dar, B stellt den reflektierten Strahl dar, und C stellt den übertragenen Strahl dar, und D = 0 . Als Randbedingungen benötigen wir die Wellenfunktion ψ und sein Derivat D ψ D X kontinuierlich sein bei X = 0 , was dazu führt:

A Ö + B = C
Und
ich k 1 ( A Ö B ) = ich k 2 C

Die obigen Bedingungen erlauben uns, Koeffizienten zu bestimmen B Und C in Bezug auf unsere bekannt A Ö :

B = A Ö ( k 1 k 2 ) k 1 + k 2
C = A Ö 2 k 1 k 1 + k 2

Beachten Sie, dass der zeitliche Teil von Ψ ( X , T ) geht in der Gleichung für auf Null J , wir gebrauchen ψ ( X ) = A Ö e ich k X um den Wahrscheinlichkeitsstrom des einfallenden Strahls zu finden, J ich N C

J ich N C = ich 2 M { A Ö e ich k 1 X ( A Ö ich k e ich k 1 X ) A Ö e ich k 1 X A Ö ich k e ich k 1 X }

= 2 ich 2 | A Ö | 2 k 1 2 M

J ich N C = | A Ö | 2 k 1 M

Folglich, J R e F l Und J T R A N S sind ähnlich gegeben durch:

J R e F l = | B | 2 k 1 M
J T R A N S = | C | 2 k 2 M = ( A Ö 2 k 1 k 1 + k 2 ) 2 k 2 M

Schließlich der Transmissionskoeffizient T finden Sie unter

T = | J T R A N S | | J ich N C | = ( A Ö 2 k 1 k 1 + k 2 ) 2 k 2 | A Ö | 2 k 1 = 4 k 1 k 2 ( k 1 + k 2 ) 2


Ich entschuldige mich für die Länge bis zu diesem Punkt, aber jetzt frage ich mich, ob ich einen Übertragungskoeffizienten schätzen kann T für ein lineares Potential (mit E > 0 ) wie im Diagramm, indem mehrere Stufenpotentiale entlang der Länge verwendet werden L . Das Einzelschrittpotential ist keine gute Annäherung, da es keine gibt L Abhängigkeit. Wenn ich also jetzt das durch die rote Linie im Diagramm angegebene Potential untersuche, erhalte ich

v ( X ) = { 0 für  X < 0 | v Ö | / 3 für  0 < X < L / 3 2 | v Ö | / 3 für  L / 3 < X < 2 L / 3 | v Ö | für  2 L / 3 < X < L X > L

Wenn ich jetzt Randbedingungen auferlege ψ Und D ψ D X Bei jeder der Abwärtsfunktionen wiederhole ich im Wesentlichen das oben Gezeigte, aber dreimal mit drei verschiedenen Übertragungsstrahlen, und wenn ich sehe, wie sich der Übertragungskoeffizient zwischen den einzelnen Bereichen ändert, kann ich einen Ausdruck für den Übertragungskoeffizienten selbst finden.

Gibt es einen einfacheren Weg, dies zu tun? Offensichtlich, wenn ich eine unendliche Anzahl von Stufenpotentialen durch Aufbrechen modelliere L Und v Ö in unendliche Teile, würde dies eine ideale Näherung für ergeben T . Gibt es eine Möglichkeit, dies mathematisch/numerisch oder sogar grafisch zu tun?

Antworten (2)

Haben Sie darüber nachgedacht, eine analytische Lösung zu finden ?

Im Intervall 0 X L SE ist:

2 2 M ψ | v 0 | X L ψ = E ψ
ψ + 2 M 2 [ E + | v 0 | X L ] ψ = 0
Diese Art von DE hat exakte Lösungen mit Airy-Funktionen A ich Und B ich aber da letzteres gegen unendlich strebt, nimm an C 2 = 0 .

Dann mit der Randbedingung ψ ( 0 ) = A 0 es könnte möglich sein, zu finden ψ In X = L und damit der Übertragungsfaktor.

Könnten Sie erläutern, wie ich die Randbedingungen für die von Ihnen beschriebene Methode verwenden würde? Nach meinem Verständnis erfordert der Übertragungskoeffizient das Verhältnis von einfallendem und übertragenem Wahrscheinlichkeitsstrom. Also sagen wir, ich habe die analytische Lösung für ψ einbeziehen A ich , würde ich dann berechnen J unter Verwendung der Grenzen auferlegt bei X = 0 , L ? Es ist nur nicht klar, wie ich es bekommen würde T von dem, was Sie im letzten Satz sagen

Wer es genau machen will, braucht Airy-Funktionen für die Wellenlösungen im linearen Potentialbereich und findet die Reflexions- und Transmissionskoeffizienten durch Angleichen der Wellenfunktionen und Ableitungen an den Knotenpunkten. Eine Näherungsmethode wäre die Wentzel-Kramers-Brillouin (WKB)-Näherung.

Siehe https://en.wikipedia.org/wiki/WKB_approximation

Ich bin daran interessiert, die WKB-Näherung auszuprobieren, bin mir aber nicht sicher, wie sie auf meine Frage angewendet werden würde. Ich habe mir den Abschnitt über Schrödinger-Gleichung angesehen. Wie könnte ich für mein Problem ansetzen?
Ist die Methode auch ähnlich der, die @Gert in seiner Antwort beschrieben hat?
@bleuofblue - Ich habe mir die Antwort von Gert angesehen. Er erwähnt auch die Airy-Funktionen, die orthogonale Lösungen der Schrödinger-Gleichung für den Fall eines linearen Potentials sind, genau wie sin und cos (oder Exponentiale) für ein konstantes Potential sind. Diese Lösungen wären exakt. Die WKB-Methode ist eine Näherung für langsam variierende Potentiale, die in QM-Lehrbüchern beschrieben sind.