Das Wasserstoffatom (relativistsich)

Aus PhysikWiki
Zur Navigation springen Zur Suche springen

__SHOWFACTBOX__

In einem rotationssymmetrischen Potenzial haben wir als Dirac- Hamiltonian:

H=(cα¯p¯+m0c2β+V(r))pr:=1r(r¯p¯−iℏ)αr:=1rα¯r¯ℏQ:=β(σ¯~L¯+ℏ)

Dabei sind pr,αr,ℏQ

hermitesche Operatoren

Man kann den Hamilton- Operator schreiben als:

H=(cαrpr+icrαrβℏQ+m0c2β+V(r))

Beweis:

αrpr+irαrβℏQ=αr[1r(r¯p¯−iℏ)+irβ2(σ¯~L¯+ℏ)]β2=1=αrr(r¯p¯+iσ¯~L¯)=1r2[(α¯r¯)(r¯p¯)+i(α¯r¯)(σ¯~L¯)]i(α¯r¯)(σ¯~L¯)=i(α¯r¯)(r¯p¯)−ir2(α¯p¯)⇒αrpr+irαrβℏQ=1r2[(α¯r¯)(r¯p¯)+i(α¯r¯)(σ¯~L¯)]=α¯p¯

Es gilt weiter:

[ℏQ,H]=0

.

Somit existieren gemeinsame Eigenzustände zu H

und ℏQ

Eigenwerte von ℏQ

(ℏQ)2=β(σ¯~L¯+ℏ)β(σ¯~L¯+ℏ)=β2(σ¯~L¯+ℏ)2[β,σ¯~]=0=([1,σ¯~]−[1,σ¯~])β2=1⇒(ℏQ)2=(σ¯~L¯)(σ¯~L¯)+2ℏ(σ¯~L¯)+ℏ2(σ¯~L¯)(σ¯~L¯)=L2+iσ¯~(L¯×L¯)(L¯×L¯)=iℏL¯⇒(σ¯~L¯)(σ¯~L¯)=L2−ℏσ¯~(L¯)

Somit:

(ℏQ)2=L2+ℏσ¯~L¯+ℏ2=(L¯+ℏ2σ¯~)2+ℏ24mit (L¯+ℏ2σ¯~)2=L2+ℏσ¯~L¯+ℏ24σ¯~2σ¯~2=3(L¯+ℏ2σ¯~)=J¯

Schließlich also

(ℏQ)2=J¯2+ℏ24

Die Eigenwerte von J¯2

sind jedoch bekannt, nämlich ℏj(j+1)

mit j=l±s=12,32,...

(ℏQ)2|j⟩=(ℏ2j(j+1)+ℏ24)|j⟩=ℏ2(j+12)2|j⟩(j+12)2:=q2

Somit:

(ℏQ)|j⟩=(ℏq)|j⟩q=±1,±2,...

Es bleibt das radiale Eigenwertproblem für

H=(cαrpr+icrαrβℏQ+m0c2β+V(r))

Geeignete Darstellung für αr

(αr)2=1r2(α¯r¯)(α¯r¯)=1r2αμανxμxν=12r2(αμαν+αναμ)xμxν(αμαν+αναμ)=2δμν12r22xμxμ=r2r2=1αrβ+βαr=1r(α¯β+βα¯)r¯(α¯β+βα¯)=0⇒1r(α¯β+βα¯)r¯=0

Für

β=(100−1)

kann dies durch die Darstellung αr=(0−ii0)

mit αr=αr+

erfüllt werden:

αrβ=(0ii0)βαr=(0−i−i0)

Es gilt:

pr=1r(r¯p¯−iℏ)r¯p¯=ℏir∂∂rpr=1r(ℏir∂∂r−iℏ)=−iℏ(∂∂r+1r)

Also

H=ℏc(0−110)(∂∂r+1r)−cℏqr(0110)+m0c2(100−1)+V(1001)

