Quantenmechanik
Quantenmechanik
Quantenmechanik
Wintersemester 2006/07
Matthias R. Gaberdiel
Institut fu
r Theoretische Physik
ETH-Honggerberg
CH-8093 Zu
rich
Email: [email protected]
Inhaltsverzeichnis
1 Historische Anf
ange
1.1 Das Plancksche Strahlungsgesetz (1900) .
1.2 Der Photoeffekt (1905) und Comptoneffekt
1.3 Die Bohrsche Quantenhypothese (1913) . .
1.4 Bohr-Sommerfeld Quantisierung (1915) . .
1.5 Teilchen als Welle (de Broglie 1923) . . . .
. . . .
(1923)
. . . .
. . . .
. . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
2 Wellenmechanik
2.1 Die Schrodinger Gleichung . . . . . . . . . . . . . . . . .
2.2 Wahrscheinlichkeitsstrom und die Kontinuitatsgleichung .
2.3 Das Ehrenfestsche Theorem . . . . . . . . . . . . . . . .
2.4 Realitat der physikalischen Observablen . . . . . . . . . .
2.5 Zeit-unabhangige Schrodinger Gleichung . . . . . . . . .
2.6 Energie Eigenzustande . . . . . . . . . . . . . . . . . . .
2.7 Energiemessung . . . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
5
5
8
9
11
13
.
.
.
.
.
.
.
14
14
16
17
18
19
21
22
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
26
26
29
32
35
37
38
4 Der
4.1
4.2
4.3
4.4
4.5
4.6
4.7
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
42
42
43
44
45
47
49
51
5 Die
5.1
5.2
5.3
Heisenbergsche Unsch
arferelation
53
Nicht-vertauschende Observable . . . . . . . . . . . . . . . . . . . . . . . . 53
Die Unscharfe einer Observablen . . . . . . . . . . . . . . . . . . . . . . . . 56
Die Heisenbergsche Unscharferelation . . . . . . . . . . . . . . . . . . . . . 58
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
61
61
62
63
66
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
68
68
69
71
72
75
78
79
80
.
.
.
.
.
82
82
83
85
85
88
8 Das Wasserstoffatom
8.1 Relativkoordinaten . . .
8.2 Coulomb-Potential . . .
8.2.1 Das Verhalten bei
8.2.2 Das Verhalten bei
8.3 Dynamische Symmetrie .
.
.
r
r
.
. . . .
. . . .
=0 .
. . . .
.
.
.
.
.
.
.
.
.
2
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
9 Drehimpulsaddition
9.1 Tensorprodukte . . . . . . . . . . . . . . . . .
9.2 Addition von Drehimpulsen . . . . . . . . . .
9.3 Addition zweier j = 12 Darstellungen . . . . .
9.4 Die Clebsch-Gordon Reihe im allgemeinen Fall
9.5 Clebsch-Gordon Koeffizienten . . . . . . . . .
9.5.1 Ein einfaches Beispiel . . . . . . . . . .
9.6 Physikalische Beispiele . . . . . . . . . . . . .
9.6.1 Magnetisches Moment . . . . . . . . .
9.7 Der anormale Zeeman Effekt . . . . . . . . . .
9.7.1 Spin-Bahn-Kopplung . . . . . . . . . .
10 Quantenmechanik und klassische
10.1 Das EPR Paradox . . . . . . . .
10.2 Die Bellsche Ungleichung . . .
10.3 Quanten Teleportation . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Physik
. . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . .
11 St
orungstheorie
11.1 Nicht-entartete zeit-unabhangige Storungstheorie . . . . . . . . .
11.1.1 Gestorter harmonischer Oszillator . . . . . . . . . . . . . .
11.2 Entartete zeit-unabhangige Storungstheorie . . . . . . . . . . . . .
11.2.1 Aufhebung einer zwei-fachen Entartung zu erster Ordnung
11.2.2 Der allgemeine Fall . . . . . . . . . . . . . . . . . . . . . .
11.2.3 Aufhebung der Entartung zu zweiter Ordnung . . . . . . .
11.2.4 Der Stark Effekt . . . . . . . . . . . . . . . . . . . . . . .
11.3 Zeit-abhangige Storungstheorie . . . . . . . . . . . . . . . . . . .
11.3.1 Der Propagator der Quantenmechanik . . . . . . . . . . .
11.3.2 Das Heisenberg-Bild . . . . . . . . . . . . . . . . . . . . .
11.3.3 Das Wechselwirkungsbild . . . . . . . . . . . . . . . . . . .
11.3.4 Erste Ordnung Storungstheorie . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
91
91
92
93
96
97
99
100
100
102
103
106
. 106
. 107
. 109
111
. 111
. 114
. 115
. 115
. 118
. 119
. 120
. 121
. 121
. 124
. 125
. 126
Dieses Skript basiert zum Teil auf den Vorlesungsskripten von Gianni Blatter und Gian
Historische Anf
ange
Ende des 19. Jahrhunderts basierte das physikalische Weltbild auf dem, was wir heute die
klassische Physik nennen: die wesentlichen Grundpfeiler waren die klassische Mechanik
(`a la Newton), die Elektrodynamik (`a la Maxwell) und die Thermodynamik (`a la Boltzmann). Diese Vorstellung wurde zu Beginn des 20. Jahrhunderts durch eine Reihe von
Experimenten in Frage gestellt, die wir im Folgenden kurz skizzieren wollen.
1.1
Nach klassischer Vorstellung ist ein Teilchen (z.B. ein Elektron oder ein Atom) durch Ort
und Geschwindigkeit charakterisiert, deren Angabe beliebig genau sein kann; eine Welle hat eine bestimmte Frequenz und Wellenzahl, oder ist eine Superposition von solchen.
Licht besteht aus elektromagnetischen Wellen und Materie aus Teilchen. Die beiden wechselwirken miteinander und gelangen, in einem verspiegelten Hohlraum eingeschlossen, zu
einem thermischen Gleichgewicht, das Planck untersuchte.
Das zunachst freie elektromagnetische Feld gen
ugt auf Grund der Maxwell-Gleichungen der Wellengleichung
~ =0
2E
und der Nebenbedingung
mit 2 =
1 2
c2 t2
~ =0,
div E
(1.1.1)
(1.1.2)
~ Der Separationsansatz
(analog f
ur B).
~ x, t) = f (t)E(~
~ x)
E(~
f
uhrt auf
1 ~
~ x)
f (t)E(~x) = f (t)E(~
c2
und weiter
f = 2 f
| {z }
2
~ = E
~
E
2
c
{z
}
|
harmonischer Oszillator
(1.1.3)
Eigenwertgleichung f
ur
f
ur eine Konstante 2 ( 0). Auf dem Rand des Hohlraums gelten die Randbedingungen
~
E
=0.
n
~k = 0 ,
E
ki = ni ,
ni ganz, 0, hochstens ein ni = 0.
(1.1.4)
L
~ ~k = 0, wobei E
~ = (E1 , E2 , E3 )
Wegen der Nebenbedingung (1.1.2) muss weiterhin gelten E
und ~k = (k1 , k2 , k3 ). Zu jedem ~k gibt es daher zwei linear unabhangige Eigenschwingungen
mit Eigenfrequenzen = c |~k|. Die Zahl der Eigenschwingungen ist nach (1.1.4)
asymptotisch
3
3
L
V
1 4
= 2 3
N() = 2
8 3
c
c
3
f
ur grosse , wobei V = L3 , bzw.
2
dN
=V 2 3 .
d
c
(1.1.5)
[Der Faktor 2 beschreibt die beiden Losungen und der Faktor 1/8 tritt auf, da nur ni 0
beitragen. Schliesslich ist das Volumen der 3d Einheitskugel 4r 3 /3, wobei der relevante
Radius r = L/c ist es gilt ~k 2 = 2 /c2 und ~k = L ~n.]
Planck stellte sich nun die Materie als aus Oszillatoren (Resonatoren) aller Frequenzen 0 bestehend vor, welche die sonst unabhangigen elektromagnetischen Schwingungen
ins Gleichgewicht bringen. Bei gegebener Temperatur T sind die Resonatoren termisch
angeregt. In der Tat ist nach Boltzmann die Wahrscheinlichkeit w(p, q), ein Hamiltonsches
System mit Phasenkoordinaten p, q bei der Temperatur T in dpdq zu finden gerade
w(p, q) dp dq =
eH(p,q)
dp dq ,
Z()
(1.1.6)
wobei H die Hamiltonfunktion ist, = (kT )1 die inverse Temperatur, k die BoltzmannKonstante und
Z
Z() = dp dq eH(p,q) .
(1.1.7)
E = dp dq H(p, q)w(p, q) =
log Z() .
F
ur einen 1-dimensionalen harmonischen Oszillator der Frequenz 0 ist,
H=
1
p2
+ m02 q 2 ,
2m 2
,
0
und damit
= log = kT ,
E
f
ur grosse (Wien, 1896). Planck bemerkte, dass ~ eine neue Naturkonstante sein musste.
u(, T ) 6
Rayleigh-Jeans
Planck
Wien
-
(in Klammern die heutigen Werte). Die grossartigste Bestatigung fand das Plancksche Gesetz (1.1.11) in der gemessenen Spektralverteilung der kosmischen Hintergrundstrahlung
bei T = 2, 73K (COBE 1992, WMAP 2003).
7
(n = 0, 1, 2, ...) .
(1.1.12)
Z() =
en~0 =
1
1 e ~0
~ 0
log Z() = ~0
,
E =
e
1
1.2
Wird eine (metallische) Oberflache durch Licht der (Kreis-)Frequenz bestrahlt, treten
Elektronen aus das ist der sogannante photoelektrische Effekt, der von Hertz 1887
erstmals entdeckt wurde.
00
11
Lenard (1902) bemerkte, dass die Energie T der emittierten Elektronen (monoton wachsend) nur von der Frequenz
abhangt, nicht aber von der Intensitat der einfallenden Strahlung. Davon
abhangig ist hingegen die Emissionsrate.
Metall
11
00
00
11
Licht der Frequenz
00
11
00
11
00
11
00
11
00
11
00
11
1111111
000000000
11
00
11
Elektronen der
00
11
00
11
kinetischen Energie T
00
11
00
11
00
11
Diese Deutung basiert darauf, dass das Licht (das ja durch elektromagnetische Wellen beschrieben wird) auch Teilchenaspekte besitzt, namlich, dass es aus Lichtquanten
zusammengesetzt zu sein scheint, die eine bestimme Energie tragen und mit einzelnen
Elektronen wechselwirken. Dieser Teilchenaspekt des Lichtes wurde spater eindr
ucklich
durch den Comptoneffekt demonstriert, bei dem Lichtstrahlen an Elektronen gestreut
werden:
Abbildung 2: Rontgen-
h k = p
Strahlen (Photonen )
treffen auf ein ruhendes
Elektron e und werden
gestreut.
pe
Dieser Effekt kann dadurch richtig beschrieben werden, dass man den Streuprozess als
Teilchen-Teilchen (Photon-Elektron) Streuung unter Energie- und Impulserhaltung interpretiert. (Da sich das Licht mit Lichtgeschwindigkeit ausbreitet, muss man dazu nat
urlich
k = |~k| = =
c
p = (~k, ~~k) ,
(1.2.2)
gegeben, wobei ~k der Wellenvektor ist, dessen Richtung die Ausbreitungsrichtung der
Welle beschreibt und die Wellenlange ist.
1.3
nm = R
m2 n2
(Balmer 1885, aus den 4 Linien m = 2, n = 3, 4, 5, 6).
Bohr nimmt an (analog zur Planckschen Quantisierung des Resonators, aber gegen
klassische Vorstellungen), dass das Atom nur in Zustanden mit diskreten Energien En
existieren kann. Strahlung (namlich ein Lichtquant) der Frequenz
1
nm = (En Em )
~
1
,
n2
n = 1, 2, . . .
(1.3.2)
mit Ry = R ~.
En
..
Als Modell des Atoms verwendet Bohr das von Rutherford (1911): Ein Elektron (Masse
m, Ladung e) im Feld eines viel schwereren Kerns (Ladung e), den wir zunachst als fest
annehmen. Langs klassischen Bahnen w
urde das Elektron strahlen und so dem Kern
stets naher kommen. Bohr wahlt die Quantenzustande unter den Kreisbahnen (Radius r,
Winkelgeschwindigkeit , Drehimpuls L, Energie E). F
ur diese gilt
mr 2 =
e2
,
r2
L = mr 2 ,
E=
L2
e2
.
2mr 2
r
Daraus folgt
L2
r=
,
me2
me4
E= 2 ,
2L
me4
= 3 .
L
(n = 1, 2, . . .)
(1.3.3)
und findet
rn = a0 n2 ,
a0 =
En = Ry
n =
1
,
n2
~2
me2
me4
Ry =
2~2
(Bohr-Radius) ,
(Rydbergkonstante) ,
2Ry
~n3
(1.3.4)
(1.3.5)
Ry = 13, 6058 eV .
Ber
ucksichtigt man die Mitbewegung des Kerns der Masse M, so ist m durch die reduzierte
Masse zu ersetzen; ebenso seine Ladung e durch Ze bei wasserstoffahnlichen Ionen wie
zum Beispiel He+ ; also Ry in (1.3.4) durch
Z2
M
Ry .
M +m
Die Wahl von ~ als Proportionalitatsfaktor in (1.3.3) ist zwingend, falls man die G
ultigkeit
Anode
Gitter
inelast. St.
I elast. St.
HgGas
0.5V
10
15 V
1.4
p
1
+ m02 q 2 = E ,
2m 2
also
pdq = 2mE
2mE
q
2E 1 2E
=
m 0
0
2E 1
m 0
Dabei sind (1 , . . . , f ) = Erhaltungsgrossen. Im 2f -dimensionalen Phasenraum verlauft die Bewegung dann auf dem Schnitt von f durch bestimmte Flachen
pk =
S
Sk
(q, ) =
(qk , ) ,
qk
qk
(k = 1, . . . , f ) .
(1.4.2)
Falls die f -dimensionalen Schnittflachen topologisch Tori sind, so ist die SommerfeldBedingung anwendbar. Sie zeichnet als erlaubt diejenigen Tori (und nicht spezielle, darin
verlaufende Bahnen) aus, f
ur welche
I
Wk := pk dqk = 2nk ~ ,
(nk = 0, ()1, . . .)
(1.4.3)
f
ur alle k = 1, . . . f , wobei pk durch (1.4.2) gegeben ist. Dies bestimmt (1 , . . . f ) als
Funktion der nk und insbesondere die moglichen Energien En1 ...nk .
12
Auf Systeme, die nicht separabel sind (und dies ist der generische Fall), kann man
Sommerfelds Bedingung nicht anwenden. Sie ist daher nicht wirklich von fundamentaler
Signifikanz.
1.5
Wie wir in Kapitel 1.1.2 gesehen haben, verhalten sich Lichtwellen manchmal so, als ob
sie Teilchen waren. De Broglie hat vorgeschlagen, dass umgekehrt auch Teilchen Wellencharakter besitzen. Konkret postulierte er, dass die Relation (1.2.2) auch f
ur Teilchen gilt,
d.h. dass ein Teilchen mit Impuls p einer Welle mit Wellenlange und Kreisfrequenz
=
h
,
p
E
~
(1.5.1)
entspricht. [E = cp0 .]
Diese Hypothese ist u
berzeugend durch das Davison-Germer Experiment (1927)
bestatigt worden. Dabei betrachtet man ein Streuexperiment an einem Festkorper (`a la
von Laue), aber ben
utzt einen Elektronenstrahl statt einer Lichtwelle. Die Reflexion von
Elektronen von einer Kristalloberflache ergibt von Laue Reflexe
Detektor
einfallender
Strahl
gestreuter
Strahl
13
Wellenmechanik
Ausgehened von dem Postulat von de Broglie wollen wir jetzt (heuristisch) die Schrodinger
Gleichung ableiten.
2.1
Die Schr
odinger Gleichung
Nach de Broglie ist einem Teilchen mit Impuls p = (E, ~p) eine Welle zugeordnet, deren
Wellenlange und Kreisfrequenz durch (1.5.1) gegeben sind
=
2 ~
p
E
~
(2.1.1)
2x
t
(2.1.2)
wobei (0, 0) eine Konstante ist, die die Amplitude der Welle bestimmt. Diese Welle ist
durch die Gleichungen charakterisiert
2
p
(x, t) =
(x, t) = (x, t)
x
~
E
(2.1.3)
(2.1.4)
Im nicht-relativistischen Grenzfall (den wir im folgenden behandeln wollen die quantenmechanische Beschreibung der relativistischen Physik f
uhrt zur Quantenfeldtheorie)
ist die Energie E eines freien Teilchens mit Impuls p durch
E=
1 2
p
2m
(2.1.5)
~2 2
(x, t) = E(x, t) =
(x, t) .
t
2m x2
(2.1.6)
Falls das Teilchen nicht frei ist, sondern sich im Einfluss des Potentials V (x) bewegt, hat
die rechte Seite von (2.1.5) zusatzlich den Term V (x); dann erhalten wir
~2 2
~ (x, t) =
+ V (x) (x, t) .
(2.1.7)
t
2m x2
Diese Gleichung ist die (zeit-abh
angige) Schr
odinger-Gleichung.
Die Funktion (x, t) in (2.1.7) ist eine komplexe Funktion von (x, t). Falls eine Lichtwelle beschreiben w
urde, dann ist die Intensitat des Lichtes bei (x, t) gerade |(x, t)|2 .
14
Im gegenwartigen Kontext interpretieren wir daher (Born 1926) |(x, t)|2 als die Wahrscheinlichkeit mit der das Teilchen sich zur Zeit t am Punkt x befindet. Damit diese
Interpretation Sinn macht muss (x, t) so normalisiert sein, dass
Z
dx|(x, t)|2 = 1
f
ur alle t.
(2.1.8)
R
Damit haben wir bereits die wichtigsten Postulate der Wellenmechanik zusammengetra
gen; zur Ubersicht
stellen wir sie noch einmal zusammen:
(i) Die Wahrscheinlichkeit das Teilchan zur Zeit t am Ort x anzutreffen, ist durch das
Normquadrat |(x, t)|2 der Wellenfunktion gegeben.
(ii) Die Wellenfunktion (x, t) ist f
ur jedes t quadrat-integrabel in x; insbesondere ist
sie wie in (2.1.8) normalisiert.
(iii) Die Wellenfunktion (x, t) ist eine komplexe Funktion, die der zeit-abhangigen
Schrodinger Gleichung gen
ugt
~
(x, t) = H(x, t) .
t
(2.1.9)
2
+ V (x)
2m x2
~2
(2.1.10)
definiert ist. Auch in der Quantenmechanik ist also die Hamiltonfunktion (bzw. der
Hamitlonoperator) die Erzeugende der Zeitentwicklung.
Wir sollten betonen, dass die obige Diskussion diese Postulate nicht ableitet, genausowenig wie man beispielsweise Newtons Kraftgesetz ableiten kann. Wir haben lediglich
versucht, sie zu motivieren. Wir wollen sie jetzt als die Regeln, die Quantenmechanik
definieren, akzeptieren, und studieren, welche Konsequenzen sie haben.
In der Definition des Hamiltonoperators (2.1.10) begegnen wir zum ersten Mal einer speziellen Eigenschaft der Quantenmechanik: klassische Observable (zum Beispiel die
Energie, d.h. die Hamiltonfunktion) werden in der Quantenmechanik durch Operatoren
(in diesem Fall den Hamiltonoperator H) dargestellt. Diese Operatoren wirken auf dem
Raum der Wellenfunktionen. Der Umstand, dass diese Operatoren im allgemeinen nicht
vertauschen, ist der entscheidende Grund f
ur viele Quanteneffekte der Quantenmechanik,
insbesondere f
ur die Heisenbergsche Unscharferelation (die wir ein wenig spater ableiten
werden).
Ein anderes Beispiel f
ur diese Operatorbeschreibung von Observablen ist (2.1.3): auf
Wellenfunktionen ist die Impulsobservable p durch den Differentialoperator
p = ~
15
(2.1.11)
dargestellt.
Die obigen Postulate sind spezifisch f
ur die Situation, wo wir ein Teilchen im Raum
analysieren wollen. Wir werden die Postulate spater in eine abstraktere Form bringen. Wir
sollten auch noch darauf hinweisen, dass es noch ein weiteres Postulat gibt, das erklart,
was genau bei einer Messung geschieht. Wir werden nat
urlich auch auf diesen Punkt spater
zur
uckkommen. F
ur den Moment wollen wir jedoch zunachst ein paar Konsequenzen dieser
Postulate studieren.
2.2
Wie wir oben postuliert haben hat |(x, t)|2 die Interpretation einer Wahrscheinlichkeitsdichte
(x, t) |(x, t)|2 .
(2.2.1)
Die Schrodinger Gleichung (und ihre komplex Kunjugierte) bestimmen die Zeitenwicklung
dieser Wahrscheinlichkeitsdichte:
d
(x, t) = (x, t) (x, t) + (x, t) (x, t)
dt
t
t
2
~
=
(x, t)
(x, t) (x, t) 2 (x, t)
2m
x2
x
(2.2.2)
(2.2.3)
(x, t) +
j(x, t) = 0
(2.2.4)
dt
x
geschrieben werden, wobei j(x, t) der Wahrscheinlichkeitstrom ist,
(x, t) (x, t) (x, t)
j(x, t) =
(x, t) .
(2.2.5)
2m
x
x
Die Gleichung (2.2.4) hat die Intepretation einer Kontinutatsgleichung. Um das genauer
zu verstehen, integrieren wir die Wahrscheinlichkeitsdichte u
ber ein kleines Intervall
Z x0 +
P (x0 ; ) =
dx (x, t) .
(2.2.6)
x0
Dann ist
Z x0 +
d
P (x0 ; ) =
dx j(x, t)
dt
x
x0
= j(x0 , t) j(x0 + , t) .
(2.2.7)
Die Anderung
in der Wahrscheinlichkeit, das Teilchen in dem Intervall zu finden, wird also
dadurch beschrieben, wieviel Wahrscheinlichket aus dem Intervall geflossen oder in das
16
Intervall geflossen ist. d.h. durch die Werte des Wahrscheinlichkeitsstroms an den Intervallgrenzen. Die Wahrscheinlichkeit verhalt sich also wie eine inkompressible Fl
ussigkeit:
sie kann weder vernichtet noch erzeugt werden; sie kann sich in einem Intervall nur dadurch andern, dass sie aus dem Intervall ausfliesst oder in es hineinfliesst.
2.3
Wie schon im vorigen Kapitel bestimmt die Schrodinger Gleichung die Zeitentwicklung
dieses Erwartungswertes. In der Tat berechnen wir
Z
d
hxi =
dx x ((x, t)(x, t) )
dt
t
ZR
~ 2
(x, t) V (x)(x, t) (x, t)
=
dx x
2
2m
x
~
R
2
~
.
(2.3.2)
+(x, t)
(x, t) + V (x)(x, t)
2m x2
~
Die Terme, die das Potential involvieren, k
urzen sich gegenseiteig. Im ersten Term integrieren wir nun partiell und erhalten
x=
Z
2
~
~
dx x (x, t)
x (x, t)
(x, t) =
(x, t)
2m R
x2
2m
x
x=
Z
~
dx (x, t)
(x, t) + x (x, t)
(x, t) .
2m R
x
x
x
Da quadrat-integrabel ist, d.h. da es (2.1.8) erf
ullt, muss der Randterm verschwinden. [Dieses Argument ist ein klein wenig schlampig!] Nun erinnern wir uns, dass der
Impulsopertor durch (2.1.11) definiert ist. Der erste Term in der zweiten Zeile ist also
Z
Z
1
1
~
dx (x, t)
dx (x, t) p (x, t) =
(x, t) =
hpi .
(2.3.3)
2m R
x
2m R
2m
1
Mal den Erwartungswert des Impulses
In Analogie zu (2.3.1) interpretieren wir dies als 2m
p. Den zweiten Term in der zweiten Zeile integrieren wir wiederum partiell
x=
Z
~
dx x (x, t)
x
(x, t) =
(x, t) (x, t)
2m R
x
x
2m
x
x=
Z
~
+
dx (x, t) (x, t) + x (x, t) 2 (x, t) .
2m R
x
x
17
~
~
dx (x, t) (x, t) =
[(x, t)(x, t) ]x=
x=
2m R
x
2m
Z
~
dx (x, t) (x, t)
(2.3.4)
2m R
x
1
hpi ,
(2.3.5)
=
2m
wobei wir wiederum den Randterm weggelassen haben [dieses Mal ist das wirklich korrekt!]
sowie (2.3.3) ben
utzt haben. Insgesamt erhalten wir also
d
1
hxi = hpi .
dt
m
(2.3.6)
Das ist die bekannte (nicht-relativistische) Relation zwischen Geschwindigkeit und Impuls.
d2
hxi = hV (x)i .
dt2
Dieses Resultat wird Ehrenfestsches Theorem genannt. Es stellt eine Beziehung zwischen
Quantenmechanik und klassischer Mechanik her: die Erwartungswerte der Quantenmechanik erf
ullen die klassischen Bewegungsgleichungen.
2.4
Realit
at der physikalischen Observablen
(x, t)
hpi = ~ dx (x, t)
x
R
(2.4.1)
definiert. Der Erwartungswert ist das mittlere Ergebnis einer Impulsmessung. Damit das
eine vern
unftige physikalische Interpretation hat, m
ussen wir verlangen, dass (2.4.1) eine
reelle Zahl ist, d.h. dass
hpi = hpi .
(2.4.2)
18
Wir wollen nun zeigen, dass dies automatisch der Fall ist. (Im wesentlichen ist das dieselbe
Rechnung, die wir oben schon durchgef
uhrt haben.)
Z
(x, t)
hpi = ~ dx (x, t)
x
R
Z
x=
(x, t)
= ~ [(x, t)(x, t) ]x= ~ dx (x, t)
x
R
= hpi ,
(2.4.3)
wobei wir ben
utzt haben, dass die Integrabilitat von (x, t) impliziert, dass der Randterm
verschwindet. Man kann entsprechend nachpr
ufen, dass der Erwartungswert der Energie
Z
hHi =
dx (x, t) H (x, t) ,
(2.4.4)
R
wobei H der Hamilton Operator (2.1.10) ist, auch reell ist (Ubungsaufgabe),
hHi = hHi .
2.5
(2.4.5)
Zeit-unabh
angige Schr
odinger Gleichung
Wir haben bisher immer angenommen, dass das Potential V (x) tatsachlich nicht von der
Zeit t explizit abhangt. Falls das der Fall ist, konnen wir die x und t Variablen separieren;
dazu machen wir den Ansatz
(x, t) = (x)(t) .
(2.5.1)
Die Schrodinger Gleichung impliziert nun, dass
~2
(x) + V (x)(x) ,
(x)~(t)
= (t)
2m
(2.5.2)
d
wobei (t)
=
(x) + V (x)(x) .
(2.5.3)
~
(t)
(x)
2m
Wir beobachten, dass die linke Seite nicht von x abhangt, wohingegen die rechte Seite
nicht von t abhangt. Damit die beiden Ausdr
ucke gleich sein konnen, m
ussen sie beide
gleich der selben Konstanten sein, die wir mit E bezeichnen. Wir haben damit die zeitabhangige Schrodingergleichung in zwei separate Gleichungen zerlegt,
~(t)
= E(t)
H(x) = E(x) ,
19
(2.5.4)
wobei H der Hamiltonoperator (2.1.10) ist. Wir konnen die erste Gleichung leicht losen
E
(2.5.5)
(t) = exp t .
~
Falls also (x) die Eigenwert Gleichung (2.5.4) mit Eigenwert E erf
ullt, dann lost
E
(2.5.6)
(x, t) = (x) exp t
~
(2.5.7)
wird als zeit-unabhangige Schrodinger Gleichung bezeichnet. Da |(x, t)| = |(x)|, ist die
Wellenfunktion korrekt normalisiert (d.h. erf
ullt (2.1.8)), falls (x) die Normalisierungsbedingung as
Z
R
dx |(x)|2 = 1
(2.5.8)
erf
ullt. Offenbar ist die Normalisierungsbedingung (2.1.8) dann f
ur alle Zeiten erf
ullt.
Die zeit-unabhangige Schrodinger Gleichung ist eng mit der zeit-unabhangigen Hamilton-Jakobi Gleichung verwandt, die letztes Jahr in der Mechanik diskutiert wurde. (In
der Tat fand mit solchen Argumenten Schrodinger urspr
unglich die nach ihm benannte
Gleichung!) Um das genauer zu verstehen machen wir den Ansatz (in 3 Dimenisonen)
(~x, t) = A(~x) e
S(~
x)
~
e ~ t ,
(~x) = A(~x) e
S(~
x)
~
(2.5.9)
wobei sowohl A(~x) als auch S(~x) reelle Funktionen sind. Dann ist die zeit-unabhangige
Schrodinger Gleichung f
ur (~x) (wir nehmen an, wie beschreiben ein Teilchen mit Potential V (~x))
(S)2
~2 A
+ (V (~x) E) =
,
2m
2m A
A
A S + S = 0 ,
2
wobei die erste (zweite) Gleichung gerade den Realteil (Imaginarteil) der zeit-unabhangigen Schrodinger Gleichung beschreiben. Im klassischen Limes, d.h. f
ur ~ 0, wird die
erste Gleichung dann
(S)2
+ V (~x) = H(~x, S) = E .
(2.5.10)
2m
Das ist gerade die zeit-unabhangige Hamilton-Jakobi Gleichung. Die Phase der Welle
kann also mit der (dimensionslosen) Grosse S/~ identifiziert werden, wobei S die erzeugende Funktion der kanonischen Transformation ist, die in der Hamilton-Jakobi Theorie
auftritt.
20
Andererseits erinnern wir uns aus der Mechanik Vorlesung daran, dass eine Losung S(~x) der HamiltonJakobi Gleichung (2.5.10) die physikalische Bahn gerade
vermoge
p~ = m~x = S .
(2.5.11)
beschreibt. Die physikalische Bahn ist also immer senkrecht zu der Flache S = const. Diese Flachen sind jedoch
gerade die Wellenfronten der Welle , da S die Phase
der Welle ist. In dieser Weise ist also die quantenmechanische Welle direkt mit der klassischen Bahn verbunden.
2.6
Energie Eigenzust
ande
(2.6.1)
In vielen Fallen ist der Satz der Eigenfunktionen abzahlbar, d.h. wir konnen annehmen,
dass sich die Eigenfunktionen durch eine positive ganze Zahl bezeichnen lassen. Weiterhin
wollen wir annehmen, dass jede quadrat-integrierbare Funktion sich als Linearkombination
dieser n schreiben lasst. (Wir werden auf diese Annahmen spater zur
uckkommen, wenn
wir die Quantenmechanik ein wenig formaler beschreiben werden.)
Angenommen wir kennen die Wellenfunktion zur Zeit t = 0, (x, 0). Wir wollen die
Losung f
ur alle t 0 bestimmen. Nach Annahme konnen wir die Funktion (x, 0) durch
die n ausdr
ucken
X
(x, 0) =
an n (x) ,
n=0
wobei die an im allgemeinen komplexe Koeffizienten sind. Es folgt dann aus (2.6.1), dass
X
En
(x, t) =
an n (x) exp
t
(2.6.2)
~
n=0
f
ur alle t. Nach Konstruktion stimmt (2.6.2) mit der gegebenen Funktion (x, 0) f
ur t = 0
u
berein; (x, t) ist ausserdem eine Losung der zeit-abhangigen Schrodinger Gleichung
da jeder einzelne Term eine Losung ist, und da die Differentialgleichung linear ist. Der
Umstand, dass die Summe zweier Losungen wieder eine Losung definiert, nennt man
manchmal das Superpositionsprinzip; wir werden es immer wieder ben
utzen.
Wir haben hier also das Problem, die allgemeine Losung der Schrodinger Gleichung (zu
vorgegebenen Anfangsbedingungen) zu finden, darauf zur
uckgef
uhrt, den vollstandigen
Satz der normierbaren Losungen von (2.5.7) zu finden. Wir sollten schon an dieser Stelle
darauf hinweisen, dass die Bedingung, dass n quadrat-integrabel ist, im wesentlichen f
ur
die Quantisierung der Eigenwerte verantwortlich sein wird.
21
2.7
Energiemessung
Wir wollen nun annehmen, dass wir einen vollstandigen Satz normalisierter Eigenfunktionen des Hamilton-Operators n , n N gefunden haben, wobei die zugehorigen Eigenwerte
Hn = En n sind. Der Einfachheit halber wollen wir weiterhin annehmen, dass alle Eigenwerte distinkt sind, d.h. dass En 6= Em f
ur n 6= m. Wir wollen zeigen, dass dies impliziert,
dass die normalisierten Eigenfunktionen tatsachlich orthonormal sind, d.h. dass
Z
1 falls m = n
(2.7.1)
dx n (x) m (x) = m,n =
0
falls m 6= n.
R
Dazu beobachten wir, dass
Z
Z
Em dx n (x) m (x) =
dx n (x) H m (x)
(2.7.2)
R
R
Z
Z
~2
x=
dx n (x) m (x) =
dx n (x) m
(x)
[n (x) m (x)]x= +
2m R
2m
2m R
Z
x=
~2
~2
n (x) m (x) x=
dx n (x) m (x) ,
=
2m
2m R
wobei wir wiederum die Randterme weggelassen haben (da n quadrat-integrabel ist).
Zusammen erhalten wir also
Z
Z
Em dx n (x) m (x) =
dx (H n ) m (x)
R
R Z
(2.7.4)
= En dx n (x) m (x) ,
R
(2.7.5)
Dies zeigt, dass alle Eigenwerte Em des Hamilton-Operators reell sind. (Wie wir gleich
sehen werden, impliziert das insbesondere auch, dass der Erwartungswert reell ist.)
22
Da die Eigenwerte reell sind, konnen wir nun die beiden Seiten von (2.7.4) voneinander
abziehen, und erhalten
Z
(Em En ) dx n (x) m (x) = 0 .
(2.7.6)
R
Nach Annahme ist Em 6= En (falls m 6= n), und daher muss das Integral verschwinden;
das beweist dann (2.7.1).
We wir gerade gezeigt haben, konnen wir die allgemeinste Losung der Schrodinger
Gleichung immer in der Form (2.6.2) schreiben. Wir wollen nun den Erwartungswert der
Energie f
ur diese Wellenfunktion berechnen. Wegen (2.4.4) m
ussen wir berechnen
Z
hHi =
(x, t) H (x, t)
R
Z
X
X
(En Em )
am an exp
=
t
dx m (x) H n (x)
=
m=0 n=0
X
n=0
|an |2 En ,
(2.7.7)
dx n (x) (x, t) .
(2.7.8)
R
dx n (x) (x, t) =
m=0
am exp
Em
~
En
= an exp
t
23
Z
t
dx n (x) m (x)
R
(2.7.9)
zusammen mit dem Umstand, dass das Exponential eine reine Phase ist (da En reell ist).
Falls die Eigenwerte nicht alle voneinander verschieden sind, gilt eine ahnliche Formel, die
wir im Rahmen der formaleren Diskussion der Quantenmechanik besprechen werden.
Damit die obige Interpretation Sinn macht, muss gelten
X
n=0
|an |2 = 1 .
(2.7.10)
Hierbei ist A eine Konstatne mit |A| = 1, und wir haben wiederum (wie schon
zuvor) angenommen, dass die Eigenwerte von H alle unterschiedlich sind.
Dieses Axiom ist noch weniger klassisch als die anderen; der Kollaps der Wellenfunktion
als Konsequenz einer Messung ist intuitiv nur schwer verstandlich. Die physikalische Interpretation dieses Axioms ist, dass es unmoglich ist, eine Messung durchzuf
uhren, ohne
dabei das System zu beeinflussen. Das Axiom sagt in genauer Weise, was der Effekt dieser Beeinflussung ist. Man kann die Aussage des Axioms auch noch anders formulieren:
nachdem die Energie des Systems einmal als En gemessen wurde, wird jede zuk
unftige
Energiemessung immer wieder En finden und zwar mit Wahrscheinlichkeit 1. (Hierbei
ist angenommen worden, dass das System zwischendurch nicht andersweitig beeinflusst
worden ist.)
Es gibt ein (offensichtliches) konzeptuelles Problem (das Messproblem der Quantenmechanik), das mit diesem Axiom in Verbindung steht. Dies hangt damit zusammen, dass
es nicht klar ist, wie man prazise definieren soll, was genau eine Messung ist. In dieser
Vorlesung wollen wir nicht versuchen, auf diese Fragen naher einzugehen eine moderne
Diskussion mancher dieser Aspekte findet sich in dem Buch von Isham [I]. Im folgenden
wollen wir die Axiome einfach als die Regeln ansehen, die die Wellenmechanik definieren,
ohne ihre konzeptionellen Probleme zu hinterfragen.
24
25
Bevor wir damit beginnen wollen, die Quantenmechanik ein wenig formaler zu analysieren,
wollen wir zunachst ein Gef
uhl f
ur die Theorie bekommen. Dazu wollen wir die zeitunabhangige Schrodinger Gleichung f
ur ein paar einfache Systeme losen. Die einfachsten
Systeme sind 1-dimensionale Systeme, in denen das Teilchen sich nur in einer Richtung
bewegen kann.
3.1
Wir betrachten ein Teilchen, das sich im Intervall [0, a] frei bewegen kann, aber das Intervall nicht verlassen kann. Durch ein Potential V (x) ausgedr
uckt bedeutet dies, dass
0 0xa
(3.1.1)
V (x) =
sonst.
Da das Potential ausserhalb der Box unendlich ist, muss ausserhalb der Box (x) = 0,
x 6 [0, a] gelten. Da das Teilchen im Innern der Box, d.h. f
ur x [0, a], frei ist, m
ussen
wir dort die Differentialgleichung losen
d2
= E .
2m dx2
~2
(3.1.2)
Weiterhin muss u
berall (und insbesondere an den Randpunkten x = 0, a) stetig sein;
falls das nicht der Fall ware, dann ware die zweite Ableitung von bei x = 0, a proportional zu der Ableitung der Deltafunktion (mehr dazu spater), und w
urde daher nicht die
Schrodinger Gleichung losen. Da aber ausserhalb des Intervalls verschwindet gilt also
(0) = (a) = 0 .
Die allgemeine Losung von (3.1.2) ist
p
p
x
x
2m|E|
A
cosh
2m|E|
+
B
sinh
falls E < 0
~
~
A + Bx
falls E = 0
(x) =
p
p
x
x
A cos
2m|E| ~ + B sin
2m|E| ~
falls E > 0,
(3.1.3)
(3.1.4)
wobei A und B Konstanten sind, die durch die Randbedingungen bestimmt werden. In
jedem der drei Falle verlangt (0) = 0, dass A = 0. Falls E 0, impliziert (a) = 0 dann
B = 0; alle Losungen mit E 0 sind also trivial. Andereseits gibt es eine nicht-triviale
Losung, falls E > 0 von der Form
n2 2 ~2
En =
2ma2
26
(3.1.5)
p
ist, wobei n N, 2m|En | = na ~ . Denn dann ist das Argument des Sinus in (3.1.4) gleich
n f
ur x = a. In diesem Fall haben wir also (a) = 0 selbst wenn B nicht verschwindet.
Wir finden also die nicht-triviale Losungen
nx
.
(3.1.6)
n = Bn sin
a
Die Konstanten Bn werden nun dadurch bestimmt, dass wir verlangen, dass n wie in
(2.5.8) normiert ist. Also wollen wir Bn so finden, dass
Z a
1
2
2 nx
1 = |Bn |
dx sin
= a|Bn |2 .
(3.1.7)
a
2
0
q
Eine Losung ist Bn = a2 . Die normierten Eigenfunktionen sind also
n =
nx
2
sin
a
a
mit
En =
n2 2 ~2
.
2ma2
(3.1.8)
X
En t
(x, t) =
cn n (x)e ~ ,
(3.1.9)
n=0
schreiben, wobei
cn =
a
0
dx n (x) (x, 0) .
(3.1.10)
Diese Formel ist eine direkte Folge der Orthonormalitatsrelation (2.7.1) (falls die Funktionen vollstandig sind). Zum Beispiel konnen wir jetzt die folgende Frage beantworten:
27
Problem: Angenommen das System wird zur Zeit t = 0 durch die konstante Wellenfunktion beschrieben, (x, 0) = 1a . Was ist die Wahrscheinlichkeit, dass bei einer Energiemessung En gemessen wird?
Wie wir zuvor erklart haben ist die Wahrscheinlichkeit, dass die Energie En gemessen
wird einfach |cn |2 in der Entwicklung (3.1.9). Wegen (3.1.10) berechnen wir nun
r Z a
1
2
nx
dx sin
cn =
a
a
a
h0
i
a
a
2
nx
=
cos
n
a 0
a
2
[1 (1)n ] .
(3.1.11)
=
n
Die Wahrscheinlichkeit, die Energie En zu messen ist also Null falls n gerade ist; f
ur
ungerade n gilt andererseits
!2
8
2 2
= 2 2.
(3.1.12)
Pn =
n
n
Als Nebenprodukt dieser Rechnung finden wir auch
X
n>0 odd
8
n2 2
= 1,
(3.1.13)
2
1
=
.
(2m 1)2
8
m=1
(3.1.14)
Ausgedr
uckt durch die Riemannsche -Funktion, deren Reihenentwicklung
(s) =
ns ,
(3.1.15)
n=1
f
ur s > 1 konvergiert, ist die linke Seite
X
X
1
1
2
=
n
(2n)2 = (2) (2) .
2
(2m 1)
4
m=1
n=1
n=1
(3.1.16)
2
4 2
=
.
3 8
6
28
(3.1.17)
3.2
Teilchen im Topf
Als nachstes wollen wir das Teilchen im Topf diskutieren, d.h. den Fall, f
ur den das
Potential durch
V falls a x a,
(3.2.1)
V (x) =
0 sonst
(x) = E(x) ,
(3.2.2)
(3.2.3)
2m
wobei
E=
~2 2
2mE
.
(3.2.4)
2m
~
Damit f
ur x abfallt muss reell gewahlt werden, d.h. E muss negativ sein. Im
folgenden werden wir das Vorzeichen der Wurzel immer so wahlen, dass positiv ist. Mit
dieser Konvention ist die Losung mit den richtigen asymptotischen Randbedingungen
(x) = Aex
if x < a,
(x) = Dex if x > a.
(3.2.5)
In der mittleren Region, a < x < a, ist die zeit-unabhangige Schrodinger Gleichung
~2
2m
(3.2.6)
p
2m(V + E)
~
(3.2.7)
(3.2.8)
Bevor wir mit der detaillierten Diskussion der daraus zusammengesetzten Losung beginnen, ist es n
utzlich zu beobachten, dass falls (x) eine gerade Funktion ist, (x) = (x),
dass dann das auch f
ur H(x) gilt, und entsprechend f
ur ungerade Funktionen. Wir
29
falls x < a
Aex
2B cos(kx) falls a < x < a
(x) =
(3.2.9)
x
Ae
falls x > a.
Als nachstes wollen wir beweisen, dass sowohl als auch bei x = a stetig sein
m
ussen.1 Die Analyse ist im wesentlichen f
ur die beiden Falle x = a gleich; wir werden
daher nur die Klebebedingung bei x = a analysieren. Zunachst beobachten wir
Z a+
~2
~2
dx k (x)
(3.2.10)
(k (a + ) k (a )) =
2m
2m a
Z a
Z a+
=
dx (E + V ) k (x) +
dx E k (x) .
a
Da E und E +V endliche Konstanten sind und da in der Nahe von x = a beschrankt ist,
folgt, dass die rechte Seite der zweiten Zeile im Limeas 0 gegen Null strebt. Daraus
folgt, dass bei x = a stetig sein muss. Insbesondere muss daher auch stetig sein.
Da (x) gerade ist, sind die beiden Klebebedingungen bei x = a gleich, und es
gen
ugt die Bedingung bei x = a zu betrachten. Die relevanten Gleichungen sind dann
Aea = 2B cos(ka) ,
and
(3.2.11)
(3.2.12)
2mV
~2
(3.2.13)
(3.2.15)
Im Fall des Teilchens in der Box (Kapitel 3.1) musste nur stetig sein, nicht aber . Der Umstand,
dass dort nicht stetig sein musst kann man nun auch verstehen!
30
nur eine Losung existiert siehe Abbildung 6. Entprechend zeigt man, dass die Anzahl
der geraden gebundenen Losungen gerade N ist, falls
(N 1)2 2 <
2ma2 V
~2
N 22 .
(b)
(3.2.16)
(c)
(a)
3
2
1
1
1
/2
3 /2
Abbildung 6: Die Losungen in den zwei Fallen V = 4E0 und V = 64E0 . (a) Graphische
Losung der transzendentalen Gleichung mit tan (durchgezogene Linie) und mit cot
(gestrichelte Linie). (b) Losung im Fall des flachen Potentials V = 4E0 . (c) Losungen f
ur
2
2
das tiefe Potential V = 64E0 . [Hier ist a = w/2 und E0 = ~ /(8ma ).]
Die Analyse f
ur den ungeraden Fall ist ahnlich. In diesem Fall haben wir A = D und
B = C, und der Ansatz ist
falls x < a
Aex
2B sin(kx) falls a < x < a
(x) =
(3.2.17)
x
Ae
falls x > a.
and
(3.2.18)
Wiederum teilen wir diese beiden Gleichungen durcheinander, und finden dann
= k cot(ka) .
(3.2.19)
= cot() ,
2ma2 V
~2
(3.2.20)
Diese Gleichungen konnen wiederum numerisch (oder graphisch) gelost werden. Insbesonder kann man leicht sehen, dass es keine ungerade gebundene Losung gibt, falls
2ma2 V
~2
31
2
.
4
(3.2.21)
3.3
Das Stufenpotential
(3.3.1)
wobei V0 > 0. Wie zuvor wollen wir die Losungen der zeit-unabhangigen Schrodinger
Gleichung
H = E
(3.3.2)
bestimmen. Wir sind weiterhin an den Losungen, die die folgende physikalische Konfiguration beschreiben: ein Teilchen kommt von links und wir an der Stufe gestreut. Wir
wollen die Wahrscheinlichkeit bestimmen, mit der das Teilchen reflektiert bzw. transmittiert wird.
F
ur x < 0 ist die zeit-unabhangige Schrodinger Gleichung
d2
= E ,
2m dx2
~2
(3.3.3)
und daher
2mE
d2
.
dx2
~2
Die allgemeinste Losung dieser Gleichung ist
kx
kx
(3.3.4)
k (x) = Ae ~ + Be ~ ,
(3.3.5)
+ V0 = E
(3.3.6)
2m dx2
ist, und daher
d2
2m
= 2 (E V0 ) .
(3.3.7)
2
dx
~
Die Eigenschaften der Losungen zu dieser Gleichung hangen nun stark von dem Vorzeichen
von E V0 ab, namlich ob die Energie der einfallenden Welle kleiner oder grosser als die
Hohe der Potentialstufe ist. Wir diskutieren die beiden Falle separat:
32
lx
(3.3.8)
k (x) = Ce ~ + De ~ ,
p
ur die Konfiguration, in der wir
wobei C und D Konstanten sind und l = 2m(E V0 ). F
interessiert sind, gibt es kein Teilchen das von rechts aus dem Unendlichen eingeschossen
wird; daher interessieren wir uns f
ur die Losung mit D = 0.
Mit denselben Argumenten wie oben sieht man leicht, dass die richtige Klebebedingung
ist, dass k und k bei x = 0 stetig sind. Die Stetigkeit von k bei x = 0 impliziert, dass
A+B =C,
(3.3.9)
(A B)k = Cl .
(3.3.10)
l
k
l
k
und daher
2
C
=
A
1+
1
B
=
A
1+
l
k
l
k
l
k
(3.3.11)
k (x) = C e ~ + D e ~ ,
33
(3.3.12)
p
wobei eine positive reelle Zahl ist = 2m(V0 E). Damit (3.3.12) f
ur x + nicht
(3.3.13)
und daher
2
C
=
A
1 + i
k
1
B
=
A
1+
i
k
i
k
Insbesondere ist also C 6= 0. Die Wahrscheinlichkeit, das Teilchen bei x > 0 anzutreffen
ist also nicht gleich Null! Die Wahrscheinlichkeit nimmt jedoch mit x > 0 exponentiell
ab, und das Teilchen erreicht deshalb niemals x = . (Dies passt dann damit zusammen,
dass in diesem Fall die Reflektions-Wahrscheinlichkeit R = |B/A|2 = 1 eins ist.)
Tunnelstrom
Spannung entlang
x Piezo
IT
lx (Vx )
z
x
Piezokristalle
IT
y
metallische
Spitze
Metall
Abbildung 7: Tunnelmikroskop (Binnig & Rohrer, IBM): Der Tunnelstrom IT fliesst von
der Metallspitze des Mikroskops zur untersuchten Oberflache. Der Strom IT hangt empfindlich (exponentiell) von der Distanz Spitze-Oberflache ab und kann damit als Steuergrosse genutzt werden: Der Strom IT wird in eine Spannung Vz (x, y) verwandelt die den
Piezokristall entlang z manipuliert, so dass IT konstant bleibt. Die Kristalle entlang der
x, y Ebene dienen der Verschiebung der Spitze (scanning). Das Potentialrelief Vz (x, y)
dient dann als Hohenkarte der Oberflache.
Der Umstand, dass das Teilchen eine Region erreichen kann, die klassisch nicht erlaubt
ist, wird manchmal Tunneln genannt. [Streng genommen spricht man nur von Tunneln
im Fall einer endlichen Stufe, d.h. falls V (x) = V0 f
ur 0 x a mit V (x) = 0 f
ur x > a; in
diesem Fall ist die Transmissionswahrscheinlichkeit endlich, obgleich klassisch die Region
34
x > a nicht erreichbar ist siehe spater.] Dieser Tunneleffekt hat wichtige Anwendungen; zum Beispiel ist er der physikalische Mechanismus des Rastertunnelmikroskop (see
Abbildung 7).
3.4
Die obigen Wellenfunktionen sind streng genommen nicht normierbar. Zum Beispiel ist
f
ur den Fall eines freien Teilchens mit festem k (bzw. fester Energie E) die Wellenfunktion k (x) = eikx offenbar nicht quadrat-integrabel, da |k (x)|2 = 1 f
ur alle x. Dieses
Problem tritt immer dann auf, wenn die Menge der Losungen der zeit-unabhangigen
Schrodingergleichung durch einen kontinuierlichen Parameter (wie zum Beispiel k hier)
parametrisiert wird.
Es ist dennoch n
utzlich, diese Losungen auf diese Weise zu beschreiben. Obgleich diese
Losungen nicht quadrat-integrabel sind, kann man dennoch die allgemeinste Losung der
zeit-abhangigen Schrodinger Gleichung durch sie ausdr
ucken. In der Tat gilt weiterhin,
dass die allgemeinste Losung von der Form
Z
Ek t
1
(x, t) =
dk f (k) exk e ~
(3.4.1)
2 R
ist, wobei f (k) die Superposition dieser fundamentalen Losungen beschreibt. Diese Funktion ist dann quadrat-integrabel, falls wir f (k) geeignet wahlen. Eine Superposition dieser
freien Wellen wird oft als Wellenpaket bezeichnet.
Um zu sehen, unter welchen Bedingungen dieses Wellenpaket quadrat-integrabel ist,
erinnern wir uns daran, dass
Z
dx eikx = (k) ,
(3.4.2)
R
definiert. Formal kann man sich die -Distribution als Funktion vorstellen, die u
berall
(mit Ausnahme des Ursprungs) verschwindet. Das obigen Integral kann daher nur von
f (0) abhangen. Der Wert der -Funktion bei 0 ist gerade so, dass das obige Integral
nicht verschwindet sondern gerade f (0) ergibt. [(0) ist daher Unendlich, aber in einem
sehr prazisen Sinn!]
Distributionen kann man in vieler Hinsicht wie Funktionen behandeln; zum Beispiel
kann man ihre Ableitung dadurch definieren, dass die u
bliche partielle Integrationsformel
weiterhin gilt. Genauer gesagt ist (x) die Distribution, die durch
Z
Z
x=
= f (0)
35
(3.4.4)
definiert ist. In mancher Hinsicht verhalten sich daher Distributionen wie normale Funktionen; der Hauptunterschied besteht darin, dass das Produkt zweier Distributionen im
allgemeinen nicht wohl-definiert ist.
Mit diesen Vorbereitungen konnen wir nun die Normalisierung der Funktion (3.4.1)
analysieren. Dazu berechnen wir
Z
Z
Z
Z
E t
Ek t
1
k
2
dx
dk
dk f (k )f (k)k (x)k (x)e ~ e ~
dx |(x, t)| =
2
R
R
Z RZ R
E t
Ek t
k
=
dk
dk f (k )f (k)e ~ e ~ (k k)
ZR
ZR
Ek+l t
Ek t
=
dk
dl f (k + l)f (k)e ~ e ~ (l)
R
ZR
=
dk |f (k)|2 ,
(3.4.5)
R
haben; in der vorletzten Zeile haben wir ausserdem substituiert k k = l. Das Wellenpaket
(3.4.1) ist also genau dann korrekt normiert, falls die Funktion f quadrat-integrierbar (mit
Integral 1) ist.
Als Beispiel betrachten wir das Wellenpaket, das durch
a
2
2
(3.4.6)
f (k) = 1/4 ea (kk0 ) /2
definiert wird; dies beschreibt eine um k0 zentrierte Impulsverteilung deren Breite umgekehrt proportional zu a ist. In diesem Fall kann man das Wellenpaket direkt berechnen:
Z
2
~k
(x, t) =
dk eikx ei 2m t f (k)
Z
k2
a
ikx i ~2m
t a2 (kk0 )2 /2
=
dk
e
e
e
1/4
2
~t
a a2 k02
2
k 2 a2 +i 2m
ek(ix+a k0 )
=
e 2
dk e
1/4
Z
c 2
a a2 k02 c22
(bk 2b
) ,
2 e 4b
=
(3.4.7)
e
dk
e
1/4
wobei wir
~t
a2
+i
,
c = ix + a2 k0
(3.4.8)
b2 =
2
2m
definiert haben. Das Gausssche Integral kann nun leicht durchgef
uhrt werden, und wir
finden
Die Funktion ist also wieder von Gaussschem Typ, d.h. sie ist proportional zu ea x ,
wobei nun die Breite proportional zu 2b a ist. Da |b| mit der Zeit anwachst, |b| t,
lauft das Wellenpaket langsam auseinander.
36
3.5
Das Delta-Funktionspotential
~2
2m
(3.5.1)
F
ur x < 0 ist die allgemeinste Losung
(x) = Ae
kx
~
+ Be
kx
~
(3.5.2)
k
wobei, wie zuvor, E = 2m
. F
ur x > 0 ist die Losung im wesentlichen dieselbe; da wir uns
aber nur f
ur die Losungen interessieren, die ein von links einfliegendes Teilchen beschreiben, tragt nur einer der beiden Terme bei, namlich
(x) = Ce
kx
~
(3.5.3)
( () ()) + V (0) = E
(x) .
(3.5.4)
2m
Nun bilden wir den Limes 0. Da (x) auf dem Intervall [, ] beschrankt ist, verschwindet die rechte Seite und wir erhalten
(0+ ) (0 ) =
2mV
~2
(0) ,
(3.5.5)
x0
x0
(3.5.6)
x<0
x>0
ben
utzt haben. Damit die beiden Grenzwerte (0+ ) und (0 ) existieren, muss bei
x = 0 stetig sein; in der obigen Notation bedeutet dies
(0+ ) (0 ) = 0 .
(3.5.7)
F
ur den obigen Ansatz (3.5.7) muss also gelten
A+B =C,
37
(3.5.8)
[C (A B)] = 2
k
~
2mV (A + B)
B=
~2
(3.5.9)
k ~
mV
1
1
(3.5.10)
k ~
mV
k ~
mV
(3.5.11)
1
1
+
2
2
mV
mV
m V
und
2
C
T = =
A
2 2
k ~
~ k
mV
m2 V 2
=
2 2 .
~
1 k
1
1 + m~ 2kV 2
mV
k ~
mV
k ~
mV
(3.5.12)
(3.5.13)
R+T =1
(3.5.14)
gilt. Wie erwartet wachst die Transmissionswahrscheinlichkeit T mit der Energie, wohingegen die Reflektionswahrscheinlichkeit in gleichem Mass mit wachsender Energie abnimmt.
3.6
Resonanzen
Schliesslich wollen wir nochmals das Potential aus Kapitel 3.2 analysieren, aber dieses Mal
f
ur den Fall positiver Energie. In der klassischen Physik w
urde man erwarten, dass ein
Teilchen mit positiver Energie nicht von dem Potential beeinflusst wird. Wie wir jedoch
nun zeigen wollen, ist diese Intuition in der Quantenmechanik nicht richtig.
Wir machen nun den Ansatz, dass f
ur x < a die Wellenfunktion durch
2mE
x < a :
(x) = ekx + rekx ,
wobei
k=
> 0.
(3.6.1)
~
gegeben ist. Diese Wellenfunktion beschreibt ein von links einfliegendes Teilchen der Energie E, sowie das entsprechend reflektierte Teilchen. (Die Reflektionswahrscheinlichkeit ist
also R = |r|2 .) F
ur a < x < a ist die Losung
p
2m(E + V )
lx
lx
> 0 . (3.6.2)
a < x < a :
(x) = Ae + Be
,
wobei
l=
~
38
F
ur x > a haben wir schliesslich
x>a:
kx
(x) = te
kx
+ De
wobei
k=
2mE
~
>0.
(3.6.3)
Wie zuvor wahlen wir D = 0, da kein Teilchen von rechts einfliegen soll. Die Transmissionswahrscheinlichkeit ist dann T = |t|2 .
Wie zuvor konnen wir dieses Problem dadurch losen, dass wie die verschiedenen
Losungen gemass der richtigen Klebebedingungen zusammenf
ugen. Bei x = a haben wir
die beiden Bedingungen
t = Ae(lk)a + Be(l+k)a
(k/l) t = Ae(lk)a Be(l+k)a ,
die sich leicht nach A und B auflosen lassen
k
1 (lk)a
e
t 1+
A =
2
l
1 (l+k)a
k
B =
.
e
t 1
2
l
Die Klebebedingungen bei x = a sind nun
eka + reka = Aela + Bela
eka reka = (l/k) Aela Bela .
e2ika
.
2 +k 2
cos(2la) i l 2lk
sin(2la)
(3.6.4)
Schliesslich dr
ucken wir l und k durch E und V aus, und erhalten die Transmissionswahrscheinlichkeit
1
(3.6.5)
T = |t|2 =
2 (2la)V 2 .
1 + sin
4E(E+V )
Wir bemerken, dass im allgemeinen T < 1. Dies bedeutet daher, dass in der Quantenmechanik selbst Teilchen mit E > 0 von dem Potential beeinflusst werden!
Im Rahmen dieser Analyse konnen wir nun auch den Tunneleffekt (den wir schon zuvor
angesprochen haben siehe Kapitel 3.3) im Detail verstehen. Der Tunneleffekt tritt dann
auf, falls V negativ ist (so dass V positiv ist) und 0 < E < V . In der obigen Analyse
ging nicht ein, dass V positiv war V spielt nur f
ur die Analyse der Wellenfunktion f
ur
a < x < a eine Rolle, f
ur die die Frage, ob der Exponent reell oder imaginar ist, irrelevant
ist. Wir konnen also einfach in dem Ausdruck f
ur die Transmissionswahrscheinlichkeit
39
~2 2
8ma2
n2 V = E0 (n)2 V .
(3.6.7)
ik
l
.
l
ik
(3.6.8)
ik
,
l
oder
tan(la/2) =
ik
.
l
(3.6.9)
En
l d l
4 k
(3.6.12)
2a .
=
+
wobei
l
k
dE En
40
Als Funktion von E verhalt sich also t(E) in der Nahe einer Resonanz wie
e2ika t(E)
=
i/2
E (En i/2)
(3.6.13)
und das Absolutquadrat T = |t(E)|2 wird durch ein Lorentzkurve beschrieben (siehe
Abbildung 8(a))
2 /4
T = |t|2
.
(3.6.14)
=
(E En )2 + 2 /4
F
ur das weitere ist es bequem, die Amplitude t(E) als
t(E) = |t| ei(E)
(3.6.15)
zu schreiben, wobei (E) eine Phase ist. Man sieht leicht, dass diese Phase jedes Mal,
wenn eine Resonanz durchlaufen wird, um wachst siehe Abbildung 8(b).
|t |2
(a)
E n2
E n1
r (E)
(b)
E n1
E n2
41
Wir wollen nun versuchen, die Axiome der Quantenmechanik ein wenig abstrakter zu
formulieren. Bisher haben wir die verschiedenen Systeme, die wir besprochen haben,
durch komplexe Wellenfunktionen (x, t) beschrieben, deren Zeitentwicklung durch die
Schrodinger Gleichung gegeben ist. Da die Schrodinger Gleichung eine lineare Differentialgleichung ist, gilt das Superpositionsprinzip, d.h. die Summe zweier Losungen und
das Vielfache einer Losung sind immer auch Losungen. Der Raum der Losungen ist daher ein komplexer Vektorraum. (Im Fall der Wellenfunktionen ist der Losungsraum ein
Unterraum des Raums der komplexwertigen Funktionen auf dem Konfigurationsraum.)
Damit die Wellenfunktion eine vern
unftige Wahrscheinlichkeitsinterpretation besitzt,
war weiterhin wichtig, dass sie normierbar war, d.h. dass
Z
d~x |(~x, t)|2 = 1
f
ur jedes t.
Rf
Formal gesehen bedeutet dies, dass der komplexe Vektorraum der Losungen ein inneres
Produkt besitzen muss. Einen solchen Vektorraum nennt man einen Hilbertraum.
4.1
Der Hilbertraum
Sei H ein komplexer Vektorraum. H ist ein Hilbertraum, falls H eine positiv-definite
sesqui-lineare Form (Skalarprodukt) besitzt, bez
uglich deren Norm H vollstandig ist. Das
Skalarprodukt h|i ist eine Abbildung
h|i : H H C ,
(4.1.1)
(4.1.2)
(4.1.3)
und
42
h|i = 0 = 0 .
(4.1.4)
(4.1.5)
Wegen (iii) ist der Ausdruck unter der Wurzel immer positiv, und von daher ist die Norm
ebenso positiv. Weiterhin zeigt man leicht mit Hilfe der Cauchy-Schwarz Ungleichung
|h|i|2 h|i h|i ,
(4.1.6)
(4.1.7)
Schliesslich ist ein komplexer Vektorraum H mit einem Skalarprodukt ein Hilbertraum,
falls H bez
uglich der Norm (4.1.5) vollstandig ist: dies bedeutet, dass jede Cauchy Folge
von Vektoren n H zu einem Element in H konvergiert
lim n = H .
[Eine Cauchy Folge n ist eine Folge, die die folgende Eigenschaft besitzt: f
ur jedes > 0
gibt es ein N N, so dass f
ur alle n, m > N,
||n m || < .]
Das einfachste Beispiel eines Hilbertraumes ist der n-dimensionale komplexe Vektorraum Cn , f
ur den das Skalarprodukt durch
* x1
y1 + X
n
.. ..
xi yi
(4.1.8)
. . =
i=1
yn
xn
4.2
Der L2 Raum
Der Hilbertraum, der der Wellenmechanik zu Grunde liegt, ist der Raum der quadratintegrierbaren Funktionen, der sogenannte L2 Raum. Im einfachsten Fall eines 1-dimensionalen Konfigurationsraums (die Verallgemeinerung f
ur f -dimensionale Systeme ist relativ
2
2
offensichtlich) ist der relevante L Raum L (R), der Vektorraum der komplex-wertigen
Funktionen f : R C, f
ur die gilt
Z
dx|f (x)|2 < .
(4.2.1)
R
43
(4.2.2)
Die L2 Bedingung (4.2.1) ist also einfach die Bedingung, dass hf |f i < .
Auf der Menge dieser Funktionen definieren wir die Vektorraumoperationen
(f + g)(x) = f (x) + g(x)
(f )(x) = f (x) .
Es ist klar, dass das obige Skalarprodukt bez
uglich dieser Operationen sesquilinear ist,
d.h. (i) und (ii) erf
ullt sind. Weiterhin ist offensichtlich, dass hf |f i 0. Um zu sehen,
dass das Skalarprodukt tatsachlich positiv definit ist (Eigenschaft (iii)), m
ussen wir genauer erklaren, wie das Integral definiert ist. Das richtige Integral-Mass ist das sogenannte
Lebesgue-Mass, f
ur das auch nicht-stetige Funktionen integrierbar sind. Der Hilbertraum
ist dann der Raum aller lebesgue-integrierbaren Funktionen, wobei wir Funktionen identifizieren, die sich nur auf einer Menge vom Mass Null unterscheiden; das Skalarprodukt ist
dann auf diesem Quotientenraum eindeutig definiert und ist ferner automatisch positivdefinit.
Mit Hilfe der Cauchy-Schwarz Ungleichung kann man nun leicht beweisen, dass die
obigen Vektorraumoperationen tatsachlich auf L2 schliessen. F
ur den Fall der Skalarmultiplikation ist das relativ offensichtlich, denn es gilt wegen der Sesquilinearitat
hf |f i = ||2 hf |f i .
(4.2.3)
(4.2.4)
Die letzten beiden Terme kann man nun mit Hilfe der Cauchy-Schwarz Ungleichung
abschatzen, und es folgt daher, dass hf + g|f + gi < . Dies beweist, dass die Vektorraumoperationen auf L2 schliessen.
Mit demselben Argument folgt u
brigens auch, hf |gi < , falls f und g in L2 liegen
(und daher (4.2.1) erf
ullen); das Skalarprodukt ist also f
ur alle Elemente von L2 wohldefiniert.
Die verbleibende Eigenschaft, namlich dass der so-definierte Vektorraum vollstandig
ist, erfordert ein wenig mehr Aufwand. Die Details spielen aber f
ur das Weitere keine
wichtige Rolle, und wir werden daher nicht darauf eingehen.
4.3
Separable Hilbertr
aume
Der L2 (R) Hilbertraum ist unendlich-dimensional und separabel. Dies bedeutet, dass er
eine abzahlbar unendliche Basis fn besitzt, d.h. linear unabhangige Vektoren fn , wobei
44
X
f=
cn fn ,
cn C .
(4.3.1)
n=1
Wann immer ein Hilbertraum eine abzahlbare Basis besitzt, kann man daraus eine orthonormale Basis konstruieren. Eine orthonormale Basis hat die Eigenschaft, dass
hfn |fm i = m,n .
(4.3.2)
Allgemeiner nennen wir zwei Vektoren f und g eines Hilbertraumes orthogonal, falls
hf |gi = 0.
Die orthonormale Basis kann mit Hilfe des Schmidtschen Orthogonalit
atsverfahren
rekursiv konstruiert werden. Dabei setzen wir
h1 = f1 /||f1 ||
g2 = f2 hh1 |f2 ih1
g3 = f3 hh1 |f3 ih1 hh2 |f3 ih2
h2 = g2 /||g2 ||
h3 = g3 /||g3 ||
usw.
4.4
Wie wir gesehen haben sind die Wellenfunktionen Elemente eines Hilbertraumes, namlich
des L2 Raumes. Wie wir weiter gesehen haben entsprechen physikalischen Observablen (so
wie zum Beispiel dem Impuls oder der Energie, d.h. der Hamiltonfunktion) DifferentialOperatoren, die auf diesen Wellenfunktionen wirken. Zum Beispiel ist der Impuls p der
Operator
f (x)
.
(4.4.1)
p : L2 (R) L2 (R)
f (x) 7 p f (x) = i~
x
[Wir ignorieren hier das Problem, dass p nicht auf dem ganzen L2 Raum definiert ist.]
Als Operator auf L2 ist p ein linearer Operator, d.h. f
ur , C und f, g H gilt
p(f + g) = p(f ) + p(g) .
(4.4.2)
Ein anderes (noch einfacheres) Beispiel ist der Ortsoperator; er ist einfach durch Multiplikation mit x definiert
x : L2 (R) L2 (R)
(4.4.3)
45
(4.4.5)
(4.4.7)
f
ur alle f, g H
(4.4.8)
eindeutig bestimmen. Eine hinreichende Bedingung, die (4.4.7) garaniertiert, ist also einfach, dass
A = A ,
(4.4.9)
d.h. dass der Operator selbst-adjungiert ist. Man pr
uft leicht nach (im wesentlichen mit
derselben Rechnung wie in Kapitel 2.4), dass der Impuls und Hamiltonoperator selbstadjungiert sind. Wir verlangen daher im allgemeinen: physikalische Observable werden
durch selbst-adjungierte Operatoren auf dem Hilbertraum dargestellt.
Im endlich dimensionalen Fall, d.h. falls H = Cn , konnen wir jeden linearen Operator
A durch eine Matrix M(A) beschreiben. (Diese Matrix hangt nat
urlich von der Wahl
der Basis ab; im folgenden wahlen wir immer eine Orthonormalbasis, z.B. die Basis e1 =
(1, 0, . . . , 0), e2 = (0, 1, 0, . . . , 0), usw.) Dann ist die Matrix des adjungierten Operators
A gerade die hermitesch konjugierte Matrix von A, d.h.
M(A ) = M(A)t .
(4.4.10)
Die Observablen entsprechen also gerade den hermiteschen Matrizen. Der allgemeine Fall
ist die unendlich-dimensionale Verallgemeinerung dieser Bedingung.
46
4.5
Als nachstes wollen wir erklaren, was die moglichen Resultate einer Messung der Observablen A sind und was die richtige Wahrscheinlichkeitsinterpretation ist, d.h. wir wollen
die Axiome (iv) und (v) aus Kapitel 2.7 auf allgemeine Observable verallgemeinern.
Der Einfachheit halber nehmen wir zunachst an, dass der Hilbertraum endlich-dimensional (von Dimension N) ist. Jeder selbst-adjungierte Operator (d.h. jede hermitesche
Matrix) kann dann diagonalisiert werden; dies bedeutet, dass wir eine Basis von H aus
Eigenvektoren von A finden konnen. Wir bezeichnen diese Eigenvektoren mit n , wobei
n = 1, 2, . . . N, und die zugehorigen Eigenwerte mit n , n = 1, 2, . . . , N,
A n = n n .
(4.5.1)
Ohne Beschrankung der Allgemeinheit konnen wir annehmen, dass die n normiert sind,
hn |n i = 1. Weiterhin zeigt man leicht (das ist die analoge Rechnung zu jener in Kapitel
2.7), dass die Eigenvektoren zu verschiedenen Eigenwerten zueinander orthogonal sind,
hm |n i = 0 falls m 6= n . In der Tat folgt dies aus
n hm |n i = hm |An i = hAm |n i
= m hm |n i .
(4.5.2)
Wenden wir diese Gleichung mit m = n an, folgt dass alle n reell sein m
ussen. Dann
konnen wir die beiden Seiten der Gleichung voneinander abziehen und deduzieren, dass
hm |n i = 0 falls m 6= n . Falls es mehrere Eigenvektoren zu demselben Eigenwert
gibt, konnen wir diese separat (mit dem Schmidt-schen Orthogonalisierungsverfahren)
orthonormieren; wir konnen daher also immer annehmen, dass die Eigenvektoren eines
selbst-adjungierten Operators (im endlich-dimensionale Fall) orthonormiert sind,
hn |m i = m,n .
(4.5.3)
F
ur das folgende ist es bequem die sogenannte Dirac-Notation einzuf
uhren. Dabei
bezeichnen wir Vektoren des Hilbertraums als kets, namlich als
|ai H .
(4.5.4)
(4.5.5)
(4.5.6)
Mit Hilfe dieser Notation kann man nun leicht Spektralprojektoren definieren. Dazu
betrachten wir
Pn = |n ihn | : H H .
(4.5.7)
Man rechnet leicht nach (wegen der Normierung von n ), dass Pn ein Projektor ist,
Pn2 = Pn . Ausserdem ist offensichtlich, dass Pn auf den Eigenraum |n i projeziert. Es ist
dann klar, dass
X
X
A=
n Pn =
n |n ihn | .
(4.5.8)
n
(iv) die moglichen Ergebnisse einer Messung von A sind die verschiedenen
Eigenwerte n .
(v) die Wahrscheinlichkeit, dass als Messergebnis auftritt ist einfach
X
|hm |i|2 .
W () =
(4.5.10)
m =
(t) = H(t)
t
(4.5.11)
beschrieben. [Das ist nun eine Differentialgleichung in H. Man sieht leicht, dass diese
4.6
(4.6.2)
(x) = (p) = R ,
(H) = R+ .
(4.6.3)
F
ur das Teilchen im Potentialtopf (siehe Kapitel 3.2) gilt andererseits
Potentialtopf:
(H) = {E1 , . . . , EN } R+ ,
(4.6.4)
(H) = {~( 12 + n) : n = 0, 1, 2, . . .} .
(4.6.5)
Die Verallgemeinerung der Wahrscheinlichkeitsinterpretation auf beliebige Observablen beruht auf dem Spektralsatz. Falls das Spektrum rein diskret ist (wie zum Beispiel
im endlich-dimensionalen Fall), dann gilt
X
f (A) :=
f (a)Pa ,
(4.6.6)
a(A)
49
wobei f : R C eine Funktion und f (A) : H H ein Operator ist. Die Zuordnung
f 7 f (A) hat die Eigenschaften
(1 f1 + 2 f2 )(A) = 1 f1 (A) + 2 f2 (A) ,
(f1 f2 )(A) = f1 (A)f2 (A) ,
f(A)
= f (A) ,
1I f
ur f (x) 1 ,
f (A) =
A f
ur f (x) = x ,
(1 , 2 C) ,
(4.6.7)
(4.6.8)
(4.6.9)
(4.6.10)
sowie Stetigkeit bzgl. f , auf die wir nicht naher eingehen. Der Spektralsatz besagt, dass
es auch im unendlich-dimensionalen Fall eine eindeutige Zuordnung gibt
f 7 f (A) ,
(4.6.11)
PI (x)
I
(I1 I2 = ) ,
Wir interpretieren dann W (I) als die Wahrscheinlichkeit, dass A im Zustand |i einen
Messwert a I annimt. Die obigen drei Bedingungen garantieren, dass das eine vern
unftige Definition ist.
Weiterhin definieren wir dann den Erwartungswert von A im Zustand durch
Z
hAi = dW ((, ]) = h|A|i ,
(4.6.12)
da
4.7
In der urspr
unglichen Beschreibung der Quantenmechanik als Wellenmechanik war unser
Hilbertraum der Ortsraum L2 (R), wobei wir R (bzw. R3 ) als den physikalischen Raum
interpretierten. In dieser Betrachtung wirkt der Ortsoperator einfach als Multiplikationsoperator
x : L2 (R) L2 (R)
(f )(y) 7 [x(f )](y) = yf (y) ,
(4.7.1)
wohingegen der Impulsoperator der Differentialoperator
p : L2 (R) L2 (R)
f (x)
x
(4.7.2)
ist. Aus der Perspektive der obigen Analyse ist nun klar, dass dies nur eine spezielle
Darstellung der Quantenmechanik ist. Zum Beispiel konnten wir genauso gut im Impulsraum arbeiten, der (f
ur den Fall eines freien Teilchens) auch wieder gerade L2 (R)
(bzw. L2 (R3 )) ist. In dieser Impulsdarstellung ist dann der Impulsoperator einfach ein
Multiplikationsoperator, namlich
p : L2 (R) L2 (R)
(4.7.3)
f (p)
.
p
(4.7.4)
Die Definition des Ortsoperators im Impulsraum ist dadurch vorgeschrieben, dass er dieselbe algebraische Struktur wie in der Ortsdarstellung besitzen muss. Zum Beispiel muss
der Kommutator von x und p in beiden Beschreibungen u
bereinstimmen. Der Kommutator zweier Operatoren A und B ist einfach
[A, B] = A B B A .
(4.7.5)
Um diesen Kommutator auszurechnen, wenden wir die rechte Seite auf ein beliebiges
Element des Hilbertraums an. Zum Beispiel gilt im Ortsraum
f (x)
f (x)
= ~x
(4.7.6)
x p f (x) = x ~
x
x
und
f (x)
(4.7.7)
(4.7.8)
(4.7.9)
52
Viele Quantenphanomene sind schlussendlich eine Folge davon, dass die Observablen in
der Quantenmechanik durch selbst-adjungierte Operatoren beschrieben werden, die im
allgemeinen nicht miteinander vertauschen. Zum Beispiel haben wir ja bereits im letzten
Kapitel gesehen, dass der Orts- und Impulsopertor nicht miteinander vertauschen, sondern
dass gilt
[x, p] = ~ .
Wie wir in K
urze zeigen wollen, ist der Umstand, dass dieser Kommutator nicht verschwindet der wesentliche Grund f
ur die Heisenbergsche Unscharferelation. Bevor wir
diese ableiten, wollen wir zunachst ein besseres Gef
uhl daf
ur bekommen, was die physikalische Konsequenz nicht-verschwindender Kommutatoren ist.
5.1
Nicht-vertauschende Observable
Um einige der Komplikationen zu illustrieren, die auftreten, wenn Observable nicht vertauschen, betrachten wir das folgende toy model. Wir arbeiten in der Energiedarstellung
und nehmen (der Einfachheit halber!) an, dass dieser Hilbertraum nur zwei-dimensional
ist. Wir bezeichnen die zwei orthonormalen Eigenvektoren des Hamiltonoperators H mit
1 und 2 , wobei
H1 = E1 1
H2 = E2 2 .
(5.1.1)
Die beiden Zustande 1 und 2 bezeichnen also zwei Losungen der zeit-unabhangigen
Schrodinger Gleichung. Wir wollen weiterhin annehmen, dass es ausserdem eine Observable S gibt, die auf diesen Zustande als
S1 = 2
S2 = 1
(5.1.2)
(5.1.3)
Beide Matrizen sind in der Tat hermitesch, denn sie entsprechen (als Observable) selbstadjungierten Operatoren. Nun berechnen wir den Kommutator:
E1 0
0 1
0 1
E1 0
[H, S] =
0 E2
1 0
1 0
0 E2
0 E2
0 E1
=
E1 0
E2 0
0 1
.
(5.1.4)
= (E1 E2 )
1 0
Der Kommutator verschwindet also nicht, vorausgesetzt, dass E1 E2 6= 0, d.h. dass die
beiden Energieeigenwerte nicht u
bereinstimmen.
53
Die Eigenvektoren des Hamiltonoperators sind 1 und 2 , aber offensichtlich sind diese
Vektoren nicht Eigenvektoren von S. Dies ist notwendigerweise der Fall, da ja S und H
nicht miteinander vertauschen. Da wir auch die Wahrscheinlichkeiten einer S-Messung
ausrechnen wollen, sollten wir die S-Eigenvektoren, sowie ihre Eigenwerte bestimmen.
Die Eigenwerte von S sind einfach die Losungen der charakteristischen Gleichung
1
det(S 1) = det
= 2 1 .
(5.1.5)
1
Wir sehen also, dass die beiden Eigenwerte von S gerade 1 sind. Weiterhin findet man
leicht, dass der (normierte) Eigenvektor mit Eigenwert s1 = +1 durch
1
1 = (1 + 2 )
2
(5.1.6)
(5.1.7)
Zur Zeit t wollen wir nun wiederum S bestimmen; um die zugehorigen Wahrscheinlichkeiten auszurechnen, ist es nun geschickt, (t) wiederum in der Eigenbasis von S auszurechnen. Um die Koeffizienten in dieser Basis zu berechnen, ben
utzen wir nun den jetzt
wohl-bekannten Trick; wir machen also den Ansatz
(t) = a1 (t)1 + a2 (t)2 ,
(5.1.10)
und ben
utzen, dass 1 und 2 orthonormal sind, um abzuleiten, dass
ai (t) = hi |(t)i
i = 1, 2 .
(5.1.11)
Wir berechnen
E1 t
E2 t
1
h1 + 2 |1 e ~ + 2 e ~ i
2
E t
1 E1 t
~2
~
e
+e
=
2
(E E )t
1 (E1 +E2 )t (E1 E2 )t
1 2~ 2
2
~
2
~
e
+e
=
e
2
(E +E )t
(E1 E2 )t
1 2~ 2
= e
cos
2~
a1 (t) =
(5.1.12)
und
E2 t
E1 t
1
h1 2 |1 e ~ + 2 e ~ i
2
E t
1 E1 t
~2
~
=
e
e
2
(E E )t
1 (E1 +E2 )t (E1 E2 )t
1 2~ 2
2
~
2
~
e
e
e
=
2
(E +E )t
(E1 E2 )t
1 2~ 2
.
= e
sin
2~
a2 (t) =
(5.1.13)
(E1 E2 )t
2~
(5.1.15)
5.2
Die Unsch
arfe einer Observablen
Wie wir zuvor erklart haben, ist der Erwartungswert einer Observablen A im Zustand
einfach
hAi = h|A i .
(5.2.1)
Wir definieren nun die Unsch
arfe der Observable A durch
rD
2 E
A hAi
.
A =
Das Quadrat der Unscharfe ist eine nicht-negative reelle Zahl, denn
D
2 E
(A)2 =
A hAi
D
E
= | A hAi A hAi
D
E
=
A hAi | A hAi 0 ,
(5.2.2)
(5.2.3)
dx x
n
a
a
0
56
Der Randterm ist einfach gleich a, und im zweiten Integral ist nur das Integral von x
von Null verschieden (das cos sin-Integral verschwindet, wie aus der Fouriertheorie wohlbekannt). Da das x-Integral gerade a/2 ist, folgt also insgesamt
hxi =
a
.
2
(5.2.7)
2 a
(~)
dx sin
sin
hpi =
a 0
a
x
a
Z a
2~ n
nx
nx
=
cos
dx sin
a a 0
a
a
~ h 2 nx ia
=
sin
a
a
0
= 0.
(5.2.8)
(5.2.9)
2
= hn | ~2 2 n i
x +
*
~2 n2 2
=
n
n
a2
=
n2 2 ~2
.
a2
(5.2.10)
(p) (x)
6
n2 2 ~2 a2
1 2 2
=
a2 12
n
2 2 2
6
n ~
1 2 2
=
12
n
2
2 2 2
~
n ~
.
=
12
2
(5.2.11)
2 6
= 0.32 ~2
12
57
(5.2.12)
unabhangig von n. Das ist ein Beispiel der so-genannten Heisenbergschen Unscharferelation, namlich dass das Produkt der Unscharfen von x und p von unten beschrankt ist.
Wie wir in K
urze zeigen werden, gilt ganz allgemein
2
(p) (x)
~2
(5.2.13)
5.3
1
|h[A, B]i| .
2
(5.3.1)
(5.3.3)
(5.3.4)
wobei C und D selbst-adjungiert sind. Wie bezeichnen den Zustand des Systems mit .
Ausserdem sei s eine beliebige reelle Zahl. Dann gilt
0 h(C + sD)|(C + sD)i = h|(C + sD) (C + sD)i = h(C + sD) (C + sD)i .
(5.3.5)
58
(5.3.6)
(5.3.7)
(5.3.8)
1 h[C, D]i
.
2 hD2 i
(5.3.9)
.
2 hD2 i
4 hD2 i
4 hD2 i
(5.3.10)
(5.3.11)
(5.3.12)
(5.3.13)
Die obige Analyse zeigt noch nicht, dass die untere Schranke, die wir f
ur das Produkt der
Schwankungsquadrate abgeleitet haben, optimal ist. Mit Hilfe eines Beispiels konnen wir
das aber nun sehen. Dazu betrachten wir das Gausssche Wellenpaket
a2 (xx0 )2
a
(x) = 1/4 e 2
.
(5.3.14)
59
= x0
(5.3.15)
a
=
dx (x x0 )2 ea (xx0 )
Z
1
2
2
dx ea (xx0 )
=
2 a
1
1
=
= 2 .
2 a a
2a
(5.3.16)
Andererseits ist
Z
a2 (xx0 )2
a2 (xx0 )2
a
hpi = i~
e 2
dx e 2
x
3 Z
a
2
2
= i~
dx (x x0 ) ea (xx0 ) = 0 .
(5.3.17)
2
2
2
dx e 2
e
(p) = hp i =
x2
Z
Z
a3 ~2
a5 ~2
2
2
a2 (xx0 )2
=
dx e
dx (x x0 )2 ea (xx0 )
1
1
(5.3.18)
= a2 ~2 a4 ~2 2 = a2 ~2 .
2a
2
F
ur dieses Beispiel ist also
(x)2 (p)2 =
1 1 2 2 ~2
a~ =
.
2a2 2
4
(5.3.19)
In diesem Fall ist also gerade die Heisenbergsche Unscharferelation minimal erf
ullt.
60
Der harmonische Oszillator ist eine Drosophila der Quantenmechanik. Viele Probleme lassen sich durch eine Abbildung auf den (verschobenen) harmonischen Oszillator
zur
uckf
uhren und werden damit exakt losbar. Eine ausf
uhrliche Diskussion drangt sich
deshalb auf.
6.1
Die Lo
sung
(6.1.1)
m
(6.1.3)
q,
p = i m~x .
x =
~
x2 + x2
(6.1.4)
2E
.
~
(6.1.5)
6.1.1
Die konventionelle L
osung
Der f
uhrende Term f
ur x ist (x2 ),
(x2 x2 ) = 0 .
Daher muss asymptotisch wie
ex
(6.1.6)
2 /2
(6.1.7)
(6.1.9)
n0
= 0,
= 0,
= 0,
= 0.
(6.1.10)
Gemass Voraussetzung ist a0 6= 0, weshalb s = 0 oder s = 1 sein muss; weiter setzen wir
ohne Beschrankung der Allgemeinheit a1 = 0. Aus diesen Bedingungen ergibt sich
f
ur s = 0 ist H(x) = a0 + a2 x2 +
gerade in x,
f
ur s = 1 ist H(x) = x(a0 + a2 x2 + )
ungerade in x.
Der Umstand, dass es Losungen mit definitiver Paritat gibt, ist wiederum eine Folge
davon, dass das Potential symmetrisch is, V (q) = V (q).
Aus den obigen Rekursionsformeln f
ur an folgt, dass an+2 /an 2/n; falls die Folge
an nicht abbricht, geht H exp(x2 ) f
ur x . Diese Funktion u
berkompensiert dax2 /2
x2 /2
her im Ansatz (x) = H(x)e
das Abfallverhalten von e
. Die Randbedingung
limx (x) 0 kann also nur befriedigt werden wenn die Reihe abbricht. Eine detaillierte Analyse zeigt, dass das nur dann der Fall ist, falls n = 2n + 1. Dann erf
ullt H(x)
die Differentialgleichung
0 = H 2xH + 2nH .
(6.1.11)
Diese Differentialgleichung hat eine polynomiale Losung, namlich das n-te Hermite Polynom Hn (x). (Zum Beispiel ist H0 = 1, H1 = 2x, H2 = 4x2 2, . . ..)
62
n (x) = Nn Hn (x) e 2 ,
mit zugehorigen Eigenwerten n = 2n + 1, also
1
En = ~ n +
.
2
(6.1.12)
(6.1.13)
Nn = N0 / 2n n! .
N0 = 1/ 1/4 ,
Die elegante L
osung
Diese Losung basiert auf einer Operatortechnik mit Auf- und Absteige Operatoren. Diese
Operatoren spielen in vielen Bereichen der theoretischen Physik eine wichtige Rolle
zum Beispiel treten sie in der Feldquantisierung freier Felder auf. Wir definieren
r
1
i
1 m
p ,
a (x + x ) =
q+
~
m~
2
2
1
a (x x ) .
(6.1.14)
2
Die Umkehrung ist dann
1
x = (a + a ) ,
2
1
x = (a a ) .
2
(6.1.15)
(x2 + x2 )
2
~
[(a a )(a a ) + (a + a )(a + a )]
=
4
~
=
(a2 + aa + a a a2 + a2 + aa + a a + a2 )
4
1
= ~ aa [a, a ] .
2
H =
63
(6.1.16)
F
ur den Kommutator finden wir (vergleiche (6.1.4))
1
[(x + x )(x x ) (x x )(x + x )]
2
1
=
(xx + x x)2 = 1.
2
[a, a ] =
(6.1.17)
(6.1.18)
(6.1.19)
(6.1.20)
(6.1.21)
Die Zustande
a |ni,
a|ni
(6.1.22)
(6.1.23)
Die Operatoren a und a erhohen und erniedrigen also den Eigenwert n eines Eigenzustandes |ni um 1; diese Eigenschaft ist eine Folge der Kommutationsbeziehungen zwischen
a, a , und N = a a (die letztlich das Spektrum definiert),
[a, a ] = 1,
[N, a ] = a ,
[N, a] = a.
(6.1.24)
(6.1.25)
(6.1.26)
(6.1.28)
d.h. der Operator a annihiliert |0i. Es folgen dann alle |ni durch Anwendung von a ,
|ni =
(a )n
|0i,
n!
n 0.
(6.1.29)
(6.1.30)
somit zahlt n die vorhandenen Energiequanten ~. Die Anwendung von a erzeugt ein
zusatzliches Energiequant, welches die Schwingungsamplitude des Oszillators erhoht. Der
tiefste Eigenzustand, der Grundzustand |0i, hat keine Energiequanten, aber E0 = ~/2 6=
0, eine Folge des Heisenbergschen Unscharfe-Prinzips: das Teilchen kann nicht scharf bei
x = 0 verweilen; durch die Unscharfe wird hHi =
6 0 sogar im Grundzustand; die Energie
~/2 ist das Resultat der Nullpunktsschwingungen des Vakuums.
Schliesslich brauchen wir noch die Wellenfunktionen, zuerst |0i, das Vakuum. Nach
Konstruktion ist a|0i = 0. In der Ortsdarstellung bedeutet das
65
2 /2
mit
1
1 x2 /2
1
2
e
= p 2x ex /2
hx|1i = hx|a |0i = (x x )
4
2
2
s
2
2
x ex /2 .
=
1
2
(x x )n ex /2 ,
hx|ni = p
2n n!
1
2
= p
Hn (x) ex /2 ,
2n n!
H0
H1
H2
H3
H4
=
=
=
=
=
1,
2x,
(2x)2 2,
(2x)3 6(2x),
(2x)4 12(2x2 ) + 12.
(6.1.32)
(6.1.33)
h
E0
q
6.2
Klassischer Limes
F
ur grosse Energien En nahert sich die quantenmechanische Losung der klassischen Losung an (klassischer Limes). Die klassische Aufenhaltswahrscheinlichkeit und Energie sind
66
gegeben durch
1
p
,
2q0 1 (x/q0 )2
1
m 2 q02 ,
=
2
Wkl =
Ekl
(6.2.1)
W kl
1
10
0
q0
q0
67
Physikalische Systeme besitzen haufig offensichtliche Symmetrien, zum Beispiel die Rotationssymmetrie eines Zentralkraftproblems. Mathematisch gesehen bedeutet dies, dass
das System invariant unter der Wirkung einer Symmetriegruppe G ist. Eine solche Symmetrie muss auch in der Quantentheorie weiterhin vorhanden bleiben; dabei muss die
Symmetrie mit den Strukturen der Quantenmechanik (namlich der Vektorraumstrukur
und dem inneren Produkt) kompatibel sein. Die einfachste Art, dies zu erreichen (wie
wir spater sehen werden, ist das jedoch nicht der allgemeinste Fall) ist dadurch realisiert,
falls die Gruppe G auf dem Hilbertraum H durch unitare Transformationen wirkt. Diese
respektieren namlich sowohl die Vektorraumstruktur, als auch das Skalarprodukt.
7.1
Unit
are Darstellungen
Eine unit
are Darstellung einer Gruppe G auf dem Hilbertraum H ist ein (Gruppen)Homomorphismus von G auf die Gruppe der unitaren Operatoren auf H. Konkret bedeutet
dies, dass wir f
ur jedes g G einen unitaren Operator U(g) haben, so dass
U(g1 ) U(g2 ) = U(g1 g2 ) .
(7.1.1)
(7.1.2)
(7.1.3)
f
ur jedes , H und g G. Wegen der Definition des adjungierten Operators folgt aus
dieser letzten Bedingung (zusammen mit der Gruppen-Homomorphismus Eigenschaft),
dass
U(g) = U(g)1 = U(g 1 ) .
(7.1.4)
Jede Gruppe besitzt die triviale Darstellung auf dem ein-dimensionalen Vektorraum
C, wobei jedes Gruppenelement g auf den Identitatsoperator abgebildet wird, U(g) = 1.
Ein interessanteres Beispiel ist die Darstellung der Rotationsgruppe auf dem Raum der
Wellenfunktionen L2 (R3 ). Die Gruppe der Rotationen des R3 ist die spezielle orthogonale
Gruppe SO(3), die durch die reellen 3 3 Matrizen R mit Determinante 1 beschrieben
wird. [Man zeigt leicht, dass die Menge dieser Matrizen eine Gruppe bilden; ferner ist klar,
dass jede dieser Matrizen auf den Vektoren des R3 eine Rotation definiert.] Eine unitare
Darstellung auf L2 (R3 ) ist dann durch
(U(R)) (~x) = (R1 x)
(7.1.5)
definiert. Man rechnet leicht nach, dass dies einen Homomorphismus der Rotationsgruppe
in die Gruppe der unitaren Operatoren auf L2 (R3 ) definiert.
68
(7.1.6)
d.h. sie entsprechen gerade den geraden bzw. ungeraden Funktionen. In vielen Fallen ist
jedoch der Hamiltonoperator eine gerade Funktion von x; dies bedeutet, dass die Operatoren H und P miteinander vertauschen. Dann kann man gemeinsame Eigenfunktionen
zu H und P finden das ist das, was wir zuvor bei der Analyse des Teilchens im Potentialtopf gemacht hatten.
7.2
Im Fall einer kontinuierlichen Gruppe (wie zum Beispiel der Rotationsgruppe) ist es oft
bequem, statt der Gruppentransformationen die infinitesimalen Transformationen zu betrachten. Eine kontinuierliche Gruppe ist (falls sie ein paar technische Bedingungen erf
ullt)
eine Lie Gruppe; die infinitesimalen Transformationen einer Lie Gruppe bilden eine Lie
Algebra, die f
ur das weitere wichtig sein wird. Um dieses Konzept zu verstehen betrachten wir wiederum die Gruppe der Rotationen. Sei R(t) eine differenzierbare Kurve von
Rotationen in SO(3) durch R(0) = 1. Eine infinitesimale Rotation ist dann
d
= R(t)
.
(7.2.1)
dt
t=0
[Wir stellen uns hier R(t) als eine orthogonale 3 3 Matrix vor; is daher auch eine 3 3
Matrix.] Die Menge dieser infinitesimalen Rotationen bilden einen (reellen) Vektorraum,
da
d
,
(7.2.2)
1 1 + 2 2 = R1 (1 t)R2 (2 t)
t=0
dt
wobei i R. Weiterhin gilt
d
1
1
R1 R = RR1 (t)R
,
(7.2.3)
dt
t=0
wobei R SO(3) liegt. Schliesslich haben wir
[1 , 2 ] =
d
R1 (t)2 R1 (t)1
.
dt
t=0
(7.2.4)
Der Kommutator zweier infinitesimaler Rotationen ist also wiederum eine infinitesimale
Rotation; der Vektorraum der infinitesimalen Rotationen schliesst daher unter der LieKlammer (und bildet damit eine Lie Algebra). Diese Lie Algebra bezeichnen wir mit
so(3).
69
Wie schon oben angedeutet stellen wir uns die infinitesimalen Rotationen wiederum
als Matrizen vor. Da R(t) orthogonal ist, d.h. R(t)T R(t) = 1 folgt durch Ableiten nach t,
dass jedes anti-symmetrisch ist, d.h.
T + = 0 .
(7.2.5)
Die Lie Algebra so(3) besteht also aus antisymmetrischen rellen 3 3-Matrizen. Jede
solche Matrix ist von der Form
0 3 2
0 1 ,
(~ ) = 3
(7.2.6)
2 1
0
d.h. (~ )~x = ~ ~x mit ~ = (1 , 2 , 2 ) R3 . Man kann leicht einsehen, dass so(3) alle
solche Matrizen enthalt; die Vektorraum-Dimension von so(3) ist also dimR so(3) = 3.
Zum Beispiel konnen wir die Basisvektoren ben
utzen
i := (~
ei ) ,
i = 1, 2, 3 ,
(7.2.7)
wobei {~
ei } die Standardbasis f
ur R3 ist.
F
ur ~ = e, mit |e| = 1, ist e(~ )t = R(e, t) (als Matrix aufgefasst) gerade die
(7.2.8)
wobei R eine beliebige Rotation in SO(3) ist. Wegen (7.2.4) folgt daher (wahle R = R(t)
und leite nach t ab), dass
[(~1 ), (~2 )] = (~1 ~2 ) .
(7.2.9)
(und zyklisch,)
(7.2.10)
Insbesondere ist dies der Fall, falls H endlich-dimensional ist, da dann U(R) wiederum
eine Matrix ist und die ganze Analyse wie zuvor durchgef
uhrt werden kann; im unendlichdimensionalen Fall erhalt man durch diese Vorschrift typischerweise unbeschrankte Operatoren, und man muss mit der Definition ein wenig aufpassen. Wie wir jedoch gleich
sehen werden sind alle interessanten Darstellungen der SO(3) endlich-dimensional.
70
Die durch (7.2.11) definierte Abbildung ist in der Tat eine Darstellung von so(3), d.h.
ein Vektorraum-Homomorphismus von so(3) in den Vektorraum der Operatoren auf H,
die die Lie Algebra Klammer erhalt, d.h.
U([1 , 2 ]) = [U(1 ), U(2 )] .
(7.2.12)
Jene letzte Eigenschaft folgt aus (7.2.4). Die Bedingung, dass die Darstellung U(R) unitar
ist, impliziert nun, dass
U() = U() ,
(7.2.13)
3
X
(7.2.14)
Mi i ,
i=1
(und zyklisch)
(7.2.15)
lauten. Wir nennen eine Darstellung der Lie Algebra unitar, wenn sie (7.2.13) erf
ullt, d.h.
wenn die zugehorige Darstellung der Lie Gruppe unitar ist.
Die unitare Darstellung der Lie Gruppe SO(3) ist eindeutig durch die Darstellung
ihrer Lie Algebra so(3), d.h. durch die Darstellung der Drehimpulsoperatoren, bestimmt.
In der Tat kann man aus letzterer Darstellung durch Exponenzieren die Darstellung der
gesamten Lie Gruppe rekonstruieren. Wir werden uns daher im folgenden haufig auf die
Darstellung der Lie Algebra (d.h. der Drehimpulsoperatoren) beschranken.
7.3
Sei H der Hilbertraum, auf dem eine unitare Darstellung der Gruppe G definiert ist. [D.h.
wir haben f
ur jedes g G einen unitaren Operator U(g).] Wir nennen diese Darstellung
irreduzibel, falls die einzigen Unterraume K von H, die unter der Wirkung von U(g)
f
ur alle g G auf sich selbst abgebildet werden, U(g)K K, einfach nur K = {0} und
K = H sind. Andernfalls nennen wir H reduzibel.
F
ur unitare Darstellungen kann man jede reduzible Darstellung in irreduzible Darstellungen zerlegen. Sei K also ein nicht-trivialer Unterraum von H, der unter der Wirkung
von G invariant ist. Dann gilt
H = K K ,
(7.3.1)
wobei
K = { H : h|i = 0 f
ur alle K} .
71
(7.3.2)
Wegen der Linearitat des Skalarproduktes ist dann auch K ein Unterraum von H. Weiterhin impliziert die Unitaritat von U nun, dass auch K unter der Wirkung von G
invariant ist: f
ur jedes K gilt namlich
h|U(g)i = hU(g 1 )|i = h |i = 0
(7.3.3)
falls K . Daher liegt also auch U(g) in K . Wir konnen also sukzessive die Darstellung H (falls sie reduzibel ist) in kleinere Unterdarstellungen zerlegen; diese sind entweder
irreduzibel oder konnen auf gleich Weise weiterbehandelt werden. Zumindest im endlichdimensionalen (im unendlich-dimensionalen Fall braucht man so etwas wie das axiom of
choice) terminiert dieser Prozess und es folgt, dass man jede reduzible unitare Darstellung
als direkte Summe von irreduziblen Darstellungen schreiben kann. Im folgenden wollen
wir daher irreduzible Darstellungen untersuchen.
7.4
Wir wollen nun die Struktur der irreduziblen endlich-dimensionalen Darstellungen von
so(3) beschreiben. [Alle diese Darstellungen sind automatisch unitar. Ferner gilt, dass
es keine unendlich-dimensionale unitaren Darstellung der SO(3) gibt; alle unitaren Darstellungen von SO(3) sind daher also direkte Summen der endlich-dimensionalen Darstellungen, die wir jetzt beschreiben wollen.] Die Lie Algebra von SO(3) wird durch die
Generatoren Mi , i = 1, 2, 3 aufgespannt, wobei die Kommutatoren durch
[M1 , M2 ] = iM3
(und zyklisch)
(7.4.1)
(7.4.2)
Ausgedr
uckt durch diese Basis sind dann die Vertauschungsregeln
[M3 , M ] = M ,
[M+ , M ] = 2M3 .
(7.4.3)
(7.4.4)
(7.4.5)
Da dim H < , kann dieses Argument nicht beliebig wiederholt werden: Es gibt also
einen Eigenwert j C mit Eigenvektor j , derart dass
M3 j = jj ,
M+ j = 0 .
(7.4.6)
(7.4.7)
f
ur m = j, j 1, . . .; somit ist
M3 m = mm .
(7.4.8)
M jk = 0 .
(7.4.9)
Falls
M+ m = m m+1
(7.4.10)
(was f
ur m = j zutrifft mit j = 0), so gilt auch
M+ m1 = M+ M m = [M+ , M ]m + M M+ m
= (2m + m )m m1 m .
Es folgt induktiv
m = j(j + 1) m(m + 1) = (j m)(j + 1 + m) .
(7.4.11)
(7.4.12)
Jeder irreduziblen Darstellung entspricht damit ein solches j. Umgekehrt verifiziert man,
dass M3 , M , durch (7.4.7, 7.4.8, 7.4.10) auf Basisvektoren j , . . . , j definiert, (7.4.3)
erf
ullen und somit eine Darstellung Dj der so(3) liefern.
Die Darstellung kann auch durch den Wert des Casimiroperators charakterisiert
werden. F
ur den Fall von so(3) vertauscht der Operator
~ 2 = M12 + M22 + M32 = M32 + 1 (M+ M + M M+ )
M
2
(7.4.13)
(7.4.14)
(7.4.15)
(7.4.16)
Wir haben also damit gezeigt: die endlich dimensionalen irreduziblen Darstellung der
so(3), sind parametrisiert durch (7.4.12) mit dim Dj = 2j + 1. Es gilt (7.4.16).
73
|j, ji :=
gegeben. Diesbez
uglich ist
j
,
kj k
m |j, mi := M |j, m + 1i
(7.4.18)
(7.4.19)
F
uhrt man umgekehrt ein Skalarprodukt ein, indem man die Basis {|j, mi}jm=j als orthonormiert erklart, so ist die Darstellung Dj unitar.
Die triviale Darstellung entspricht nat
urlich gerade D0 . Die fundamentale (oder definierende) Darstellung ist auf H = R3 (bzw. auf H = C3 ) definiert, wobei U(R) = R und
U() = ist. Sie ist irreduzibel, hat Dimension d = 3 und ist daher isomorph zu D1 .
Dasselbe gilt f
ur die adjungierte Darstellung auf H = so(3) (bzw. auf ihrer Komplexifizierung), die durch
U(R) = R R1
bzw.
U(1 )2 = [1 , 2 ]
(7.4.20)
definiert ist.
Eine andere Klasse von Darstellungen sind die Unterraume der Funktionen auf der
Kugel, die durch Kugelfunktionen erzeugt werden. Dazu betrachten wir den Raum aller
(komplex-wertigen) Funktionen, die auf der Kugeloberflache definiert sind. Jede solche
Funktion konnen wir als Funktion von und schreiben, wobei den Breitengrad und
den Langengrad beschreibt. Auf dem Vektorraum dieser Funktionen wirkt die Drehgruppe
wie durch (7.1.5) beschrieben; in den obigen Kugelkoordinaten gilt dann
i
.
(7.4.21)
+ i cot
,
M = e
M3 = i
[In der Tat rechnet man leicht nach, dass diese Differentialoperatoren den obigen so(3)
Vertauschungsregeln gen
ugen.] Der Vektorraum aller auf der Kugel definierten Funktionen
ist jedoch bez
uglich dieser Wirkung nicht irreduzibel. Um die invarianten Unterraume zu
beschreiben f
uhren wir die Kugelfunktionen Yl,m durch
s
2l + 1 (l m)! m
Yl,m (, ) =
Pl (cos ) eim
(7.4.22)
4 (l + m)!
ein, wobei m = l, l + 1, . . . , l 1, l ist. Hier ist Plm die assoziierte Legendrefunktion,
die die verallgemeinerte Legendre Gleichung
d
m2
2 dP (z)
(1 z )
+ l (l + 1)
P (z) = 0
(7.4.23)
dz
dz
1 z2
74
z
)
(z 2 1)l
l
l+m
2 l!
dz
Plm (z) =
(7.4.24)
gegeben.
Man kann nun leicht sehen, dass
M3 Yl,m = mYl,m
(7.4.25)
gilt. Weiterhin zeigt man leicht, dass Yl,l tatsachlich von M+ vernichtet wird; in der Tat
ist
Yl,l Pll (cos )eil
d2l
(1 z 2 )l/2 2l (z 2 1)l eil
dz
l
il
sin e .
(7.4.26)
il
M+ sin e
= e
+ i cot
sinl eil
(7.4.27)
Mit ein wenig mehr Aufwand kann man schliesslich zeigen, dass
M Yl,m = Cl,m, Yl,m1 .
(7.4.28)
7.5
Uberlagerungsgruppe
SU(2). Diese Subtilitat ist f
ur das Weitere nicht unwichtig; sie soll
deshalb ein wenig genauer erklart werden.
Die Gruppe der unitaren (komplexen) 22 Matrizen mit Determinante +1 wird SU(2)
genannt; jedes Element von SU(2) kann man also schreiben als
a
b
,
wobei
|a|2 + |b|2 = 1 .
(7.5.1)
g=
b a
75
Wie zuvor kann man daraus ihre Lie Algebra bestimmen (die wir mit su(2) bezeichnen
wollen tatsachlich stimmt su(2) wie wir gleich sehen werden mit der Lie Algebra von
SO(3) u
ur die
berein): sie besteht aus alle komplexen 2 2 Matrizen A, f
A + A = 0 ,
sp(A) = 0 ,
(7.5.2)
d.h. aus allen anti-hermiteschen, spurlosen Matrizen. Ein beliebiges Element von su(2)
kann man mit Hilfe der Pauli-Matrizen entwicklen; diese sind die spurlosen hermiteschen
Matrizen, die explizit durch
1 0
0 i
0 1
(7.5.3)
,
3 =
,
2 =
1 =
0 1
i 0
1 0
definiert. In der Tat kann nun jedes A su(2) als
i
A A(~a) =
2
a3
a1 ia2
a1 + ia2
a3
iX
i
j aj ~ ~a
2 j=1
2
(7.5.4)
geschrieben werden. Insbesondere hat also su(2) (reelle) Dimension 3, und eine nat
urliche
Basis ist
j
j = 1, 2, 3 .
(7.5.5)
Aj = i ,
2
Die Matrizen (7.5.3) erf
ullen
i j = ij + iijk k ,
(7.5.6)
wobei 123 = +1 und ijk total antisymmetrisch ist. In Vektorschreibweise gilt also
(~ ~a)(~ ~b) = (~a ~b)1I + i~ (~a ~b) .
(7.5.7)
(7.5.8)
Damit ist
bzw.
[A1 , A2 ] = A3
(und zyklisch).
(7.5.9)
Die Lie-Algebren su(2) und so(3) sind also isomorph, wobei der Isomorphismus durch
su(2) so(3) ,
A(~ ) 7 (~ ) ,
(7.5.10)
d.h. durch Aj 7 j gegeben ist. Die irreduziblen Darstellungen der su(2) sind damit die
Dj aus Kapitel 7.4
Obgleich die Lie Algebren von SU(2) und SO(3) u
bereinstimmen, sind die beiden
Gruppen verschieden. Um ihre Relation genauer zu verstehen definieren wir f
ur jeden
Vektor ~x R3 die Matrix
3
X
x3
x1 ix2
j
,
(7.5.11)
x =
x j =
x1 + ix2
x3
j=1
76
1
sp(
x j ) ,
2
(7.5.12)
wie man leicht nachrechnet. Falls ~x ein reeller Vektor ist, dann ist x eine hermitesche
Matrix; umgekehrt f
uhrt jede hermitesche Matrix x vermoge (7.5.12) zu einem reellen
Vektor. Schliesslich rechnet man leicht nach, dass die Determinante von x gerade mit dem
Skalarprodukt von x u
bereinstimmt,
det x = ~x ~x .
(7.5.13)
F
ur jedes Element A von SU(2) betrachten wir nun die Abbildung
x 7 x = A x A .
(7.5.14)
Als Abbildung von ~x 7 ~x aufgefasst, ist dies eine lineare Abbildung, die reelle Vektoren
auf reelle Vektoren abbildet. [Falls x hermitesch ist, dann ist auch A x A hermitesch.]
Wegen (7.5.13) lasst diese Transformation die Lange von ~x invariant, und definiert daher
eine Drehung. Diese Konstruktion definiert daher einen Gruppenhomomorphism
SU(2) SO(3) .
(7.5.15)
Es ist relativ offensichtlich, dass der Kern dieser Transformation gerade 12 SU(2) ist.
Ferner ist der Homomorphismus surjektiv, und daher haben wir
SU(2)/{1} SO(3) .
(7.5.16)
Offensichtlich wirkt in diesem Fall die Matrix 1 nicht trivial, und daher definiert diese
Darstellung nicht eine Darstellung auf der Quotientengruppe SO(3)=SU(2)/Z2 .
Im allgemeinen gilt in SO(3)
0 1 0
R(~e3 , ) = exp i 1 0 0 = eiM3 ,
(7.5.17)
0 0 0
und daher ist insbesondere auf jedem Darstellungsvektor |j, mi
(7.5.18)
Wir sehen also, dass die Rotation um 2 (die in SO(3) die Identitatsabbildung sein sollte)
nur in Darstellungen, f
ur die j ganzzahlig ist, trivial dargestellt ist: die Darstellungen f
ur
die j ganz-zahlig ist, sind umgekehrt durch die Kugelfunktionen explizit realisiert. Diese
sind also tatsachlich alle Darstellungen von SO(3)!
Falls andererseits j halb-ganz ist, dann wird die Rotation um 2 (die ein nicht-triviales
7.6
Projektive Darstellungen
Bis anhin haben wir nur unitare Darstellungen der Symmetriegruppe studiert; im allgemeinen ist aber nicht klar, ob eine physikalische Symmetrie tatsachlich auf diese Weise
auf dem Hilbertraum realisiert sein muss. Eine mogliche Abschwachung ist relativ direkt
einsichtig: wie wir zuvor in Kapitel 4 erklart haben, wird ein physikalische Zustand nicht
durch einen Vektor im Hilbertraum beschrieben, sondern lediglich durch einen Strahl,
(7.6.1)
wobei c(g1 , g2 ) eine Phase ist. Diese Phasen konnen nicht beliebig gewahlt sein: da die Operatorprodukte immer noch assoziativ sein m
ussen, gilt notwendigerweise die sogenannte
Kozykelbedingung
c(g1 , g2 g3 ) c(g2 , g3 ) = c(g1 g2 , g3 ) c(g1 , g2 ) .
(7.6.2)
Es gibt eine relativ offensichtliche Art, solche Phasen k
unstlich zu erzeugen: sei U(g)
eine echte Darstellung von g, dann definiere
U(g)
= U(g) c(g) ,
(7.6.3)
wobei c(g) f
ur jedes g G eine Phase ist. Man kann leicht nachrechnen, dass U(g)
dann
eine projektive Darstellung definiert, wobei
c(g, h) =
c(g) c(h)
.
c(gh)
(7.6.4)
Umgekehrt, kann man versuchen aus einer projektiven Darstellung durch diese Modifikation wiederum eine echte Darstellung zu erhalten. Manchmal wird das aber nicht moglich
sein: man spricht dann von einer (echt) projektiven Darstellung der Gruppe.
78
Man kann zeigen, dass die projektiven Darstellung einer Lie Gruppe G in eins zu eins
von
7.7
ohne Storung
mit Storung
Aus der Beobachtung der Spektren findet man Falle mit 2j +1 gerade, also j halbganz.
Um diese Entartung zu erklaren f
uhrt man einen zusatzliche Freiheitsgrad ein, den man
den Spin nennt. Der Hilbertraum eines einzelnen Elektrons ist dann nicht einfach nur
L2 (R3 ), sondern (Pauli 1924)
H = L2 (R3 ) C2 .
(7.7.1)
Auf diesem Raum wirkt dann V SU(2) gemass
U(V ) = U0 (R(V )) V ,
(7.7.2)
wobei hier U0 (R) die Darstellung von R SO(3) auf L2 (R3 ) ist (siehe (7.1.5). F
ur den
Fall des Elektrons transformiert sich der interne Freiheitsgrad also in der Darstellung D 1 .
2
79
Mj =
Damit ist
M+ =
0 1
0 0
j
.
2
M =
(7.7.3)
0 0
1 0
(7.7.4)
(7.7.5)
7.8
Wigners Theorem
Wie wir in Kapitel 7.6 gesehen haben, muss die Darstellung einer Symmetriegruppe auf
dem Hilbertraum der Quantenmechanik nicht notwendigerweise eine echte Darstellung
sein; es gen
ugt, wenn die Darstellung projektiv ist. Bis anhin haben wir jedoch immer
implizit angenommen, dass die Darstellung unitar ist. Vom Gesichtspunkt der Quantenmechanik kommt es jedoch lediglich darauf an, dass die Wahrscheinlichkeiten
|h|i|2
h|i h|i
(7.8.1)
unter der Symmetrieoperation invariant sind. Offensichtlich gilt dies, falls die Symmetrie
durch unitare Operatoren dargestellt ist, da dann jedes Skalarprodukt h|i invariant ist.
Im allgemeinen gen
ugt es jedoch auch, wenn die Symmetrie durch anti-unit
are Operatoren U(g) dargestellt wird. Ein anti-unitarer Operator hat die Eigenschaft, dass
hU(g)|U(g)i = h|i = h|i .
(7.8.2)
In der Tat hat man damit alle Moglichkeiten erschopft; das ist der Inhalt des Wignerschen Theorems: Sei G eine Gruppe die auf den Zustanden des Hilbertraums H so
wirkt, dass die Wahrscheinlichkeiten (7.8.1) erhalten bleiben. Dann kann man die Wirkung von g G als 7 U(g) schreiben, wobei jedes U(g) ein unitarer oder anti-unitarer
Operator ist, und
U(gh) = c(g, h) U(g) U(h)
(7.8.3)
80
(T )(t, x) = (t,
x) .
(7.8.4)
Die komplexe Konjugation der Wellenfunktion ist hier notig, damit weiterhin Losungen
der Schrodingergleichung auf Losungen der Schrodingergleichung abgebildet werden. Da
T diese Konjugation involviert, ist der Operator auf dem Raum der Wellenfunktionen
anti-unitar und nicht unitar. Ein anderes wichtiges Beispiel ist die Ladungskonjugation,
die auch durch einen anti-unitaren Operator auf dem Hilbertraum dargestellt wird.
Die meisten Symmetrien werden jedoch durch unitare Operatoren dargestellt; dies ist
zum Beispiel eine Folge des Weylschen Lemmas: der Operator U(g 2 ) ist immer unitar
f
ur alle g G. In der Tat folgt aus (7.8.3), dass
U(g 2 ) = c(g, g)U(g) U(g) .
(7.8.5)
81
Das Wasserstoffatom
Ein anderes wichtiges Problem, das exakt gelost werden kann, ist das Wasserstoffatom.
Das Wasserstoffatom ist ein 2-Korperproblem (Atomkern plus Elektron) mit rotationssymmetrischem 2-Korperpotential. Es ist also einfach die quantenmechanische Version des
Keplerproblems.
8.1
Relativkoordinaten
p21
p2
+ 2 + V (|~r1 ~r2 |) ,
2m1 2m2
(8.1.1)
wobei m1 und m2 die beiden Massen sind, und die Wechselwirkung durch das Potential
V (|~r1 ~r2 |) beschrieben ist. Wir definieren die Schwerpunkts- und Relativ-Koordinaten
~ und ~x
X
~ = m1~r1 + m2~r2 .
(8.1.2)
~x = ~r1 ~r2 , X
m1 + m2
Entsprechend definieren wir die konjugierten Impulse
p~ =
m2 ~p1 m1 p~2
,
m1 + m2
P~ = ~p1 + p~2 .
(8.1.3)
Ausserdem definieren wir die reduzierte Mass sowie die totale Masse M durch
=m=
m1 m2
,
m1 + m2
M = m1 + m2 .
(8.1.4)
[Xi , Pj ] = i~ ij .
(8.1.6)
Auf dem Raum der Wellenfunktionen (L2 (R6 )) konnen diese Operatoren also durch die
Differentialoperatoren
~ ~x ,
~ ~
~p = i~
P~ = i~
(8.1.7)
X
82
p2
P2
+
+ V (|~x|) .
2M
2
(8.1.8)
Die zeit-abhangige Schrodingergleichung kann nun dadurch separiert werden, dass wir den
Ansatz machen
~ X
~ i~K 2 t/2M i E t/~
~ = eiK
K~ (~x, X)
e
e
(~x) .
(8.1.9)
Das resultierende Relativproblem ist dann
2
p
+ V (r) (~x) = E(~x) ,
2
(8.1.10)
wobei wir r = |~x| geschrieben haben. Die totale Energie ist die Kombination Etot =
E + ~2 K 2 /2M.
8.2
Coulomb-Potential
~ = i~x ~x .
M
(8.2.2)
Wie wir in Kapitel 7.2 gesehen hatten, sind die infinitesimalen Rotation von der Form
~x 7 ~x + i~ ~x. Auf dem Raum der Funktionen wirken diese also wie (siehe (7.1.5))
f (~x) 7 f (~x i~ ~x) = f (~x) i(~ ~x) ~x f (~x) = f (~x) i~ (~x ~x f ) .
(8.2.3)
~ ist daher also wie (8.2.2) definiert. Bis auf einen Faktor stimmt
Der Rotationsgenerator M
~ mit dem Drehimpulsgoperator L
~
also M
~ = ~x p~ = i ~ ~x ~x = ~ M
~
L
(8.2.4)
u
berein.
Wie wir in Kapitel 7.4 gesehen haben, spielt in der Analyse der Darstellungen der
~ 2 eine wichtige Rolle. Auf dem Raum der WelRotationsgruppe der Casimir-Operator M
lenfunktionen wirkt dieser Operator nun wie
~2 =
M
=
=
=
(8.2.5)
(8.2.6)
(8.2.7)
(8.2.8)
In Polarkoordinaten ist
xi i = ~x ~x = r
(8.2.9)
~ = ~x + r
M
= ~x2 + r 2 r ,
+r
r
r
r
2
(8.2.10)
bzw.
1 ~2
1 2
r 2M
.
(8.2.11)
2
r r
r
~ 2 gerade durch l(l+1) gegeben (siehe Kapitel 7.4),
Andererseits sind die Eigenwerte von M
wobei f
ur die Darstellungen auf dem Funktionenraum L2 (R6 ) nur l = 0, 1, . . . auftreten
konnen. (Diese Funktionen definieren ja immer eine echte Darstellung der Rotationsgruppe
SO(3).) In der Tat sind die Funktionen mit Eigenwert l(l + 1) gerade die Kugelfunktionen
Ylm (, ) mit m = l, l + 1, . . . , l 1, l. Wenn wir also den Ansatz
=
(8.2.12)
(8.2.13)
einfach
~2
2m
1 2
1
r(r) 2 l(l + 1)(r) + V (r)(r) = E(r) .
r r 2
r
(8.2.14)
gesetzt haben. Wir diskutieren nun das Verhalten der Losung bei r 0 und r .
84
8.2.1
F
ur r 0 f
uhrt (8.2.16) auf
l(l + 1)
u=0
r2
(8.2.18)
u(r) = ar l+1 + br l .
(8.2.19)
u +
(8.2.20)
(k =
) .
(8.2.21)
Insbesondere ist zu erwarten, dass die regulare Losung u(, r) von (8.2.16) f
ur r die
asymptotische Form
u(, r) a()eikr + b()eikr
(8.2.22)
besitzt. Damit u(, ) L2 (0, ) muss erstens < 0: wir legen dann k fest durch
(8.2.23)
k = i ,
= > 0 ,
also eikr = er . Zweitens muss b() = 0 sein, da die zugehorige Losung exponentiell
anwachst: die Eigenwerte ergeben sich als Nullstellen der Funktion b(). Dann reduziert
sich (8.2.22) auf
u(, r) a()er ,
(r ) .
(8.2.24)
Die bisherige Analyse war in der Tat unabhangig von der speziellen Form des Potentials
V (r). Wir behandeln nun den Fall (8.2.1) des Coulombpotentials
V (r) =
,
r
2mZe2
~2
(8.2.25)
u(r) = e
k=l+1
85
ck r k .
(8.2.26)
2
,
k+1
also ck C
(2)2k
,
k!
(8.2.28)
was auf
u(r) er Ce2r = Cer
(8.2.29)
f
uhrt. Die zugehorige Losung u(r) er e2r ist dann nicht quadratintegrierbar. [Die
Situation ist also genauso wie in der konventionellen Losung des harmonischen Oszillators
siehe Kapitel 6.1.1.] Falls (8.2.27) hingegen abbricht, d.h. falls f
ur ein n
cn 6= 0 ,
cn+1 = 0 ,
(8.2.30)
so ist die Losung eine Eigenfunktion. [Die resultierenden Polynome sind die sogenannten
zugeordneten Laguerre Polynome.] Die Bedingung daf
ur ist
n =
,
2n
d.h.
(n = l + 1, l + 2, . . .) ,
(8.2.31)
m(Ze2 )2 1
2 ,
(8.2.32)
2m
2~2
n
(Schrodinger 1926). Dies ist die von Bohr im Rahmen der alten Quantentheorie hergeleitete Formel f
ur die Energieniveaux siehe Kapitel 1.3.
Da es f
ur jedes l 2l + 1 verschiedene Eigenfunktionen Ylm gibt (die durch m = l, . . . , l
parametrisiert werden) tritt jedes dieser Energieniveaux mit der Multiplizitat
En =
Dn =
n1
X
~2
2n =
(2l + 1) = 2
l=0
n(n 1)
+ n = n2
2
(8.2.33)
auf. Die Quantenzahlen n, l, m heissen Haupt-Quantenzahl n, Neben- oder Bahndrehimpuls-Quantenzahl l, und magnetische Quantenzahl m. Die zusatzliche Entartung in l
ist eine Konsequenz der dynamischen Symmetrie im Coulombpotential. Wird der Spin
ber
ucksichtigt, (siehe spater) verdoppelt sich die Entartung, Dn = 2n2 .
Fundamental ist die Existenz eines energetisch tiefsten Zustands (l = 0, n = 1): die
Energie des H-Atom ist nach unten beschrankt und es ist damit stabil! Dies im Gegensatz zum klassischen H-Atom, wo das (beschleunigte) Elektron beliebig viel Energie durch
Ausstrahlung abgeben w
urde (vgl. Elektrodynamik). Die Wellenfunktion des Grundszustandes ist nach (8.2.26)
u(r) = e1 r r ,
1 =
86
me2
= 2 ,
2
~
(8.2.34)
a=
~2
me2
(8.2.35)
Der Bohr-Radius a ist der Radius des Atoms in der Bohrschen Theorie. Durch (8.2.32)
sind alle Eigenwerte von H gefunden (ohne Beweis). Damit ist aber (im Unterschied zum
harmonischen Oszillator) noch nicht das gesamte Spektrum (H) ausgeschopft, welches
(wie beim freien Teilchen) auch einen kontinuierlichen Anteil [0, ) besitzt:
0
Die Uberg
ange zwischen den Niveaus involvieren l = 1 (siehe spater) und gehoren
zu einer der Sequenzen in der Tabelle 1; siehe auch Abbildung 12 f
ur eine graphische
Darstellung der involvierten Zustande.
Endniveau n
n=1
n=2
n=3
n=4
n=5
Name
Energiebereich
Lyman
Balmer
Paschen
Bracket
Pfund
UV
sichtbar
IR
IR
IR
Tabelle 1: Uberg
ange zwischen Niveaus
n
5
4
3
Bracket
Pfund
Paschen
2
Balmer
1
Lyman
87
8.3
Dynamische Symmetrie
p2
,
2 r
(8.3.1)
r
2H 2
J~2 =
L + 2 .
~ = e ,
|J|
~ =0,
J~ L
(8.3.2)
Der Umstand, dass J~ konstant ist, impliziert, dass das Perihel sich nicht bewegt; dies ist
eine Besonderheit des 1/r-Potentials.
Aphel
L
Perihel
Abbildung 13: Der Runge-Lenz Vektor J~ zeigt vom Brennpunkt zum Perihel und ist eine erhaltene Grosse
im Keplerproblem. Storungen des
1/r-Potentials verursachen eine Dre~ die Beobachtung der Pehung von J;
riheldrehung des Merkurs war eine
der ersten experimentellen Bestatigungen der ART.
F
ur V 1/r 1+ resultiert eine Periheldrehung; ein ber
uhmtes Beispiel sind die Korrekturen zum Kepler Problem infolge der Allgemeinen Relativitats Theorie (ART); die
Periheldrehung des Merkurs war eine der ersten experimentellen Bestatigungen der ART.
88
(8.3.3)
Diese Definition ist so gewahlt, dass J~ weiterhin mit dem Hamiltonoperator vertauscht; die
Eigenzustande des Hamiltonoperators zu vorgegebenem Eigenvektor (d.h. die Losungen
der zeit-unabhangigen Schrodingergleichung zu festem E) transformieren sich also inein~
ander unter der Wirkung von J.
~ = 0, wobei L
~ der u
Es gilt noch immer J~ L
bliche Drehimpulsvektor
~ = ~x ~p
L
(8.3.4)
(8.3.5)
2H 2
i~ ijk Mk .
(8.3.8)
Wir betrachten zunachst den Fall, wo der Eigenwert E von H negativ ist. Dann reskalieren
wir Ji als
r
1
Ji
(8.3.9)
Ki =
2H ~
(da E negativ ist, ist der Vorfaktor reell, und Ki ist immer noch ein selbst-adjungierter
Operator), und erhalten dann die kombinierten Vertauschungsregeln
[Ki , Kj ] = iijk Mk ,
[Mi , Kj ] = iijk Kk ,
[Mi , Mj ] = iijk Mk .
89
(8.3.10)
(8.3.11)
(8.3.12)
und beide bilden gerade jeweils die Lie Algebra von su(2):
[Si , Sj ] = iijk Sk ,
[Di , Dj ] = iijk Dk .
(8.3.13)
Die Lie Algebra dieser Operatoren ist also gerade su(2) su(2), was tatsachlich zu der Lie
H =
~ 2 + 2D
~2+1
2~2 2S
.
(8.3.14)
~ 2 und D
~ 2 folgen wiederum aus der Darstellungstheorie
Die moglichen Eigenwerte von S
von su(2):
~ 2 = s(s + 1) ,
~ 2 = d(d + 1) ,
S
D
(8.3.15)
~ noch D
~ haben eine direkte Interpretation
wobei s und d halb-ganze Zahlen sind. [Weder S
als Rotationen im Raum; von daher ist nicht notwendigerweise klar, dass nur ganz-zahlige
~ M
~ = 0, und daher, dass
Spins erlaubt sind.] Weiterhin wissen wir, dass K
~ D)
~ (S
~ + D)
~ =S
~2 D
~2 = 0 .
(S
(8.3.16)
Dies impliziert daher also, dass s = d. Die moglichen Energieeigenwerte sind also
2
E = 2
2~ (4s(s + 1) + 1)
2
= 2
2~ (2s + 1)2
2
mit
n = 2s + 1 .
(8.3.17)
= 2 2
2~ n
Die moglichen Werte von n sind also n = 1, 2, 3, . . .. Die Dimension der Eigenraume ist
dimEn = (2s + 1)(2d + 1) = (2s + 1)2 = n2
(8.3.18)
in Ubereinstimmung
mit dem, was wir zuvor gesehen hatten. Wir haben also das gesamte
Spektrum (zusammen mit seinen Multiplizitaten) aus Symmetrie
uberlegungen bestimmen
konnen!
90
Drehimpulsaddition
In vielen Fallen ist ein Quantensystem aus mehreren Teilsystemen zusammengesetzt. Wir
wollen verstehen, wie sich der Zustandsraum eines sochen Quantensystems durch jene seiner Teilsysteme beschreiben lasst. Das relevante Konzept ist jenes eines Tensorproduktes.
9.1
Tensorprodukte
Der Hilbertraum eines Systems, das aus zwei Teilsystemen zusammengesetzt ist, ist das
Tensorprodukt der Hilbertraume der Teilsysteme
H = H(1) H(2) .
(9.1.1)
Ein einfaches Beispiel davon haben wir gerade gesehen: der Hilbertraum des 2-Korperproblems war L2 (R6 ) = L2 (R3 ) L2 (R3 ). Anders ausgedr
uckt: der Hilbertraum der 2Teilchen-Wellenfunktionen (~x1 , ~x2 ) wird aufgespannt durch die (Tensor-)produkte
(1) (~x) (2) (~x2 ) = (1) (2) (~x1 , ~x2 ) .
(9.1.2)
Eine einfache Art, sich das Tensorprodukt vorzustellen, ist wie folgt. Das Tensorprodukt
ist der Vektorraum der endlichen Summen
l
X
j=1
vj wj
vj H(1) ,
wj H(2) .
(9.1.3)
vj wj =
(vj ) wj =
vj (wj )
(9.1.4)
j=1
j=1
j=1
(2)
ei ej ,
i = 1, 2, . . . ,
j = 1, 2, . . .
(9.1.5)
gegeben ist. Insbesondere ist also im endlich-dimensionalen Fall die Dimension des Tensorproduktes gerade das Produkt der Dimensionen der beiden Vektorraume.
Im Kontext der Quantenmechanik interessieren wir uns f
ur Hilbertraume. Das Skalarprodukt f
ur das Tensorprodukt wird erklart durch
hv1 w1 |v2 w2 i = hv1 |v2 i hw1|w2 i .
91
(9.1.6)
Dieses Produkt wird dann linear bzw. anti-linear auf beliebige endliche Summen solcher
Ausdr
ucke fortgesetzt. Man zeigt leicht, dass das resultierende Skalarprodukt positivdefinit ist, falls die Skalarprodukte der beiden Faktortheorien diese Eigenschaft besitzen.
(1)
(2)
Falls die ei und ei eine Orthonormalbasis (ONB) f
ur H(1) und H(2) definieren, dann
definiert (9.1.5) gerade auch eine ONB f
ur das Tensorprodukt.
9.2
Wir betrachten nun den Fall, dass beide separaten Hilbertraume eine (projektive) Darstellung der Drehgruppe tragen. Zum Beispiel ist das f
ur das 2-Korperproblem der Fall, da
2
3
der Hilbertraum H = L (R ) eine (reduzible) Darstellung der Drehgruppe tragt diese
ist durch (7.1.5) definiert. In vielen Situationen ist jedoch der Hamiltonoperator des Gesamtsystems nicht unter diese separaten Rotationen invariant (d.h. der Hamiltonoperator
vertauscht nicht mit den separaten infinitesimalen Rotationsgeneratoren), sondern lediglich unter der gemeinsamen Rotationen. In diesem Fall vertauscht der Hamiltonoperator
also nur mit den infinitesimalen Generatoren der gemeinsamen Rotationen. Wir wollen
nun verstehen, wie diese explizit definiert sind.
Wir betrachten die Situation, wo die beiden Hilbertraume H(1) und H(2) Darstellungen
(9.2.1)
(9.2.2)
wobei i = 1, 2. Auf dem Tensorprodukt H(1) H(2) haben wir deshalb eine Darstellung
der Produktgruppe SU(2) SU(2) durch
(g, h) 7 U (1) (g) U (2) (h) ,
wobei die rechte Seite auf ein Element 1 2 als
U (1) (g) U (2) (h) (1 2 ) = (U (1) (g)1 ) (U (2) (h)2 )
(9.2.3)
(9.2.4)
wirkt. [Die Fortsetzung auf ein beliebiges Element des Tensorproduktraumes erfolgt unter Ben
utzung der Linearitat. Man rechnet leicht nach, dass dadurch tatsachlich eine
Darstellung der Produktgruppe definiert wird.]
Wir interessieren uns f
ur die gemeinsamen Rotationen, also f
ur die Darstellung der
Gruppe SU(2) auf dem Tensorprodukt durch
g 7 U (1) (g) U (2) (g) .
(9.2.5)
Die zugehorigen infinitesimalen Generatoren werden wir zuvor definiert: wir betrachten
eine Kurve g(t) mit g(0) =id, und nehmen die Ableitung nach t am Punkt t = 0, wobei
d
= g(t)
(9.2.6)
dt
t=0
92
der entsprechende Lie-Algebra Generator ist. Auf dem Tensorprodukt gilt dann
d
(1)
(2)
(1)
(2)
U (g(t)) U (g(t))
U U
() =
dt
t=0
d (1)
(2)
=
U (g(t)) U (g(t)) t=0
dt
t=0
d
(2)
(1)
+ U (g(t))t=0 U (g(t))
dt
t=0
= U (1) () 1 + 1 U (2) () .
(9.2.7)
Die Lie-Algebra wirkt also dadurch auf dem Tensorprodukt, dass man die separaten Wirkungen auf den einzelnen Faktoren addiert. [Manchmal fasst man das als Komultiplikation auf der Lie-Algebra g auf
() = 1 + 1 ,
(9.2.8)
Wie wir gesehen haben konnen wir alle unitaren Darstellungen von SU(2) als direkte
Summe der irreduziblen unitaren Darstellungen beschreiben; diese sind die Darstellungen Dj mit j = 0, 21 , 1, . . . der Dimension (2j + 1). Falls wir zwei solche Darstellungen
zusammentensorieren erhalten wir nach dem, was wir gerade erklart haben, wiederum
eine Darstellung der Lie-Algebra su(2). Im allgemeinen ist jedoch diese Darstellung nicht
irreduzibel. Es ist daher von einigem Interesse zu verstehen, wie sich diese Darstellung
wiederum in irreduzible Darstellungen zerlegen lasst. Dieses Problem kann in volliger Allgemeinheit gelost werden (Clebsch-Gordon Reihe). Bevor wir aber den allgemeinen Fall
behandeln wollen, ist es instruktiv, einen Spezialfall zu studieren.
9.3
Addition zweier j =
1
2
Darstellungen
Der einfachste nicht-triviale Fall tritt auf, falls beide Darstellungen gerade die 2-dimenur D 1 mit
sionale j = 12 Darstellung ist. Wir wahlen die u
bliche Basis f
2
M3 | i = 12 | i
M+ | i = 0
M | i = | i
M3 | i = 12 | i
M+ | i = | i
M | i =
0.
93
(9.3.1)
p
[Zur Erinnerung: der Koeffizient in der M Relation ist j(j + 1) m(m 1) siehe
(7.4.19).] Das Tensorprodukt ist also 4-dimensional, und eine Basis besteht aus
| , i,
| , i,
| , i,
| , i .
(9.3.2)
Wir analysieren nun die Wirkung der Lie-Algebra Generatoren auf diesen Vektoren. Aus
der Definition (9.2.7) folgt sofort, dass alle vier Vektoren Eigenvektoren von M3 sind,
wobei
M3 | , i
M3 | , i
M3 | , i
M3 | , i
=
=
=
=
| , i
0
0
| , i .
(9.3.3)
(9.3.4)
Der Vektor | , i ist also der Hochstgewichtszustand einer Darstellung D1 : diese besteht
aus den Vektoren
|1, 1i = | , i ,
1
(9.3.5)
|1, 0i = M | , i
2
1
|1, 1i =
M M | , i .
2
p
j(j + 1) m(m 1)|j, m 1i ist. Dies bestimmt die
[Zur Erinnerung: M |j, mi =
obigen Normierungsfaktoren.] Explizit gilt daher also
1
|1, 0i = (| , i + | , i) ,
2
sowie
(9.3.6)
1
(| , i + | , i) = | , i .
(9.3.7)
2
Schliesslich gilt dann auch automatisch, dass M |1, 1i = 0.
In dem 4-dimensional Tensorproduktraum haben wir also bisher die 3-dimensionale
Darstellung D1 (9.3.5) wiedergefunden. Offensichtlich spannen diese 3 Vektoren aber noch
nicht das ganze Tensorprodukt auf, da in D1 nur eine bestimmte Linearkombination der
beiden Vektoren mit M3 = 0 auftritt. Der zu |1, 0i orthogonale Vektor muss daher in einer
anderen Darstellung liegen. [Wie wir in Kapitel 7.3 gesehen haben ist das orthogonale
Komplement einer Unterdarstellung im unitaren Fall auch eine Unterdarstellung. Mit der
Definition des Skalarproduktes auf dem Tensorprodukt kann man leicht einsehen, dass
das Tensorprodukt zweier Darstellungen unitar ist, falls die beiden Darstellungen es sind.
Das ist der Fall, der hier f
ur uns relevant ist.]
|1, 1i =
94
Man rechnet leicht nach, dass der zu D1 (bzw. zu |1, 0i) orthogonale Vektor gerade
1
|0, 0i = (| , i | , i)
2
(9.3.8)
ist. Der Vektor M+ |0, 0i hat dann M3 Eigenwert 1, aber da er orthogonal zu D1 liegt
kann er nicht mit dem Vektor |1, 1i u
bereinstimmen. Da es in dem Tensorprodukt keinen
anderen Vektor mit M3 Eigenwert +1 gibt, muss also M+ |0, 0i verschwinden. Nat
urlich
kann man das auch direkt nachpr
ufen
1
M+ |0, 0i = (M+ | , i M+ | , i)
2
1
= (| , i | , i) = 0 .
2
(9.3.9)
Der Vektor |0, 0i ist also der Hochstgewichtszustand der Darstellung D0 , der trivialen
Darstellung. In der Tat findet man nat
urlich auch, dass M |0, 0i = 0.
Wir haben also jetzt das Tensorprodukt der beiden j = 12 Darstellungen vollstandig
zerlegt
D 1 D 1 = D1 D0 ,
(9.3.10)
2
wobei D1 durch die Vektoren |1, mi mit m = 1, 0 aufgespannt wird, und D0 aus dem
Vektor |0, 0i besteht.
Man kann sich das auch geometrisch vorstellen: die beiden Spins der beiden Darstellungen sind entweder parallel (S = 1) oder antiparallel (S = 0).
S =1
1, 1
S=0
1, 0
0, 0
1,1
Abbildung 14: Singlett und Triplett Eigenzustande des Gesamtspins zweier Spin-1/2
Teil
chen: links, Spin-Singlett zum Gesamtspin S = 0, |0, 0i = [| , i | , i]/ 2, rechts,
95
9.4
Es sollte nun relativ klar sein, wie man die obigen Argumente auf den Fall zweier beliebiger
Darstellungen Dj1 Dj2 verallgemeinern kann. Das Resultat ist einfach
Dj1 Dj2 = Dj1 +j2 Dj1 +j2 1 . . . D|j1j2 | .
(9.4.1)
(9.4.2)
Wenn man die verschiedenen Vektoren nach ihrer Eigenwerte auftragt erhalt man das
Diagramm:
m2
j2
j1
j1 m1
(ji mi ji , i = 1, 2)
m1 + m2 = m
j2
m = j1 + j2
m = j1 + j2 1
..
.
: 1
: 2
..
.
f
ur
f
ur
..
.
m = |j1 j2 |
: 2 min(j1 , j2 ) + 1
m = |j1 j2 | 1 : 2 min(j1 , j2 ) + 1
..
..
.
.
ist , wobei wir angenommen haben, dass m 0. (Das Resultat ist symmetrisch unter
m m.)
Wir schliessen daraus, dass keine irreduzible Darstellung Dj mit j > j1 +j2 in Dj1 Dj2
auftritt. Die Darstellung Dj1 +j2 kommt genau einmal vor und enthalt je einen Eigenvektor
mit Eigenwert mit m = j, . . . , j. Der verbleibende Eigenvektor mit m = j1 + j2 1
generiert eine Darstellung Dj1 +j2 1 , und so weiter bis D|j1 j2 | .
96
j1 +j2
j=|j1 j2 |
(2j + 1) =
X
j=0
(2j + 1)
j1 j2 1
(2j + 1)
j=0
(j1 + j2 )(j1 + j2 + 1)
2
(j1 j2 1)(j1 j2 )
(j1 j2 ) 2
2
= (2j2 + 1) + (j1 + j2 )(j1 + j2 + 1) (j1 j2 1)(j1 j2 )
= (j1 + j2 )2 (j1 j2 )2 + 2j1 + 2j2 + 1
= 4j1 j2 + 2j1 + 2j2 + 1 = (2j1 + 1)(2j2 + 1) ,
(9.4.3)
= j1 + j2 + 1 + 2
9.5
Clebsch-Gordon Koeffizienten
Wie wir oben gesehen haben konnen wir das Tensorprodukt der beiden irreduziblen SU(2)
Darstellungen Dj1 Dj2 vollstandig nach irreduziblen SU(2) Darstellungen zerlegen. Oft
interessiert uns aber nicht nur, welche Darstellungen in dem Tensorprodukt vorkommen,
sondern wir wollen auch wissen, wie die verschiedenen Vektoren |j, mi der Darstellungen Dj durch die Vektoren des Tensorproduktes ausgedr
uckt werden konnen, d..h. wir
interessieren uns f
ur die Verallgemeinerung der Formeln (9.3.5), (9.3.6) und (9.3.7). Wir
bezeichnen die Vektoren des Tensorproduktes durch
|j1 , j2 , m1 , m2 i |j1 , m1 i |j2 , m2 i Dj1 Dj2
(9.5.1)
und die Vektoren der Darstellung Dj , die in der Clebsch-Gordon Reihe auftaucht durch
|j1 , j2 , j, mi. Offenbar konnen wir diese Basisvektoren durch die Tensorproduktbasis ausdr
ucken
X
hj , j , m , m |j , j , j, mi |j1 , j2 , m1 , m2 i .
|j1 , j2 , j, mi =
(9.5.2)
| 1 2 1 {z2 1 2
}
m ,m
1
m1 +m2 =m
C-G Koeffizient
Die dabei auftretenden Koeffizienten werden die Clebsch-Gordon Koeffizienten genannt. Wegen der Eigenwerte bez
uglich M3 folgt sofort, dass nur die Terme mit m = m1 +
m2 auftreten. Die Notation suggeriert, dass die Clebsch-Gordon Koeffizienten tatsachlich
Skalarprodukte der Vektoren sind; das ist nat
urlich richtig, da beide Basen (die Tensorproduktbasis sowie die Basis, die durch |j1 , j2 , j, mi beschrieben sind) Orthonormalbasen
(n ) sind, f
ur die gilt
X
1=
|n ihn | .
(9.5.3)
n
Clebsch-Gordon Koeffizienten sind tabelliert (siehe zum Beispiel A.R. Edmonds, Angular Momentum in Quantum Mechanics (Princeton University Press, New Jersey, 1957)
97
(9.5.6)
|m1 1, m2 i
@
@
@
@
@
@
@
@
@
@
@
@
|m1 , m2 i
|m1 + 1, m2 i
|m1 , m2 i
|m1 , m2 1i
fixed j
m2
m = m1 + m2 = j
j2
X
j1
m1
Abbildung 16: Schema zur Berechnung der Clebsch-Gordan Koeffizienten: Wir starten
bei X und ben
utzen eine Dreiecksrekursion um die Koeffizienten entlang der Kante zu
finden. Benutzung der zweiten Dreiecksrekursion liefert die Kante , usw.
folgen der Kante unter Verwendung der Dreiecksrelation (9.5.6) (beachte, dass m2 <
(j j1 )) und ben
utzen, dass f
ur m1 > j1 der Koeffizient hj1 j2 j1 +1m2 |j1 j2 jj1 +1+m2 i = 0
verschwindet; damit finden wir die Rekursion
p
j(j + 1) m(m 1)hj1 , j2 , j1 , m2 |j1 , j2 , j, m 1i
:
p
=
j2 (j2 + 1) m2 (m2 + 1)hj1 , j2 , j1 , m2 + 1|j1 , j2 , j, mi.
(9.5.7)
Anschliessend erzeugen wir mit (9.5.5) die Werte auf der Kante und arbeiten uns alsdann
iterativ durch die Matrix. Schliesslich finden wir durch die Bedingung, dass der Zustand
|j1 j2 jji normiert ist, d.h. durch
X
|hj1 j2 m1 m2 |j1 j2 jji|2 = 1 .
(9.5.8)
m1 ,m2
m1 +m2 =j
9.5.1
(9.5.9)
l m + 12
,
2l + 1
s
l m + 21
: hl, 12 , m 21 , 21 |l, 12 , l 21 , mi =
.
2l + 1
f
ur Dl+ 1 : hl, 21 , m 12 , 12 |l, 21 , l + 12 , mi =
2
f
ur Dl 1
2
99
(9.5.10)
9.6
9.6.1
Physikalische Beispiele
Magnetisches Moment
m, q
Den Ursprung dieses magnetischen Momentes kann man wie folgt verstehen. Wie aus
der Elektrodynamik bekannt beschreiben wir das externe elektro-magnetische Feld durch
seine Potentiale
~
~ = 1 A ,
~ =A
~.
E
B
(9.6.2)
c t
Die Hamiltonfunktion des geladenen Teilchens in diesem externen Feld ist dann
q ~ 2
1
~p A + q .
(9.6.3)
H=
2m
c
~ der zu ~x konjugierte Impuls. Gemass des Korrespondenzprinzips
Hierbei ist ~p = m~x + qc A
wird dieser Impuls durch
~p 7 i~
(9.6.4)
ersetzt. Damit erhalten wir den Hamiltonoperator
q2
~q
~2
2
~
~
~
A+A +
+i
A +q .
H=
2m
2mc
2mc2
(9.6.5)
~=0.
A
~ vorbeiziehen, und erhalt insgesamt
Dann kann man den Gradientoperator an A
2
~2
q ~2
~q
~
H=
+i A+
A +q .
2m
mc
2mc2
100
(9.6.6)
(9.6.7)
~ Das zugehorige
Nun betrachten wir den Fall eines reinen konstanten Magnetfeldes B.
Vektorpotenial kann man dann als
~
~ = 1 ~x B
A
2
(9.6.8)
1
= ijk j klm xl Bm
2
1
= kij kl jl Bm = Bi .
2
(~x B)
i A
mc
mc
2
~q
~
= i
(~x ) B
2mc
q ~ ~
LB ,
=
2mc
(9.6.9)
(9.6.11)
me4 1
.
2~2 n2
101
(9.6.13)
ml
1
0
1 0 1
H0
ml
ms
1/2
0
1
1
1/2
1/2
1/2
1/2
1/2
HZeemann
(9.6.14)
e~
.
2me c
(9.6.15)
Zum Beispiel spalten sich die Niveaux eines p-Zustandes (l = 1) gemass Abbildung 18
auf.
9.7
Der obige (normale) Zeemann Effekt ist jedoch nicht der einzige Effekt. Der Spin eines
Teilchens erzeugt namlich auch ein magnetisches Moment, sogar f
ur ungeladene Teilchen
(z.B. Neutronen, siehe unten). F
ur ein Elektron findet man
~
~ el = g e S,
M
2mc
(9.7.16)
wobei g der Lande Faktor ist. In der Quantenelektrodynamik kann man g berechnen und
findet
e2
1
g = 2 + + O(2 )
wobei
=
=
(9.7.17)
~c
137.04
102
~ und S
~ antiparallel, da q = e < 0.
die Feinstrukturkonstante ist. In diesem Fall sind M
F
ur ein Proton ist andererseits
~ prot = gp
M
e ~
S,
2Mp c
gp 5.59 ,
(9.7.18)
und f
ur Neutronen
~ neut = gn
M
e ~
S,
2Mn c
gn 3.83 .
(9.7.19)
qBz
(Lz + gSz ) .
2mc
(9.7.20)
Unter Ber
ucksichtigung des Spins sind die Quantenzahlen des Elektrons durch (n, l, m, ms )
gegeben, wobei ms = 12 die Spinquantenzahl ist. Mit g = 2 erhalt man also dann den
Zusatzterm zu (9.6.13) namlich (Zeemann Effekt)
EZ = B Bz (m + 2ms ) .
(9.7.21)
Dieser Effekt f
uhrt dann zu der Auspaltung (siehe Figur 18), aus der auf die Existenz
des Elektronspins geschlossen wurde (Pauli). Die Einf
uhrung des Spins war konzeptionell
ein gewaltiger Schritt: der Spin eines Quantensystems beschreibt einen Freiheitsgrad, der
kein klassisches Analogon besitzt. Insofern stosst hier unsere klassische Anschauung an
ihre Grenzen.
9.7.1
Spin-Bahn-Kopplung
Es gibt noch einen weiteren interessanten Effekt (die Spin-Bahn-Kopplung), der dazu
f
uhrt, dass wir zur Beschreibung des Problems die zuvor entwickelte Clebsch-Gordon
~
Technologie einsetzen konnen. Das sich im E-Feld
des Kernes bewegende Elektron sp
uhrt
~
~
ein Magnetfeld B = ~v E/c woran das Spin-Moment koppelt. Die resultierende Energie
ist
e ~
~ .
~v E
(9.7.22)
2S
mc
Die Thomas Prazession f
uhrt eine weitere Korrektur von derselben Form ein: diese ist ein
relativistisch kinematischer Effekt und ber
ucksichtigt, dass das Ruhesystem des Elektrons
sich relativ zum Laborsystem dreht, siehe zum Beispiel Jackson (Elektrodynamik), Seite
541 ff. Die Korrektur ergibt sich gerade zum 1/2-fachen von (9.7.22) und wir erhalten
103
folgenden Ausdruck f
ur die Spin-Bahn Kopplung:
e ~
~
S (~v E)
2mc2
~ = V
~
eE
= rr V (r)
1 ~
1
=
S (~r ~v) r V
2
2mc
r
1 1 dV ~ ~
LS .
=
2m2 c2 r dr
HSO =
(9.7.23)
~ S
~ ist nur unter der gemeinsamen Rotation von L
~ und S
~ invariant; in der
Der Term L
Tat gilt
~ S
~ = 1 (L
~ + S)
~ 2L
~2 S
~2 .
L
(9.7.24)
2
In dem Zustand mit den Quantenzahlen (n, l, m, ms ) gilt
~ 2 = ~2 l(l + 1)
L
1
3
1
2
2
~ = ~
+ 1 = ~2
S
2 2
4
2
2
~ + S)
~
(L
= ~ J(J + 1) ,
wobei J = l
1
2
(9.7.26)
haben. Da
1
e2
r V = 3
r
r
konnen wir den Beitrag zur Energie durch
e2 ~2 1
l
J =l+
ESB =
2
2
3
(l + 1) J = l
2m c 2 r nl
abschatzen. Mit
1
r3
=
nl
1 1
1
a30 n3 l l + 12 (l + 1)
104
(9.7.27)
1
2
1
2
(9.7.28)
(9.7.29)
und
e4 ~4
e2 ~2
=
= 2 e2 a20 ,
m2 c2
~2 c2 e2 m2
e4 m
e2
=
= ER = 13.6 eV,
2~2
2a0
(9.7.30)
e
~
[zur Erinnerung: a0 = me
2 ist der Bohr-Radius und = ~c die Feinstrukturkonstante]
ergibt sich die Aufspaltung
2 ER
1
l
J = l + 21
ESB =
.
(9.7.31)
2n3 l 1 + 12 (l + 1) (l + 1) J = l 12
F
ur kleine externe Magnetfelder ist der Spin-Bahn-Kopplung Effekt relativ zum ZeemanEffekt dominant.
105
10
Wie wir gesehen haben macht die Quantenmechanik nur Aussagen dar
uber, mit welcher
Wahrscheinlichkeit ein bestimmtes Messergebnis eintreffen wird. In der klassischen Physik
gibt es nat
urlich auch Situationen, wo man effektiv nur Wahrscheinlichkeitsaussagen machen kann: zum Beispiel kann man die Lotto-Zahlen nicht wirklich vorhersagen, sondern
lediglich (hoffentlich) die Aussage machen, dass alle Zahlen mit derselben Wahrscheinlichkeit gezogen werden werden. Nat
urlich interpretiert man das nicht dahingegen, dass die
klassische Physik nicht deterministisch ware: der Grund daf
ur, dass man die Lottozahlen nicht vorhersagen kann, liegt einfach darin, dass man die Anfangsbedingungen nicht
genau genug kennt.
Man konnte also denken, dass der Grund daf
ur, dass wir in der Quantenmechanik nur
Wahrscheinlichkeitsaussagen machen konnen ahnlicher Natur ist. Insbesondere konnte
man sich vorstellen, dass der wirkliche Zustand eines (quantenmechanischen) Systems
tatsachlich alle Messwerte eindeutig festlegt, dass wir aber den wirklichen Zustand nicht
genau genug kennen, und deshalb nur Wahrscheinlichkeitsaussagen machen konnen. Diese
Idee wird u
blicherweise als Theorie verborgener Variablen (hidden variables) bezeich
net. Sie ware auch aus abstrakten (philosophischen) Uberlegungen
sehr attraktiv, da eine
solche Theorie die physikalische Wirklichkeit direkt abbilden w
urde. [Falls die Quantenmechanik keine Theorie verborgener Variablen ist, scheint sie zu suggerieren, dass es keine
physikalische Realitat gibt: der Wert bestimmter Variablen ist dann (bevor sie gemessen
werden) nicht eindeutig bestimmt!]
Wir wir im folgenden zeigen wollen, scheint die Quantentheorie jedoch keine lokale
Theorie verborgener Variablen zu sein.
10.1
Der Unterschied zwischen der Quantenmechanik und einer Theorie verborgener Variablen
kann im Kontext des sogenannten EPR Paradoxes verstanden werden. Dieses Paradox
geht auf Einstein, Podolsky und Rosen (1935) zur
uck. Dazu betrachten wir zwei Spin 1/2
Teilchen, die sich in dem Spin-0 Zustand des Gesamtsystems befinden
1
|0, 0i = (| i | , i) .
2
(10.1.1)
Dieser Zustand ist invariant unter gemeinsamen Rotationen der beiden Spin-Freiheitsgrade; er ist deshalb von der urspr
unglichen Wahl der 3-Richtung (bez
uglich derer | i und
| i definiert sind) unabhangig. Man nennt einen solchen Zustand oft ein EPR-Paar.
Er beschreibt zum Beispiel eine der Zerfallsmoglichkeiten des Pions, e e+ , in ein
Elektron und ein Positron.
Die Teilchen laufen danach frei auseinander, was ihren Spinzustand (10.1.1) nicht
(i)
andert. Je eine Spinkomponente Sj , (j = 1, 2, 3), der beiden Teilchen, (i = 1, 2), wird
dann gemessen. Die beiden Messungen seien raumartig getrennte Ereignisse, so dass nach
106
vorher:
nachher:
(1)
Sj 1
(2)
Sj 2
|i
der speziellen Relativitatstheorie eine kausale Beinflussung ausgeschlossen ist. Die Mes(1)
(2)
sung von Sj und Sj (gleiche Richtung j) ergibt entweder (~/2, ~/2) oder (~/2, ~/2),
je mit Wahrscheinlichkeit 1/2. (Diese Wahrscheinlichkeit ist unabhangig davon, in welcher
(1)
Richtung der Spin gemessen wird.) Also kann S1 vorausgesagt werden, ohne das Teilchen
(2)
(1) zu storen, namlich durch Messung von S1 . Dies klingt paradox, da die beiden Messungen beliebig weit entfernt voneinander stattfinden konnen. Die dieser Beobachtung zu
Grunde liegende Idee ist das Prinzip auf dem die moderne Quanten Teleportation beruht
(siehe Kapitel 10.3).
10.2
Die Bellsche Ungleichung bezieht sich nicht auf die Quantenmechanik, sondern gilt f
ur
beliebige lokale Theorien verborgener Variablen. In diesem Rahmen kann eine Ungleichung
abgeleitet werden (Bellsche Ungleichung), die in der Quantenmechanik nicht erf
ullt ist.
In einer Theorie verborgener Variablen gibt es einen Wahrscheinlichkeitsraum , auf
dem die verborgenen Variablen leben. Da wir diese Variablen nicht kennen (bzw.
nicht direkt messen konnen) beschrankt sich unsere Kenntnis des Zustands des System
auf ein Wahrscheinlichkeitsmass d auf . Der Erwartungswert einer Observable A ist
dann einfach
Z
hAi =
A()d() .
(10.2.1)
Im weiteren wollen wir annehmen, dass die Theorie lokal ist: falls A und B zwei raumartig
getrennte Observable sind, dann soll gelten
(AB)() = A() B() .
(10.2.2)
Hierbei bedeutet (AB) die Oberservable, bei der zuerst B und dann A (raumartig getrennt) gemessen wird. Insbesondere bedeutet diese Annahme, dass die beiden Messungen
A und B miteinander vertauschen. Die physikalische Motivation hinter (10.2.2) ist dass
in der mikroskopischen Beschreibung (durch alle relevanten Variable, inklusive jener auf
) sich raumartig getrennte Messungen nicht gegenseitig beeinflussen konnen.
Unter diesen Annahmen kann man nun die Bellsche Ungleichung [Bell (1964),
Clauser, Horne, Shimony, Holt (1969)] ableiten: Seien A, A von B wie von B raumartig
getrennte Observable, die alle nur die Werte 1 annehmen konnen. Dann gilt
|hABi + hA Bi + hAB i hA B i| 2 .
107
(10.2.3)
(10.2.4)
In der Tat gilt dass der Ausdruck (10.2.4) gerade = 2 ist, denn
A() = A ()
A() = A ()
A() + A () = 2 , A() A () = 0 ,
A() + A () = 0 , A() A () = 2 .
=
=
(1)
A = ~n = ~n1 ~ (1) ,
1
(2)
(2)
B = ~n = ~n2 ~ (2) ,
2
wobei sich (1) und (2) auf die beiden Teilchen bezieht. Die Messung von AB zum Beispiel
entspricht dann der Spin-Spin Korrelation
(1)
(2)
(10.2.6)
Falls das System durch eine lokale Theorie verborgener Variablen beschrieben wird, so
folgt nun aus (10.2.3), dass
|C(~n1 , ~n2 ) + C(~n1 , ~n2 ) + C(~n1 , ~n2 ) C(~n1 , ~n2 )| 2 .
(10.2.7)
Diese Vorhersage soll nun mit der QM verglichen werden. Dazu betrachten wir den Zustand |i, der das EPR Paar (10.1.1) beschreibt. F
ur jeden Vektor ~v gilt dann (wir konnen
die 3-Achse bez
uglich derer das System geschrieben ist, jeweils in Richtung des Vektors ~v
wahlen)
(~ (1) + ~ (2) ) ~v |i = 0
h|~ (1) ~v|i = 0 .
(10.2.8)
(10.2.9)
Damit ist
C(~n1 , ~n2 ) = h|(~(1) ~n1 )(~(2) ~n2 )|i = h| (~(1) ~n1 )(~(1) ~n2 ) |i = ~n1 ~n2 , (10.2.10)
{z
}
|
~
n1 ~
n2 1I+i~
(1) (~
n1 ~
n2 )
108
(10.2.11)
~n1
~n2
~n2
/4 /4
/4
~n1
Diese Analyse zeigt, dass die Quantenmechanik (so wie wir sie formuliert haben) nicht
zu einer lokalen Theorie verborgener Variablen aquivalent ist! Dies ist auch experimentell
bestatigt worden. Die Experimente [Freedman-Clauser (1972); Aspect, Dalibard & Roger (1982)] ben
utzen dabei polarisierte Photonen statt Spin 12 -Teilchen. Die gefundenen
Ergebnisse stimmen mit den Vorhersagen der Quantenmechanik sehr gut u
berein.
10.3
Quanten Teleportation
Wie wir schon in Kapitel 10.2 erwahnt haben, impliziert das EPR Paradox, dass es eine
Art Fernwechselwirkung zwischen den Spinfreiheitsgraden gibt. Diese Wechselwirkung
kann dazu ben
utzt werden, Quanteninformation zu teleportieren. Dazu betrachten wir die
folgende Anordnung.
Alice und Bob besitzen je ein Spin 12 -Teilchen eines EPR Paars (Teilchen 1 und 2)
1
|i = (| i | i) .
2
(10.3.1)
Alice soll den unbekannten Zustand |i eines weiteren Spin-1/2 Teilchens (Teilchen 0) an
Bob u
bermitteln, und zwar unter Verwendung bloss klassischer Information, d.h. endlich
vieler Bits. Insbesondere ist es deshalb nicht zugelassen, einfach das Teilchen von Alice
zu Bob zu bewegen. Naiverweise scheint das unmoglich zu sein, da Alice nicht direkt
|i messen kann, sondern lediglich den Spin in einer vorgegebenen Richtung. Dabei geht
jedoch die genaue Information u
ber den Vektor |i verloren (Kollaps der Wellenfunktion).
Mit Hilfe von (10.3.1) konnen Alice und Bob den Zustand des Teilchens 0 auf das
Teilchen 2 dennoch u
bertragen (Teleportation). Dazu betrachtet man den Zustand aller
drei Teilchen
1
(10.3.2)
|i |i = (|, i | i |, i | i) .
2
109
Alice kann Messungen an ihren Teilchen 0 und 1 vornehmen, wie z.B. die mit Zerlegung
1I =
4
X
i=1
Pi
4
X
i=1
Pi 1I ,
(10.3.3)
,
,
,
,
Alice u
bermittelt das Ergebnis i = 1, 2, 3, 4 (2 Bits) an Bob. Je nach Ergebnis u
bt er
folgende unitare Operatoren auf sein Teilchen 2 (realisierbar durch Spinprazession) und
erhalt dessen Zustand:
Alices Ergebnis Bobs Operation Zustand
1
2
3
4
Hier haben wir ausgen
utzt, dass
1I2
3
1
2
3 | i = | i ,
1 | i = | i ,
2 | i = i| i ,
|i
|i
|i
i|i
3 | i = | i ,
1 | i = | i ,
2 | i = i| i .
In allen Fallen ist der Zustand |i (Phase ohne Bedeutung) wiederhergestellt! Quanten
Teleportation ist mit Photonen experimentell realisiert worden [Zeilinger et al. (1997)].
110
11
Sto
rungstheorie
Die meisten interessanten physikalischen Systeme lassen sich nicht exakt losen. Man kann
sie jedoch haufig als Storung eines losbaren Systems auffassen. Dann kann man das Problem in Storungstheorie behandeln. Wir wollen zunachst den Fall betrachten, wo der
Hamiltonoperator zeit-unabhangig ist.
11.1
Nicht-entartete zeit-unabh
angige Sto
rungstheorie
(11.1.1)
H = H0 + H
(11.1.2)
schreiben lasst, wobei H0 exakt losbar ist. Wir bezeichnen den vollstandigen Satz von
orthonormierten Eigenfunktionen von H0 durch n , n N, wobei
H 0 n = E n n
(11.1.3)
(11.1.4)
(11.1.5)
:
:
:
:
(H0 E0 )0 =
(H0 E0 )1 =
(H0 E0 )2 =
(H0 E0 )3 =
hochstes n
0,
(E1 H )0 ,
(E1 H )1 + E2 0 ,
(E1 H )2 + E2 1 + E3 0 .
hochstes En
111
(11.1.6)
E0 = En
(11.1.7)
f
ur ein geeignetes n N. Die weiteren k f
ur k 1 konnen wir nun iterativ finden,
da auf der rechten Seite der obigen Gleichungen jeweils nur l mit l < k auftreten.
Allerdings bestimmen diese Gleichungen k nur bis auf Vektoren, die von (H0 E0 )
vernichtet werden. Da wir angenommen haben, dass die verschiedenen Eigenwerte von H0
nicht-entartet sind, ist jeder solche Vektor zu n proportional. Ein entsprechender Term
w
urde also lediglich die Normierung des ersten Terms 0 = n verandern. Wir legen diese
Freiheit dadurch fest, dass wir verlangen, dass
h0 |k i = 0k .
(11.1.8)
(11.1.9)
[Diese Normierungsbedingung f
uhrt dann aber dazu, dass nicht richtig normiert sein
wird; am Ende der Storungsrechnung werden wir also das resultierende noch richtig
normieren m
ussen.]
Nach diesen Vorbemerkungen wollen wir nun (11.1.6) systematisch losen. Das Skalarprodukt mit 0 der Ordnung s -Gleichung gibt
s1 E1 h0 |H |s1 i
Es = hn |H |s1 i ,
s
X
j=2
Ej h0 |sj i
| {z }
s1.
(11.1.10)
sj
(11.1.11)
Die Berechnung von E in Ordnung s setzt also die Berechnung von in der Ordnung
s 1 voraus. Die Multiplikation von (11.1.11) mit s und Addition der Terme s ergibt
E = E0 +
s Es
s=1
= E0 +
s=1
s1 hn |H |s1i
= En + hn |H |i .
Um s zu finden entwickeln wir nun
X
|s i =
hl |s i |l i ,
l6=n
112
(11.1.12)
s1.
(11.1.13)
s
X
j=1
Ej hl |sj i hl |H |s1 i .
(11.1.14)
s
X
1
Ej hl |sj i ,
hl |s i =
hl |H |s1 i
En El
j=1
(11.1.15)
Es = hn |H |s1i
Ps1
X hl |H |s1 i j=1
Ej hl |sj i
|s i =
|l i ,
En El
l6=n
(11.1.16)
wobei wir in der ersten Gleichung (11.1.11), und in der zweiten Gleichung (11.1.13) und
(11.1.15) ben
utzt haben. In der letzten Summe haben wir den Term j = s weggelassen,
da hl |0 i = 0 nach Annahme.
F
ur s = 0 bilden E0 = En und |0 i = |n i den Startpunkt f
ur die Iteration. Wir
beobachten, dass (11.1.16) nur f
ur den nicht-entarteten Fall Sinn machen, da ansonsten
der Nenner der zweiten Gleichung nicht ungleich Null ist; wir werden auf den Spezialfall,
wo diese Bedingung nicht erf
ullt ist, im nachsten Kapitel zur
uckkommen.
Schliesslich geben wir die Resultate f
ur s 2 explizit an, wobei wir die Notation
|l i = |li verwenden:
0 :
1 :
E0 = En ,
|0 i = |ni .
E1 = hn|H |ni ,
X hl|H |ni
|li .
|1 i =
E
n El
l6=n
E2 =
2 :
|2 i =
X |hl|H |ni|2
,
En El
l6=n
XhX hl|H |kihk|H |ni
l6=n k6=n
(En El )(En Ek )
(En El )2
11.1.1
Gest
orter harmonischer Oszillator
(11.1.18)
mit
m
p2
+ 2 q2 ,
H = F q .
2m
2
Wie wir in Kapitel 6 gesehen haben, hat das ungestorte Problem die Losung
1
H0 |ni = En |ni ,
E n = ~ n +
.
2
H0 =
(11.1.19)
(11.1.20)
a |li = l + 1 |l + 1i
(11.1.22)
a|li = l |l 1i ,
folgt dann f
ur die Matrixelemente
r
~(n + 1)
hl|n 1i F
hl|n + 1i
2m
2m
r
r
~n
~(n + 1)
= F
l,n1 F
l,n+1 .
2m
2m
hl|H |ni = F
~n
Damit erhalten wir die folgenden Resultate bis zur 2-ten Ordnung
1
te
~ ,
0 : E0 = n +
2
|0 i = |ni .
te
1 : E1 = hn|H |ni = 0 ,
X hl|H |ni
|1 i =
|li
E
n El
l6=n
hn 1|H |ni
hn + 1|H |ni
|n 1i +
|n + 1i
~
~
r
F
~
n|n 1i n + 1|n + 1i
=
~
2m
r
m~ 1
iF
p|ni mit p =
(a a ) .
=
2
m~
2 i
(11.1.23)
(11.1.24)
114
(11.1.25)
2te :
E2
X |hl|H |ni|2
l6=n
En El
F2
~n
~(n + 1)
=
= F2
,
2m(~) 2m(~)
2m 2
2
F
m~
1
2
p |ni
(2n + 1)|ni . .
|2 i =
2 m~ 2
2
=
(11.1.26)
[Der Korrekturterm in der letzten Zeile kommt daher, dass nach Voraussetzung |2 i
orthogonal zu |ni ist der Beitrag von p2 |ni zu |ni muss daher abgezogen werden.]
Andererseits beschreibt der Hamiltonian
p2
1
+ m 2 q 2 F q
2m 2
F 2
m 2
F2
p2
q
+
=
2
2
2m
2
|m
{z } | 2m
{z }
H =
q0
E2
|ni
|ni .
= |ni
m~ 2
2 m~ 2
(11.1.27)
Wir sehen, dass die Storungsreihe gerade den verschobenen harmonischen Oszillator als
Potenzreihe in F ergibt der Unterschied in dem F 2 -Term proportional zu |ni folgt
daher, dass |0 i und das durch die Storungsreihe definierte |i unterschiedlich normiert
sind: |0 i ist kanonisch normiert, h0 |0 i0 = 1 wohingegen die Norm von |i durch
(11.1.17) gegeben ist, und sich zur Ordnung F 2 von 1 unterscheidet.
11.2
Entartete zeit-unabh
angige St
orungstheorie
Bisher haben wir angenommen, dass alle Eigenwerte des ungestorten Hamiltonoperators
nicht-entartet waren. Diese Annahme ging explizit in die obige Analyse ein, da wir zum
Beispiel in (11.1.15) durch En El geteilt haben, was nur Sinn macht, falls En 6= El f
ur
alle l 6= n. Falls jedoch der Eigenwert E0 = En , um den wir storen entartet ist, bricht
diese Analyse zusammen. Es ist jedoch nicht schwierig zu verstehen, wie dieses Problem
umgangen werden kann.
11.2.1
Bevor wir den allgemeinen Fall diskutieren wollen ist es instruktiv, den Fall zu betrachten,
bei dem der relevante Eigenwert zweifach entartet ist. Sei also E0 = En = Em , wobei alle
115
anderen Eigenwerte Ek 6= En f
ur k 6= n, m sind. Wir nehmen weiterhin an, dass l wie
zuvor eine Orthonormalbasis der zugehorigen Eigenvektoren beschreibt.
Es ist klar, dass die Formeln (11.1.6) weiterhin gelten. Wir konnen diese Gleichungen
aber nun nicht mehr so einfach invertieren, da H0 E0 nun einen zwei-dimensionalen
Vektorraum von Vektoren vernichtet. Insbesondere ist 0 nun nicht mehr (bis auf Normierung) eindeutig festgelegt; statt dessen wissen wir lediglich, dass
0 = am |mi + an |ni ,
(11.2.1)
(11.2.2)
und entsprechend f
ur das Matrixelement mit hn|. Die beiden Gleichungen lassen sich in
die folgende Matrix-Form bringen,
hm|H |mi E1 am + hm|H |ni an = 0
hn|H |mi am + hn|H |ni E1 an = 0 .
(11.2.3)
Wir definieren nun
Em
= hm|H |mi
En = hn|H |ni
= hm|H |ni ,
= hn|H |mi .
=0.
an
En E1
(11.2.4)
(11.2.5)
Eine Losung 0 6= 0 existiert falls die Determinante gleich Null ist, also falls
(Em
E1 )(En E1 ) | |2 = 0 ,
(11.2.6)
Em
+
2
En
s
E
Em
n
2
2
+ | |2 .
(11.2.7)
m
=
a
E1 Em
n
116
(11.2.8)
wobei a
achste Ordnung ergibt die Korrektor
m , an durch (11.2.8) bestimmt ist. Die n
E1 = E1 gemass (11.2.7). Falls Em 6= En oder 6= 0 ist E1+ 6= E1 , und die Entartung ist
aufgehoben; das ist der generische Fall.
Um die Korrektur |
ur den Eigenvektor von H zu finden betrachten wir zuerst
1 i f
die Ordnung s = 1 von (11.1.6) und berechnen das Matrixelement mit hl|, l 6= m, n,
hl|H0 E0 |
1 i = hl|E1 H |0 i
(El E)hl|
1 i = E1 hl|0 i hl|H |0 i ,
| {z }
(11.2.9)
wobei wir E = Em = En geschrieben haben. Diese Gleichung legt also die Komponenten
von |
ur die Vektoren im orthogonalen Komplement von span(|mi, |ni) fest,
1 i f
|P
1i =
X hl|H | i
0
|li ,
E El
(11.2.10)
l6=m,n
wobei P = 1I |mihm| |nihn| der Projektor auf diesen Unterraum ist. Um die Komponenten von 1 im Unterraum span(|mi, |ni) zu bestimmen, wahlen wir zunachst die
m
ussen wir daher nur h1 |0 i bestimmen. Dazu ben
utzen wir die Ordnung 2 Gleichung
von (11.1.6). F
ur |+
1 i finden wir
+
+
+
+
+
h
0 | (H0 E0 )|2 i = (E1 H )|1 i + E2 |0 i
+
+
+
+
h
|H E |+ i = h
0 |E1 H |1 i + E2 h0 |0 i ,
| {z }
| 0 0 {z 0 2 }
0
woraus
h+
1 |E1 H |0 i = 0
(11.2.11)
(E1+ H )|
0 i = (E1 P H )|0 i (1l P )H |0 i
= (E1+ E1 )|
0 i P H |0 i .
(11.2.12)
(11.2.13)
da (1I P )H |
0 i = E1 |0 i wie man durch direktes Nachrechnen zeigt. [Dies ist gerade
die Eigenwertbedingung f
ur
utzung von (11.2.12) nehmen wir jetzt das
0 .] Unter Ben
+
innere Produkt mit h1 | und finden dann (nach Konjugation)
+
+
0 = (E1+ E1 )h
0 |1 i h0 |H P |1 i
(11.2.14)
+
und damit, unter Verwendung von (11.2.10), die gesuchte Komponente h
0 |1 i. Wir
beobachten, dass die
+
h
0 |1 i
X h |H |lihl|H |+ i
0
0
=
+
O(H )
(E El )(E1 E1 )
l6=m,n
117
(11.2.15)
+
+
|0 i = am |mi + an |ni ;
|0 i = a
m |mi + an |ni ;
= E1+
E1
i
X h h |H |lihl|H |+ i
hl|H |+
0
0
0i
+
|li ;
| i =
| i +
1
E El
(E1+ E1 )(E E l ) 0
l6=m,n
= E1
E1
i
X h h+ |H |lihl|H | i
hl|H |
0
0
0i
+
(11.2.16)
|
i
=
|li
.
|
i
+
0
+
1
E
E
(E
E
)(E
E
)
l
l
1
1
l6=m,n
11.2.2
Im allgemeinen Fall, d.h. falls der Eigenraum zu En k-dimensional ist, bezeichnen wir mit
{1n , , kn } eine Orthonormalbasis dieses Eigenraumes. F
ur 0 machen wir nun den
Ansatz
k
X
0 =
ajn jn .
(11.2.17)
j=1
h1n |H |kn i
2
1
2
2
2
k
a2n
h
|H
|
i
(h
|H
|
i
E
)
h
|H
|
i
1
n
n
n
n
n
n
. = 0 .
..
..
..
.
.
k
1
k
k
akn
hn |H |n i
(hn |H |n i E1 )
(11.2.18)
Wir bezeichnen die Eigenwerte mit E1 = hn |H |n i, = 1, . . . , k, wobei die zugehorigen
Eigenvektoren durch
k
X
|n i =
aln () |ln i
(11.2.19)
l=1
gegeben sind. Die weitere Rechnung funktioniert analog wie das, was wir oben besprochen
haben. Zum Beispiel findet man zu zweiter Ordnung f
ur E
E = En + hn |H |n i +
118
X |hl|H |n i|2
,
En El
E 6=E
l
(11.2.20)
| i = |n i +
6=
E1
X hn |H |lihl|H |n i
1
|n i
En El
E1 E 6=E
l
X hl|H |n i
+
|li .
E
n El
E 6=E
l
(11.2.21)
In der zweiten Gleichung haben wir angenommen, dass die Entartung bereits zu erster
ur 6= .
Ordnung aufgehoben wurde, d.h. dass E1 6= E1 f
11.2.3
Im allgemeinen kann es jedoch passieren, dass die Entartung zu erster Ordnung noch nicht
aufgehoben wird. Zum Beispiel wird die Entartung f
ur den Fall der zwei-fachen Entartung
dann nicht aufgehoben, falls
hm|H |mi = hn|H |ni ,
hm|H |ni = 0 ,
(11.2.22)
siehe (11.2.7) mit (11.2.4). In diesem Fall muss man die Rechnung wiederum entsprechend anpassen. Dazu betrachten wir die s = 2 Ordnung von (11.1.6) und nehmen das
Matrixelement mit hm| und mit hn|,
hm|H0 E0 |2i = hm|E1 H |1 i + E2 hm|0 i .
(11.2.23)
Ferner ben
utzen wir die 1-te Ordnung Resultate
|0 i = am |mi + an |ni ,
E0 = Em = En = E
hn|H |mi = E nm ,
X am hl|H |mi + an hl|H |ni
|li ,
P |1i =
E
E
l
l6=m,n
(11.2.24)
(11.2.25)
wobei die letzte Gleichung gerade (11.2.10) ist. Die 2-te Ordnung Gleichung (11.2.23) lasst
sich nun als
0 = hm|E1 H |P 1 i + hm|E1 H |(1l P )1 i +E2 hm|0 i
|
{z
}
(11.2.26)
119
Em
(11.2.28)
(En E2 )an + am = 0 ,
(11.2.29)
wobei En analog zu Em
definiert ist,
En =
X |hl|H |ni|2
.
E El
l6=m,n
(11.2.30)
Die Gleichungen (11.2.3) und (11.2.29) haben genau dieselbe Form wie (11.2.5). Die Entartung wird also zu zweiter Ordnung aufgehoben, falls
entweder
11.2.4
En 6= Em
oder 6= 0 .
(11.2.31)
Als Anwendung betrachten wir das Wasserstoffatom im elektrischen Feld E~ = (0, 0, E).
Zustatzlich zu dem Hamiltonoperator des Wasserstoffatoms haben wir also jetzt die
Storung H = eEz. Wir untersuchen die Aufhebung der Entartung im n = 2 Niveau
mit den Zustanden |n, l, mi = |2s0 i, |2p1 i, |2p0i, |2p1 i. [Wir haben 4 = n2 entartete
Zustande; hier bezeichnet |2s0 i den Zustand mit l = m = 0, und |2pm i die Zustande mit
l = 1.] Im Prinzip m
ussten wir also eine 4 4 Matrix betrachten, da aber [Lz , z] = 0
vertauscht H mit Lz und die Matrixelemente von H zwischen Zustanden mit unterschiedlichem Lz Eigenwert (= m) verschwinden
0 = hm |[H , Lz ]|mi = (m m )hm |H |mi .
(11.2.32)
Somit bleiben nur die Matrixelemente h2s0 |H |2p0 i = , ihr komplex konjugiertes ,
sowie die diagonalen Matrixelement. Die Kugelfunktionen sind Eigenfunktionen des Paritatoperators P mit Paritat (1)l . Andererseits hat der Operator z Partitat 1, und es
folgt
(1)l hYlm |H |Ylmi = hYlm |H P |Ylm i = hYlm |P H |Ylmi = (1)l hYlm |H |Ylm i ,
(11.2.33)
und daher verschwinden all diagonalen Matrixelemente von H . Die Matrixelemente von
H auf dem entarteten Unterraum sind also durch die Matrix beschrieben
0 0 0
s
p1
0 0 0 0 .
(11.2.34)
H :
0 0 0
p0
0 0 0 0
p1
120
Schliesslich berechnet sich zu = h2s0 |H |2p0 i = 3a0 e E, wobei a0 der Bohrradius ist.
Es bleibt die 2 2 Matrix zu diagonalisieren
1
0 1
1
3a0 e E
= EW = 3aB e E , EV =
.
(11.2.35)
1 0
2 1
Die Aufhebung des 4-fach entarteten Zustandes bei e2 /8a0 (n = 2) in zwei nicht-entartete
und einen 2-fach entarteten Zustand ist in Abb. 19 skizziert.
H0
H0 + H
( 2 s 0 2p 0 )/ 2
deg = 4
11.3
deg = 2
2p+1 , 2p1
( 2 s 0 + 2p 0 )/ 2
Zeit-abh
angige St
orungstheorie
Bisher haben wir angenommen, dass sowohl der ungestorte Hamiltonoperator H0 , wie
auch die Storung H zeit-unabhangig sind. Wir wollen nun den Fall betrachten, bei dem
H (t) von der Zeit abhangt; wir wollen jedoch weiter annehmen, dass der ungestorte
Hamiltonoperator zeit-unabhangig ist. Um dieses Problem zu losen sollten wir ein wenig
weiter ausholen und zunachst das sogenannte Heisenberg- und das Wechselwirkungsbild
besprechen.
11.3.1
(11.3.1)
121
(11.3.2)
Es ist dann leicht einzusehen, dass U(t, s) die Zeitentwicklung generiert: sei der Zustandsvektor zur Zeit s als (s) gegeben, dann ist der Zustandsvektor zur Zeit t gerade
(t) = U(t, s)(s) .
(11.3.3)
Falls H nicht von der Zeit abhangt (das ist das Analogon der autonomen Systeme in
der Mechanik), kann man U(t, s) direkt durch
(t s)
(11.3.4)
U(t, s) = exp iH
~
X
U(t, s) = 1I +
U (n) (t, s) ,
(11.3.6)
n=1
wobei
(n)
n Z t Z t1
Z tn1
1
(t, s)=
dt1
dt2 . . .
dtn H(t1 )H(t2 ) . . . H(tn ) .
i~
s
s
s
spatere
fr
uhere
(11.3.7)
Zeiten.
Diese Zeitordnung ist unangenehm, aber formal behebbar: die Zeitordnung impliziert eine
Integration u
ur n = 2
ber den Teil t t1 tn s des Hyperkubus im Rn , vgl. Abb. 20 f
Als formales Werkzeug definieren wir den Zeitordnungs-Operator T , der eine Sequenz von
zeitabhangigen Operatoren zeitlich ordnet,
A(t1 ) B(t2 ) t1 > t2
T A(t1 ) B(t2 ) =
(11.3.8)
B(t2 ) A(t1 ) t2 > t1 .
122
t2
t
Abbildung
20:
Integrationsbereich
(grau
schraffiert) u
ber den zeitgeordneten Sektor f
ur den Fall
n = 2.
t0
t1
(11.3.9)
i=1
dt1
t1
dt2
tn1
n Z
1 Y t
dti T .
dtn =
n! i=1 t0
(11.3.11)
(n)
1
(t, s) =
n!
n Z t
Z t
1
dt1 . . .
dtn T (H(t1 ) H(tn ))
i~
s
s
(11.3.12)
bzw. f
ur den Propagator
U(t, s) = T exp
i
~
dt H(t )
(11.3.13)
Wir bemerken, dass sich diese Formel vereinfacht, falls [H(t1 ), H(t2)] = 0 f
ur alle t1 , t2 ,
da wir dann die Zeitordnung T weglassen konnen.
123
11.3.2
Das Heisenberg-Bild
Nach diesen Vorbemerkungen konnen wir nun das Heisenberg-Bild erklaren. Im Heisenberg-Bild walzt man die Zeitevolution von den Zustanden auf die Operatoren ab. Diese
beiden Beschreibungen sind durch den unitaren Propagator U(t, s) miteinander verbunden
und daher aquivalent. Konkret definieren wir die Zustande im Heisenberg-Bild durch
H (t0 ) ,
(11.3.14)
d.h. der Zustand ist einfach der Zustand im Schrodingerbild zu einer festen Zeit. Anders
ausgedr
uckt bedeutet dies, dass
H = U(t0 , t) (t) .
(11.3.15)
Die beiden Beschreibungen hangen also mit dem (zeit-abhangigen) unitaren Operator
U(t0 , t) miteinander zusammen. Entsprechend werden dann auch die Observablen der
urspr
unglichen Beschreibung (d.h. des Schrodinger-Bildes) umtransformiert:
AH = U(t0 , t) A U(t, t0 ) .
(11.3.16)
(11.3.17)
Insbesondere sind daher auch die Erwarungswerte in den beiden Beschreibungen gleich
hAi = h(t)|A|(t)i
= h(t0 )|U (t, t0 ) A U(t, t0 ) |(t0 )i
= hH |AH (t)|H i ,
(11.3.18)
(11.3.19)
d
AH (t) = (i~t U(t, t0 ) ) A U(t, t0 ) + U(t, t0 ) A i~t U(t, t0 ) + U(t, t0 ) i~t A U(t, t0 )
| {z } |
{z
}
{z
}
|
dt
(HU (t,t0 ))
H U (t,t0 )
124
i~t AH
(11.3.20)
Falls t H = 0 so ist [H, U(t, t0 )] = 0 und damit HH (t) = HH = H. Dann ist also HH
zeit-unabhangig, wahrend f
ur die anderen Operatoren die Zeitabhangigkeit durch
AH (t) = eiH(tt0 )/~ A eiH(tt0 )/~
(11.3.21)
[H, A] = 0 .
(11.3.22)
Dies ist das quantenmechanische Analogon der Aussage, dass eine Funktion des Phasenraumes und der Zeit A(p, q, t) in der klassischen Physik dann erhalten ist, falls t A = 0 und
die Poissonklammer {A, H} = 0 verschwindet. In der Tat stimmt die Bewegungsgleichung
im Heisenbergbild gerade mit der klassischen Formel
dA
= {A, H} + t A
dt
(11.3.23)
u
berein siehe Mechanikskript Kapitel 6.3. Diese Relation ist kein Zufall: eine Art die
Quantisierung ein wenig formaler zu verstehen besteht darin, dass man die Poissonalgebra
(die die Algebra der Observablen eines klassischen mechanischen Systems beschreibt)
quantisiert. Dabei wird die Poissonklammer durch den Kommutator ersetzt (Dirac):
i
{A, B} = [A, B] .
~
(11.3.24)
Zum Beispiel kann man in dieser Weise auch die kanonischen Vertauschungsregeln zwischen Ort und Impuls erhalten, denn
i
{q, p} = 1 = [q, p] = 1 ,
~
(11.3.25)
Das Wechselwirkungsbild
F
ur die Beschreibung der zeit-abhangigen Storungstheorie ist weder das Schrodinger- noch
das Heisenberg-Bild, sondern das sogenannte Wechselwirkungsbild (oder Dirac-Bild) am
angemessensten. Wir betrachten also den Fall, wo der Hamiltonoperator von der Form
H = H0 + H (t)
(11.3.26)
ist, wobei H0 zeit-unabhangig und exakt losbar ist, und H (t) eine kleine Storung beschreibt. Die Idee des Wechselwirkungsbildes ist, dass wir die (trivale) Zeitentwicklung,
die von H0 herkommt von den Zustanden auf die Observablen hin
uberschieben; das Wechselwirkungsbild ist also ein Kompromiss, wo sowohl die Zustande als auch die Observablen
zeit-abhangig sind.
125
Konkret betrachten wir also die unitare Transformation mit U0 (t, t0 ), wobei U0 der
Propagator des ungestorten Problems ist, also (11.3.2) mit H0 erf
ullt. Dann definieren
wir
D (t) = U0 (t0 , t)(t) = U0 (t0 , t) U(t, t0 ) (t0 ) ,
(11.3.27)
wobei U der Propagator f
ur H = H0 + H ist, und entsprechend
AD (t) = U0 (t0 , t)A U0 (t, t0 ) .
(11.3.28)
(11.3.29)
(11.3.30)
beschrieben, wobei
ist. Der zugehorige Hamiltonoperator ist
i~t UD = i~t eiH0 (tt0 )/~ U(t, t0 )
= U0 (t0 , t) (H0 + H) U0 (t, t0 ) U0 (t0 , t)U(t, t0 )
{z
}
{z
}
|
|
H
UD
HD
UD
(11.3.31)
mit
HD
= U0 (t0 , t) H U0 (t, t0 ) = eiH0 (tt0 )/~ H eiH0 (tt0 )/~ .
(11.3.32)
UD (t, t0 ) = T exp
11.3.4
i
~
t0
dt
HD
(t )
(11.3.33)
Erste Ordnung St
orungstheorie
Wir sind nun in der Lage, diese Ideen auf das Problem der zeit-abhangigen Storungstheorie
anzuwenden. Wie zuvor betrachten wir ein System, bei dem der Hamiltonoperator von
der Form (11.3.26) ist. Wir nehmen an, dass f
ur t < t0 H = 0 ist, d.h. dass die Storung
erst zur Zeit t = t0 eingeschaltet wird. Die Zeitentwicklung f
ur t < t0 ist also durch H0
beschrieben. Weiterhin konnen wir ohne Beschrankung der Allgemeinheit annehmen, dass
das System zur Zeit t = t0 in dem Eigenzustand |ii (initial) von H0 mit Eigenwert Ei
prapariert ist. Wir wollen die Wahrscheinlichkeit berechnen, mit der das System zu einem
spateren Zeitpunkt t > t0 im Zustand |f i (final) ist, wobei wir annehmen wollen, dass
hf |ii = 0.
Wie wir im vorigen Unterkapitel erklart haben ist der zugehorige Zustand im Wechselwirkungsbild durch
(t) = UD (t, t0 )|ii
(11.3.34)
126
dt HD
(t )|ii
hf |ii hf |
|{z} ~
t0
=0
Z
i t
(11.3.35)
t0
F
ur die Ubergangswahrscheinlichkeit
erhalten wir daher also zu dieser Ordnung
2
Z
1 t i(Ef Ei )t /~
(11.3.36)
Pif = 2 dt e
hf |H (t )|ii .
~
t0
Als einfaches Beispiel betrachten wir eine Storung, die abrupt eingeschaltet wird,
H (t) = V (t t0 ) ,
(11.3.37)
Pi
t2
2 /t
2 / t
if
Der zeitabhangige Vorfaktor von (11.3.38) hat Nullstellen bei if t = 2n und verhalt
sich f
ur if 0 wie
2
(if t/2)2
t
,
2
(~if /2)
~
127
(11.3.39)
Das Gewicht dieser energieerhaltenden Prozesse nimmt linear mit der Zeit zu,
Z
Pif ()d t .
(11.3.40)
2/t
In der Folge betrachten wir ein Kontinuum von Zustanden oder eine Gruppe von Finalzustanden welche dicht liegen. Die Dichte dieser Zustande im Energieraum sei beschrieben
durch (Ef ), das heisst, wir konnen Summen u
ber Zustande |f i ersetzen durch Integrale
u
ber Ef ,
Z
X
dEf (Ef ) .
(11.3.41)
f
Die Wahrscheinlichkeit f
ur einen Ubergang
in einen Zustand um |f i herum ist dann
gegeben durch
Z
X
sin2 (if t/2)
|hf |V |ii|2
t =
Pif (t) =
dEf (Ef )
(~if /2)2
f
2
=
t (Ef ) |hf |V |ii|2
.
(11.3.42)
~
Ei =Ef
|hf |V |ii|2 in der Umgebung von Ef stetig sind. Damit erhalten wir f
ur die Ubergangsrate
dPif /dt =
2
=
|hf |V |ii|2 (Ef ) ,
Fermis Goldene Regel,
(11.3.43)
~
Pi
2 h/t
129
Literatur
[D]
[H]
[DB]
P.C.W. Davies, D.S. Betts, Quantum Mechanics, Chapman & Hall (1994).
[M]
[LMR]
[I]
130