Dreiecksbarriere in unendlicher Potentialmulde

Angenommen, ich möchte die Wellenfunktion für das folgende 1D-Potential lösen:

(1) U ( x ) = { v 0 a | x | a zum | x | < a zum | x | > a
Da unser Potential symmetrisch ist, haben wir gerade und ungerade Lösungen und können das System damit lösen x 0 und konstruieren Sie anschließend die vollständige Lösung für beide ψ e v e n oder ψ Ö d d . Für die Region innerhalb des Brunnens können wir die zeitunabhängige Schrödinger-Gl. in eine Form, die die Airy-Differentialgleichung ist (beachten Sie, dass es sich um ein lineares Potential handelt), nämlich
(2) d 2 ψ d z 2 z ψ = 0 ,
wo
(3) z Q ( 1 η ) ,

wo

(4) η = x a 0 ,       Q = 2 m v 0 a 0 2 a ,       a 0 a a v 0 E .

Daher können wir ausdrücken ψ wie

(5) ψ ( z ) = C 1 EIN ich ( z ) + C 2 B ich ( z ) .

Unsere Randbedingungen sagen uns:

(6) ψ ( ζ a ) = 0
(7) ψ ( ζ 0 ) = 0 zum ψ Ö d d
(8) ψ ' ( ζ 0 ) = 0 zum ψ e v e n .

wo ich bezeichne ζ 0 = z | x = 0 und ζ a = z | x = a .

Frage : Ist dieses System exakt lösbar?

Typischerweise haben wir bei dieser Art von Problemen nur eine Null-Randbedingung, die es uns ermöglicht, die Energieniveaus durch die Nullstellen der Airy-Funktion zu quantisieren. Aber hier müssen wir gleichzeitig zwei Null-Randbedingungen erfüllen. Zur Verdeutlichung hier unser lineares Gleichungssystem.

ψ Ö d d :

(9) C 1 EIN ich ( ζ 0 ) + C 2 B ich ( ζ 0 ) = 0
(10) C 1 EIN ich ( ζ a ) + C 2 B ich ( ζ a ) = 0

ψ e v e n :

(11) C 1 EIN ich ' ( ζ 0 ) + C 2 B ich ' ( ζ 0 ) = 0 ,
(12) C 1 EIN ich ( ζ a ) + C 2 B ich ( ζ a ) = 0.

Meine vorgeschlagene Idee wäre, die homogenen linearen Gleichungen nach der Nulldeterminante zu lösen, aber dies würde mich zu einer numerischen Lösung zwingen. Ich möchte die Energien in Bezug auf die Nullstellen der Airy-Funktionen finden. Ideen?

Ich vertraue darauf, dass Sie den Artikel Airy-Gleichung in Wikipedia kennen. Der Punkt x = a wo Sie möchten ψ = 0 , fällt auf z = Q ( 1 a / a 0 ) . Es wird Ihnen das Verhältnis geben C 2 / C 1 . Aber achten Sie darauf, wie x sinkt, z steigt. Sie können sehen, wie die Region aussieht z > 0 in den Grafiken im Artikel. Sie werden daran interessiert sein, auch das zu verlangen ψ oder seine Ableitung, bei Null sein x = 0 dh z = Q , das ist links von z = Q ( 1 a / a 0 ) . Es scheint nicht, dass Sie in dieser Region eine weitere Null oder ein Maximum (oder Minimum) für erhalten ψ ( z ) .(Ich fahre fort)
Im Allgemeinen ist ein System aus 2 homogenen Gl. mit 2 Variablen, möglicherweise inkompatibel. Was die Situation retten konnte, waren der Wert von E , die den Wert bestimmt a 0 . Aber auch hier sehe ich ein Problem darin, dass die Werte von z fallen in die Region z > 0 . Es scheint mir, dass Sie Lösungen Ihrer Gleichung erhalten könnten, einschließlich der Bedingung der Symmetrie oder Antisymmetrie, wenn z nahm Werte in der Region z < 0 wo die Funktionen EIN ich und B ich oszillieren. Aber mit z > 0 Ich bin ziemlich skeptisch.

Antworten (2)

Tatsächlich führt kein Weg daran vorbei, dieses Problem numerisch zu lösen. Abweichend von der ursprünglich geposteten Determinantenmethode wäre hier vielleicht ein alternativer Ansatz, das System nur mit der Energie zu lösen. Wir beginnen mit der vollständigen Form unserer Points of Interest am Beispiel der ungeraden Wellenfunktionen:

z ( 0 ) = x 0 a ( E v 0 + 1 )
z ( a ) = x 0 a ( E v 0 ) .
Ich bezeichne a ich und β ich als die ich te Nullen von EIN ich ( z ) und B ich ( z ) beziehungsweise. Ebenfalls a ich ' und β ich ' die Nullstellen der Ableitungen. Ebenfalls x 0 [ 2 m v 0 2 a ] 1 / 3 - Notation des ursprünglichen Beitrags wurde geändert :(

Erstens wissen wir das E > v m ich n und somit E > 0 . Daher sehen wir das z ( 0 ) > z ( a ) . Nun, die erste Null von beiden EIN ich ( z ) oder B ich ( z ) muss bei auftreten z < 0 . Für einen angemessenen Wert von E , können wir das sehen, wenn sich die Nullstellen der Airy-Funktionen um unterscheiden x 0 a , dann wäre unsere Bedingung für eine einzige Energie erfüllt. Wenn wir mit den Airy-Funktionen der ersten Art arbeiten, können wir die Bedingung auferlegen

a ich = x 0 a ( E v 0 )
oder
E = v 0 x 0 a a ich .
Das Kriterium auf der zweiten Null ist dann
a j = x 0 a + a ich .
Wir könnten alle Nullstellen der Airy-Funktion numerisch durchlaufen und für diejenigen mit einer Trennung von
a j a ich = x 0 a
wir hätten dann eine entsprechende Energie von
E = v 0 x 0 a a ich .
Wenn wir dasselbe Konzept sowohl auf gerade als auch auf ungerade Wellenfunktionen anwenden, haben wir:
ψ e v e n ( z ) = C 1 EIN ich ( z ) ; E = v 0 x 0 a a ich ; a j ' a ich = x 0 a
ψ Ö d d ( z ) = C 2 EIN ich ( z ) ; E = v 0 x 0 a a ich ; a j a ich = x 0 a
ψ e v e n ( z ) = C 3 B ich ( z ) ; E = v 0 x 0 a β ich ; β j ' β ich = x 0 a
ψ Ö d d ( z ) = C 4 B ich ( z ) ; E = v 0 x 0 a β ich ; β j β ich = x 0 a
wobei die dritte Spalte die Bedingung bezeichnet, die unsere Nullstellen erfüllen müssen, C ich die geeignete Normalisierungskonstante für jede Wellenfunktion ist, und
z [ 2 m v 0 2 a ] 1 / 3 ( a 0 x ) = x 0 ( a 0 x ) ; a 0 [ a a v 0 E ] .
Wenn wir uns einen Graphen der Airy-Funktion ansehen, können wir das bei sehen z 0 die ``Wellenlänge'' ist sehr klein und daher sollte es eine große Anzahl möglicher Nullstellen geben, die unsere Kriterien erfüllen. Im gleichen Licht, bei dieser Reichweite von z , unsere Energien werden sehr groß sein und wir haben den Anschein des unendlichen quadratischen Brunnens. Danke Sofia und Ali Moh für die hilfreichen Kommentare.

Was Sie verlangen, ist unmöglich. Die Energieeigenwerte sind durch die Nullstellen der Determinante gegeben (wie Sie vermutet haben) und können nicht auf die Nullstellen der Airy-Funktion bezogen werden.

EIN ich ( ζ a ) B ich ( ζ 0 ) EIN ich ( ζ 0 ) B ich ( ζ a ) = 0
Für die ungerade Lösung (und mit Ableitungen an EIN ich für die gerade Lösung).

Seien Sie auch vorsichtig, Ihr z ist nicht dimensionslos, wie es sein sollte, also ist Ihre Änderung der Variablen nicht ganz richtig ...

Bitte lesen Sie die Kommentare für physikalische Argumente, die unendlich viele Lösungen für diese Bedingung garantieren!

Ich bedauere, dass Sie meine Kommentare nicht gelesen haben. Durch die Determinante versucht das OP, die Konstanten zu finden C die es ermöglichen, die Randbedingungen zu erfüllen. Es ist nicht sicher, ob die Determinante zu Null gemacht werden kann . Allerdings, wenn die Determinante zu Null gemacht werden kann, aus dem Wert a 0 für die die Determinante Null wird, kann man die Energie finden E .
Ich habe Ihren Kommentar gelesen, aber ich stimme ihm nicht zu. Die Determinante kann immer zu Null gemacht werden, die Energie kann kontinuierlich geändert werden, um diese Bedingung zu erfüllen. Eine Variable mit einer Bedingung. Aber meine Garantie stammt hauptsächlich aus der Physik, es ist unmöglich, dass es nicht viele Energieeigenwerte gibt! weil dies ein Potenzialtopf mit unendlicher Höhe ist, so dass die Wellenfunktion für ausreichend hohe Energien tatsächlich blind für die dreieckige Erhebung ist, und wir wissen in diesem Fall, dass es eine unendliche Anzahl von diskreten Energieeigenwerten gibt
Daher kann ich nicht sicher sagen, ob diese Eigenwerte vergleichbar sind v 0 (die Höhe der Beule), aber ich kann aus dem obigen physikalischen Argument garantieren, dass es mindestens Eigenzustände mit gibt E >> v 0 .. Außerdem scheint es mir sehr unwahrscheinlich, dass es keine tief liegenden Energieeigenwerte gibt
Ein äquivalenter Weg besteht darin, das asymptotische trigonometrische Verhalten luftiger Funktionen zu beachten, das die physikalische Intuition bestätigt
Es besteht keine Notwendigkeit E >> v 0 , scheint ausreichend E > v 0 . Das Problem ist für E < v 0 . Sie betrachten die Determinante und ich betrachte, wie die Wellenfunktionen in der Region aussehen z > 0 .
Welche Asymptotik? Der Träger der Wellenfunktion ist endlich.
Du verstehst das v ( x ) zum | x | > a , Rechts? also irgendwelche E > 0 stellt einen gebundenen Zustand dar, der durch die lineare Überlagerung luftiger Funktionen ausgedrückt werden kann. Also alles, was ich sage, ist, dass DIE BESTIMMENDE LÖSUNG GARANTIERT IST, zumindest für E >> v 0 ... Ich verstehe nicht, was daran so besonders sein soll E < v 0 dass Sie sagen, es gibt ein Problem? Vergessen z , zeichnen Sie das Potential als Funktion von x und denken Sie an gebundene Zustandswellenfunktionen
Sie bemühen sich sehr, mich davon zu überzeugen, dass es eine Lösung für geben kann E > v 0 dort, wo es einfach ist. Was ist das Besondere an E < E 0 Das ist in diesem Fall schwierig. Und es ist genau wegen dem Verhalten in der Region schwierig z > 0 . Ich glaube, ich habe mich erklärt. Jedenfalls breche ich die Argumentation hier ab.
In Bezug auf die Änderung der Variablen, die in meinem Beitrag bearbeitet wurden, glaube ich nicht, dass sie korrekt sind. Basierend auf dem, was ich vorher hatte: d 2 d x 2 = d 2 d z 2 [ 2 m v 0 2 a ] 2 / 3 , wo z [ 2 m v 0 2 a ] 1 / 3 ( a 0 x )
@Sofia Entschuldigung, an dieser Stelle kann ich nicht verstehen, was Sie sagen. Auch in Bezug auf das asymptotische Verhalten ist klar, dass ich gemeint habe a 0 , was bedeutet es also "die Unterstützung der Wellenfunktion ist endlich"?