Ansatz für den Radialanteil

(ϕaϕb) ~1r(F(r)G(r))

Eingesetzt in die Eigenwertgleichung für H:

(ϕaϕb) ~1r(F(r)G(r))

folgt:

−ℏcrdGdr−cℏqr2G+m0c2rF+VrF=EFrℏcrdFdr−cℏqr2F−m0c2rG+VrG=EGrV=−e24πε01r

Also:

(E−m0c2−V)F+ℏcdGdr+cℏqrG=0(E+m0c2−V)G−ℏcdFdr+cℏqrF=0V=−e24πε01r

Skalentransformation:

a1=m0c2+Eℏca2=m0c2−Eℏca=a1a2=m02c4−E2ℏc

Führt man des weiteren ein:

ρ:=arγ:=e24πε0ℏc≈1137

Also einen skalierten Radius und die Feinstrukturkonstante,

wodurch sich auch das Potenzial vereinfacht zu:

Vℏca=−γρ
(ddρ+qρ)G−(a2a−γρ)F=0(ddρ−qρ)F−(a1a+γρ)G=0

Randbedingung:

F(ρ),G(ρ)

regulär bei ρ→0

F(ρ),G(ρ)→0

für ρ→∞

Betrachte |E|<m0c2⇒a1,a2>0a∈R

also gebundene Zustände

Asymptotisches Verhalten:

ρ→∞⇒G ´=a2aFF ´=a1aG⇒G ´ ´=G,F ´ ´=F⇒G=e−ρ=F=G=e−ρ

Weil e+ρ

divergiert!

ρ→0⇒G ´+qρG+γρF=0F ´−qρF−γρG=0

Ansatz:

F(ρ)=f0ρλG(ρ)=g0ρλ⇒(λ+q)g0+γf0=0(λ−q)f0−γg0=0

Es existieren nichttriviale Lösungen f0,g0 ,

falls (λ+q)(λ−q)+γ2=λ2−q2+γ2=0

Also λ=±q2−γ2>0

und regulär bei ρ→0

Ansatz:

F(ρ)=ρλe−ρf(ρ)G(ρ)=ρλe−ρg(ρ)⇒g ´−g+λ+qρg−(a2a−γρ)f=0f ´−f+λ−qρf−(a1a+γρ)g=0

Die Lösung erfolgt über einen Potenzreihenansatz:

f(ρ)=∑k=0∞fkρk⇒f ´(ρ)=∑k=1∞kfkρk−1=∑k=0∞(k+1)fk+1ρkg(ρ)=∑k=0∞gkρk⇒g ´(ρ)=∑k=1∞kgkρk−1f(ρ)ρ=∑k=0∞fkρk−1=f0ρ+∑k=0∞fk+1ρk

usw... wird dies ebenfalls für g ´(ρ),g(ρ)ρ

aufgestellt

Koeffizientenvergleich liefert:

O(1ρ):(λ+q)g0+γf0=0(λ−q)f0−γg0=0⇒f0,g0

bis auf Normfaktor

O(ρk):(λ+q+k+1)gk+1−gk+γfk+1−a2afk=0(λ−q+k+1)fk+1−fk+γgk+1−a1agk=0

k=0,1,2,.... Rekursionsformel!!

a[(λ+q+k+1)gk+1−gk+γfk+1−a2afk]−a2[(λ−q+k+1)fk+1−fk+γgk+1−a1agk]=0⇒[a(λ+q+k+1)+a2γ]gk+1=[a2(λ−q+k+1)−aγ]fk+1

Verhalten für große k:

akgk+1≈a2kfk+1⇒fk≈aa2gk

Dies kann man einsetzen in

(λ+q+k+1)gk+1−gk+γfk+1−a2afk=0

und es folgt:

(k+1)gk+1≈2gk⇒gk+1gk≈2k+1⇒gk+1≈2k+1(k+1)!g0⇒g(ρ) ~e2ρ⇒f(ρ) ~e2ρ

