QM-Streuprobleme 

1D am Kastenpotential 

Parameter: 

m = 1, hbar = 1, V[1] = 1, V[2] = -1, V[3] = -50 (1)
 

1) Zeitunabhängige Schrödinger-Gleichung (ZSG) 

Kasten-Potential: 

V = piecewise(`<`(abs(x), 1), V[0], 0) (1.1)
 

Graphische Darstellung des Potentials für V_0=1 

 

V = piecewise(`<`(abs(x), 1), -1, 0)
Plot_2d
 

Zeitunabhängige Schrödinger-Gleichung: 

`*`(E, `*`(u(x))) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(hbar, 2), `*`(diff(diff(u(x), x), x)))), `*`(m))), `*`(V, `*`(u(x)))) (1.2)
 

oder mit Potential 

`*`(E, `*`(u(x))) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(hbar, 2), `*`(diff(diff(u(x), x), x)))), `*`(m))), `*`(piecewise(`<`(abs(x), 1), V[0], 0), `*`(u(x)))) (1.3)
 

2) Lösung der ZSG 

2.1) Ansatz 

x<-1: einlaufende und reflektierte Ebene Welle 

|x|<1: Resonanzanteil der Welle 

x>1: auslaufende (transmittierte) Welle 

u(x) = piecewise(`<`(x, -1), `+`(exp(`*`(I, `*`(kappa, `*`(x)))), `*`(R, `*`(exp(`+`(`-`(`*`(`+`(I), `*`(kappa, `*`(x))))))))), `and`(`<`(-1, x), `<`(x, 1)), `+`(`*`(A, `*`(exp(`*`(I, `*`(k, `*`(x))))... (2.1.1)
 

Einsetzen in die ZSG 

`*`(E, `*`(piecewise(`<`(x, -1), `+`(exp(`*`(I, `*`(kappa, `*`(x)))), `*`(R, `*`(exp(`+`(`-`(`*`(`+`(I), `*`(kappa, `*`(x))))))))), `and`(`<`(-1, x), `<`(x, 1)), `+`(`*`(A, `*`(exp(`*`(I, `*`(k, `*`(x...
`*`(E, `*`(piecewise(`<`(x, -1), `+`(exp(`*`(I, `*`(kappa, `*`(x)))), `*`(R, `*`(exp(`+`(`-`(`*`(`+`(I), `*`(kappa, `*`(x))))))))), `and`(`<`(-1, x), `<`(x, 1)), `+`(`*`(A, `*`(exp(`*`(I, `*`(k, `*`(x...
(2.1.2)
 

Wellenzahlen 

für x>1 

`*`(E, `*`(T, `*`(exp(`*`(I, `*`(kappa, `*`(x))))))) = `+`(`/`(`*`(`/`(1, 2), `*`(`^`(hbar, 2), `*`(T, `*`(`^`(kappa, 2), `*`(exp(`*`(I, `*`(kappa, `*`(x))))))))), `*`(m))) (2.1.3)
 

für |x|<1 

`*`(E, `*`(`+`(`*`(A, `*`(exp(`*`(I, `*`(k, `*`(x)))))), `*`(B, `*`(exp(`+`(`-`(`*`(`+`(I), `*`(k, `*`(x))))))))))) = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(hbar, 2), `*`(`+`(`-`(`*`(A, `*`(`^`(k, 2), `*`... (2.1.4)
 

Auflösen nach kappa und k 

kappa = `/`(`*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2)))), `*`(hbar)), k = `/`(`*`(`^`(`+`(`*`(2, `*`(E, `*`(m))), `-`(`*`(2, `*`(V[0], `*`(m))))), `/`(1, 2))), `*`(hbar)) (2.1.5)
 

und einsetzen in das Quantenfeld 

u(x) = piecewise(`<`(x, -1), `+`(exp(`/`(`*`(I, `*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2)), `*`(x)))), `*`(hbar))), `*`(R, `*`(exp(`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`^`... (2.1.6)
 

2.2) Bestimmung der Koeffizienten B, A, R, T 

Stetigkeit links 

`+`(limit(u(x), x = -1, left), `-`(limit(u(x), x = -1, right))) = 0 (2.2.1)
 

also für obiges Quantenfeld 

`+`(exp(`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2))))), `*`(hbar))))), `*`(R, `*`(exp(`/`(`*`(I, `*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2))))), `*`... (2.2.2)
 

Stetigkeit rechts 

`+`(limit(u(x), x = 1, left), `-`(limit(u(x), x = 1, right))) = 0 (2.2.3)
 

also für obiges Quantenfeld 

`+`(`*`(A, `*`(exp(`/`(`*`(I, `*`(`^`(2, `/`(1, 2)), `*`(`^`(`+`(`-`(`*`(m, `*`(`+`(`-`(E), V[0]))))), `/`(1, 2))))), `*`(hbar))))), `*`(exp(`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`^`(`+`(... (2.2.4)
 

Differenzierbarkeit links 

`+`(limit(diff(u(x), x), x = -1, left), `-`(limit(diff(u(x), x), x = -1, right))) = 0 (2.2.5)
 

also für obiges Quantenfeld 

`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`+`(`-`(`*`(exp(`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2))))), `*`(hbar))))), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2...
`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`+`(`-`(`*`(exp(`+`(`-`(`/`(`*`(`+`(I), `*`(`^`(2, `/`(1, 2)), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2))))), `*`(hbar))))), `*`(`^`(`*`(E, `*`(m)), `/`(1, 2...
(2.2.6)
 

Differenzierbarkeit rechts 

`+`(limit(diff(u(x), x), x = 1, left), `-`(limit(diff(u(x), x), x = 1, right))) = 0 (2.2.7)
 

also für obiges Quantenfeld 

`/`(`*`(I, `*`(`^`(2, `/`(1, 2)), `*`(`+`(`*`(exp(`/`(`*`(I, `*`(`^`(2, `/`(1, 2)), `*`(`^`(`+`(`-`(`*`(m, `*`(`+`(`-`(E), V[0]))))), `/`(1, 2))))), `*`(hbar))), `*`(`^`(`+`(`-`(`*`(m, `*`(`+`(`-`(E),... (2.2.8)
 

2.3)Lösung des Gleichungssystems 

Die analytische Formen Koeffezienten B.A.R.T sind sehr komplex. Es macht wenig Sinn, sie hier anzugeben. Sie werden aber vom Maple berechnet im Folgenden verwendet. 

3) Streuphase Delta 

Die Streuphase ist die Phase zwischen einlaufender und transmiitierten Welle. Sie enthält Informaiton über das Potential 

delta = arctan(Im(T), Re(T)) (3.1)
 

 

Graphische Darstellung für verschiedene Potentialtiefen 

 

400

Plot_2d Plot_2d Plot_2d

 

Vor allem sieht man, ob das Potential repulsiv oder attraktiv ist und ob es Bindungszustände oder Resonanzen gibt. 

4) Reflektions- und Transmissionskoeffizienten 

Reflektions- und Transmissionskoeffizienten sind die Wahrscheinlichkeiten für Reflektion oder Transmission 

`*`(conjugate(R), `*`(R)), `*`(conjugate(T), `*`(T)) (4.1)
 

Graphische Darstellung 

Plot_2d Plot_2d Plot_2d

 

Beim repulsiven Potential beginnt die Kurve mit Reflektion. Es gibt jedoch auch schon für Energien kleiner als V0 eine Tranmission, ein sogenanntes Tunneln.  

Beim attraktiven Potential geht die Kurve sehr schnell in Transmission über. 

Für große Energie wird das Potential in allen Fällen eher bedeutungslos. 

5) Dichte 

Dichte 

rho = `*`(conjugate(u(x)), `*`(u(x))) (5.1)
 

Plot_2d Plot_2d Plot_2d

 

Die Wahrscheinlichkeitsdichte zeigt in Abhängigkeit der Energie und Potentialtiefe ein sehr komplexes Verhalten