Falls die Potenzreihen

f(ρ)=∑k=0∞fkρk,g(ρ)=∑k=0∞gkρk

nicht abbrechen, so divergiert F(ρ)=ρλe−ρf(ρ)G(ρ)=ρλe−ρg(ρ)

exponentiell für ρ→∞⇒F(ρ),G(ρ) ~eρ

Dies ist jedoch ein Widerspruch zu den gesetzten Randbedingungen!

Also muss es einen Abbruch bei k=n ´∈N

geben:

fn ´+1=gn ´+1=0

Setzt man dies in die Rekursionsformel ein, so folgt:

−gn ´−a2afn ´=0⇒a2fn ´=−agn ´−fn ´−a1agn ´=0⇒afn ´=−a1gn ´

Diese beiden Gleichungen stimmen jedoch für alle f,g überein, da

a2a=aa1

Setzt man a2fn ´=−agn ´

in [a(λ+q+k+1)+a2γ]gk+1=[a2(λ−q+k+1)−aγ]fk+1

ein, so folgt mit k+1=n ´

a(λ+q+n ´)+a2γa=−[a2(λ−q+n ´)−aγ]a2λ+q+n ´+a2aγ+λ−q+n ´+aa2γ=02a(λ+n ´)=(a2a2−a2)γ=2Eℏcγa2a2=a1a2(λ+n ´)2=E2ℏ2c2γ2

Weiter gilt:

a2=m02c4−E2ℏ2c2⇒(m02c4−E2)(λ+n ´)2=E2γ2

Löst man dies nach den exakten Energieeigenwerten, die sich damit ergeben, also nach E auf, so erhält man die Feinstrukturformel:

E=m0c21+(γλ+n ´)2

Mit der Feinstrukturkonstanten

γ≈1137
λ=qa2=m02c4−E2ℏ2c2⇒(m02c4−E2)(λ+n ´)2=E2γ2
λ=q2−γ2=(j+12)2−γ2j=12,32,...,n ´∈N0j=l±s

entwickelt man die Energieeigenwerte nach der Feinstrukturkonstanten bis O(γ4) ,

so folgt:
E=m0c2[1−12(γλ+n ´)2+38(γλ+n ´)4+O(γ6)]

mit

λ(γ)=|q|1−(γq)2=|q|[1−12(γq)2]+O(γ4)
(1λ+n ´)2=1[n ´+|q|−12(γ2|q|)]2+O(γ4)n=n ´+|q|n ´=0,1,2,...|q|=j+12=1,2,....(1λ+n ´)2=1n2[1−12(γ2|q|n)]−2+O(γ4)=1n2[1+(γ2|q|n)]+O(γ4)=1n2+(γ2|q|n3)+O(γ4)|q|=j+12=l±s+12

Setzt man dies in die exakten Energieeigenwerte E ein, so folgt:

E=m0c2[1−(γ22n2)−(γ42n3)(1j+12−34n)+O(γ6)]n=1,2,3j=12,32,...,n−12,wegen n=n ´+j+12j=l±s

Diskussion

O(γ0):E=m0c2

Ruheenergie

O(γ2):ΔE(2)=−m0c2(γ22n2)=−RHn2

nicht relativistisches, entartetes Energiespektrum

O(γ4):ΔE(4)=−m0c2(γ42n3)(1j+12−34n)

Feinstruktur- Aufspaltung. Eine Aufhebung der j-Entartung durch Spin- Bahn- Kopplung. Dabei bleibt die Freiheit der Ausrichtung der Achse des magnetischen Moments, also die 2(2j+1) - fache mj - Entartung+ Parität!

Spektroskopische Beziehung der Feinstrukturterme: nlj

n=1:j=12: 1s12
n=2:j=12: 2s122p12n  ´=1j=32:2p32n  ´=0

n´=0. .