Hopp til innhold

Kvantemekanisk perturbasjonsteori

Fra Wikipedia, den frie encyklopedi
Feynman-diagram som bidrar til elektronets magnetiske moment i laveste orden av perturbasjonsteori.

Perturbasjonsteori benyttes i kvantemekanikk og kvantefeltteori for å komme frem til approksimative svar på beregninger som ikke lar seg gjennomføre på en eksakt måte. Fremgangsmåten er basert på forskjellige metoder som har det til felles at nøyaktigheten til resultatet kan systematisk forbedres, men krever en tilsvarende stadig større beregningsinnsats.

Elektronets magnetisk moment er beregnet i kvanteelektrodynamikk ved bruk av perturbasjonsteori med en nøyaktighet som tilsvarer 1 til 1010. Dette er det mest presise resultat som noensinne er oppnådd i fysikk.

Fremgangsmåten har sine røtter i klassisk himmelmekanikk om planetenes bevegelser der Keplers lover gjorde det mulig å beregne eksakt hvordan én planet beveget seg alene i Solens gravitasjonsfelt. Ved en nøyaktigere beskrivelse må man også ta hensyn til at planetene påvirker hverandre med en mye svakere gravitasjonskraft. Denne sies å forstyrre eller «perturbere» den frie bevegslen, men forandringen kan beregnes så lenge denne perturbasjonen er liten.

De fleste kvantemekaniske problem kan ikke løses eksakt. Men mange ganger kan man likevel finne et lignende og enklere system som er løsbart. Hvis forskjellen mellom dette og det opprinnelige systemet er liten nok, kan man ved perturbasjonsteori finne tilnærmete resultat for egenskapene til det kompliserte problemet.

Ved beregning av det magnetiske momentet til elektronet, antar man først at man kan se bort fra dets kobling til fotoner. Da finner man det største bidraget fra løsningen av Dirac-ligningen for ett elektron. Korreksjonene som er forårsaket av det kvantiserte strålingsfeltet, er små og kan så systematisk beregnes.

Egenskapene til et kvantemekanisk system kan finnes fra løsninger av Schrödinger-ligningen. Den er bestemt av systemets Hamilton-operator som innneholder de variable som beskriver systemet. Er denne uavhengig av tiden, vil en løsning bety å finne eksakte verdier for alle egenverdiene til denne operatoren. De er systemets mulige, kvantiserte energier og bestemmer samtidig hvordan det utvikler seg med tiden. Men for de fleste, realistiske system lar disse seg ikke beregne eksakt. Med perturbasjonsteori kan approksimative verdier ofte finnes for disse størrelsene.

Hvis Hamilton-operatoren derimot varierer med tiden, vil systemet ikke ha noen stasjonære egentilstander, men hele tiden være under forandring. Det er da mindre opplagt hvordan perturbasjonsteori kan benyttes. Men hvis den ikke-konstante delen av operatoren virker bare i et begrenset tidsrom, vil perturbasjonsteori kunne benyttes til å beregne sannsynligheten for at systemet går over fra en stasjonær tilstand til en annen hvor de stasjonære tilstandene er definerte ved den konstante delen av Hamilton-operatoren.[1]

Kvantemekanisk perturbasjonsteori kan derfor deles opp i to hovedklasser. I det statiske tilfellet kan man beregne egenverdier for systemets energi som ikke forandrer seg med tiden. Men når Hamilton-operatoren varierer med tiden, finnes ikke slike permanente egenverdier. Man må da benytte tidsavhengig perturbasjonsteori. Dette er for eksempel tilfelle for problemer som omhandler kollisjoner mellom partikler og er mest benyttet ved anvendelser av kvantefeltteori.[2]

Konstante perturbasjoner

[rediger | rediger kilde]

Når den fulle Hamilton-operatoren for systemet ikke er eksplisitt avhengig av tiden, sies den å være konstant selv om den inneholder dynamiske variable som varierer med tiden. Hvis den i tillegg kan splittes opp i to deler

der den siste delen er mindre viktig, sier man at denne utgjør en konstant perturbasjon av systemet som tilhører Denne Hamilton-operatoren beskriver et enklere system, og man kan anta at dens egentilstander og tilsvarende egenverdier kan beregnes. Det vil si at den tidsuavhengige Schrödinger-ligningen

kan løses. Derimot kan ikke egentilstandene og de tilhørende egenverdiene til Schrödinger-ligningen

finnes. Med antagelsen av at effekten av perturbasjonen er liten, forventes også at energidifferansene er små. Likedan bør de «vekselvirkende» tilstandene ikke bli så forskjellige fra de «frie» tilstandene Hvis dette ikke er tilfelle, kan ikke betraktes som en perturbasjon.[2]

Da egentilstandene til Hamilton-operatoren utgjør et fullstendig sett, kan man allltid skriive en vilkårlig tilstand

hvor koeffisientene Man forventer derfor at i grensen der samtidig som at når

Første ordens korreksjon

[rediger | rediger kilde]

Hvis man projiserer den vekselvirkende Schrödinger-ligningen på den frie tilstanden kan man skrive resultatet som

Her kan man på venstre side benytte at begge tilstandene der er egentilstander for hver sin Hamilton-operator. Det gir

Når er derfor koeffisientene

I tillegg får man når at

Ut fra denne formen for energidifferansen er det hensiktsmessig å normalisere den perturberte tilstanden slik at matriseelementet Det gir det generelle resultatet

for de perturberte energiene.

Man kan nå benytte antagelsen om at perturbasjonen er liten. Da kan den perturberte tilstanden skrives som hvor de prikkete korreksjonsleddene antas å være små. Til laveste orden i perturbasjonen har man da at

Dette sies å være resultatet for denne størrelsen til første orden da det er proporsjonalt med første potens av perturbasjonen.[1]

Alle egenskaper til en kvantisert harmonisk oscillator kan beregnes eksakt. Når dette spesielle potensialet forandres litt, kan man finne korreksjoner til disse resultatene ved bruk av perturbasjonsteori. Dette kan enklest illustreres for grunntilstanden til oscillatoren. Den har den eksakte, uperturberte verdien når ω  er dens vinkelfrekvens. I en koordinatbasis er den tilsvarende , normerte egenfunksjonen

der m  er oscillatorens masse.

Som et eksempel kan man betrakte det perturberende potensialet V = λx4 hvor λ er en liiten parameter. Til første orden blir nå forandringen i grunntilstandsenergien

Integralet her er et utvidet Gauss-integral som kan uttrykkes ved gammafunksjonen som Γ(5/2) = (3/2)Γ(3/2) = (3/4)√π. Dermed blir forandringen av dette laveste energinivået

til første orden i perturbasjonen. På tilsvarende måte kan man finne forskyvningen av høyere energinivå. Beregningen er da spesielt enkel å utføre ved bruk av oscillatorens stigeoperatorer.[3]

Andre ordens korreksjon

[rediger | rediger kilde]

Noen ganger er første ordens resultat null eller man ønsker seg bedre nøyaktighet. Da må man gå till andre ordens perturbasjonsteori som betyr å ta med første ordens korreksjon till den perturberte tilstanden Fra det tidligere kravet at har man da

som er korrekt til første orden i perturbasjonen. Av samme grunn har man derfor

Dette resultatet kan nå benyttes direkte i den generelle formelen til å gi den perturberte energien

Den er nå gyldig til andre orden. Korreksjoner av enda høyere orden kan systematisk beregnes på samme vis.[2]

Høyere ordens korreksjoner innebærer summasjoner over alle tilstandene til systemet. Det dette antallet vanligvis er uendelig stort, vil dette være beregningsmessig vanskelig. Derimot for det enkleste system med bare to tilstander, kan man lett illustrere en beregning til andre orden. For eksempel kan man betrakte et slikt system hvor den frie delen er gitt ved en Hamilton-operator som på matriseform er

mens perturbasjonen er

De to egenverdiene til er og med ortonormerte egenvektorer

For å finne første ordens korreksjoner behøver man matriseelementene og Til denne orden er derfor de perturberte energiene og

Mer nøyaktige resultat forventes til andre orden. Da behøver man også matriseelementet Det gir nå

og

Denne korreksjonen er derfor liten. Med. samme nøyaktighet kan man nå også finne de tilsvarende egentilstandene.

I dette enkle eksemplet kan man finne eksakte løsninger for de samme størrelsene. Den totale Hamilton-operatoren er gitt ved matrisen

med egenverdier og Resultatene i andre ordens perturbasjonsteori har derfor en nøyaktighet som tilsvarer omtrent 5 %. Hva som kreves av perturbasjonen for at denne fremgaøngsmåten skal gi så gode resultat, er vanskelig å besvare.[4]

Degenerasjon

[rediger | rediger kilde]

Andre ordens perturbasjonsteori bryter sammen hvis noen av tilstandene man summerer over, har samme energi som den tilstanden man betrakter. Slike tilstander sies å være degenererte og opptrer ofte i forskjellige sammenhenger. For eksempel, en partikkel som beveger seg i et sentralsymmetrisk potensial, vil ha energier som avhenger av dens kvantiserte dreieimpuls gitt ved kvantetallet ℓ. Men de er uavhengig av det magnetiske kvantetallet m som angir retningen til dreieimpulsen. Dermed inneholder hvert energinivå 2ℓ +1 tilstander med samme energi, men forskjellige verdier av m.

Igjen kan man betrakte en tilstand med energi som er den samme for i alt g lignende tilstander. På grunn av perturbasjonen kan nå forandringen av denne skrives som

hvor tilstanden er ortogonal til de g degenererte tilstandene. Nå skal være en egentilstand til den totale Hamilton-operator med egenverdi Til første orden i perturbasjonen har man derfor

Ved å projisere denne ligningen på en av de degenererte tilstandene vil første ledd ikke bidra da ut fra definisjonen av en slik tilstand. Bidraget fra det andre leddet kan dermed skrives som

etter å ha definert matriseelementene Dette utgjør et homogent, lineært system med g ligninger for de ukjente koeffisientene as. Ikke-trivielle løsninger finnes da bare når determinanten

Denne ligningen har g røtter som er resultatet for energiforskyvningene Noen av dem kan fortsatt være sammenfallende. Dette gir en oppsplitting av det degenererte energinivået i maksimalt g nye nivå.[4]

I praksis kan mye av beregningsarbeidet i dette tilfellet reduseres ved å velge basistilstander som er spesielle lineærkombinasjoner av de degenererte tillstandene. Matrisen kan på den måten gjøres mest mulig diagonal før dens egenverdier bestemmes fra determinanten..[3]

Tidsavhengig perturbasjonsteori

[rediger | rediger kilde]

En generell tilstand for systemet kan for hvert tidspunkt t  alltid utvikles i egentilstander til den frie Hamilton-operatoren med egenverdier Ved tiden t 0 har man derfor

hvor kvadratet av amplitudene er sannsynligheten for å finne systemet i tilstanden I det spesielle tilfellet at perturbasjonen er systemet ved et senere tidspunkt nå beskrevet av tilstandsvektoren

Sannsynligheten for å finne systemet ved denne tiden i tilstand er derfor fremdeles det vil si den samme som ved tiden t 0. Hvis systemet er i en egentilstand til ved ett tidspunkt, er det i den samme egentilstanden ved et senere tidspunkt. Det sies å være stasjonært uten overganger mellom slike tilstander.[5]

Overgangsamplituder

[rediger | rediger kilde]

Når dette systemet er utsatt for en perturbasjon er tidsutviklingen av tilstandsvektoren styrt av den totale Hamilton-operatoren Den vil da bli ganske annerledes enn for men kan likevel formelt skrives på formen

Det betyr at sannsynligheten for at systemet befinner seg i tilstand varierer med tiden. Hvis systemet er i en bestemt tilstand ved tiden t 0, har det en viss sannsynlighet for å bli funnet i en helt annen tilstand ved et senere tidspunkt. Disse koeffisientene kalles derfor i denne sammenheng for «overgangsamplituder». Da perturbasjonen antas å være tilstrekkelig svak, er det å forvente at tidsvariasjonen til disse amplitudene er mye langsommere enn svingningene gitt ved de frie egenverdiene som opptrer i eksponensialfunksjonene.[6]

Variasjonen med tiden til tilstandsvektoren og dermed også amplitudene, er gitt ved den fulle Schrödinger-ligningen

Fra den antatte formen til tilstandsfunksjonen har man nå at

da de frie egentilstandene ikke varierer med tiden. Samtidig gir høyresiden i ligningen

Sammenligning av de to sidene fører dermed til betingelsen

Projiserer man denne ligningen på en egentilstand og benytter at disse er ortonormerte slik at kommer man frem til resultatet

Dette er et sett med koblete, første ordens differensialligninger for de ukjente overgangsamplitudene. De kan la seg beregne systematisk som en Dyson-serie ved bruk av det ekvivalente vekselvirkningsbildet.[5]

Første ordens overgang

[rediger | rediger kilde]

Differensialligningene for overgangsamplitudene er ekvivalent med den opprinnelige Schrödinger-ligningen for det koblete systemet da ingen approksimasjonener gjort under deres utledning. Men de er bedre egnet til å beskrive hva en slik perturbasjon vil gjøre med et system som befinner seg i en gitt tilstand ved et tidlig tidspunkt t 0. Da er amplitudene gitt som det vil si med null sannsynlighet for å være i andre tilstander. Så snart perturbasjonen starter å virke, vil alle disse forandre seg. Men ved å anta at denne forstyrrelsen er svak, kan man til laveste orden sette inn på høyresiden av ligningene. Amplituden for å finne systemet i en annen tilstand ved et senere tidspunkt t > t 0 er da gitt som

hvor matriseelementet er Så snart dette er beregnet, gir integralet amplituden til første orden i perturbasjonen. En mer nøyaktig beregning kan så eventuelt gjennomføres ved å benytte dette resultatet på høyre side i differensialligningene. Det vil da gi overgangsamplituden til andre orden i perturbasjonen.[7]

En enkel illustrasjon av denne fremgangsmåten kan gis for en kvantisert harmonisk oscillator. Den er beskrevet ved Hamilton-operataoren

og har egentilstander med egenverdier for n = 0, 1, 2 og så videre. Man kan anta at den opprinnelige befinner seg i grunntilstanden med kvantetall n = 0. Rundt tiden t = 0 blir den så utsatt for et lite støt som er beskrevet ved perturbasjonen

hvor k er en konstant og τ angir varigheten av støtet. Den er veldig svak ved tidlige tider, har et maksimum ved tiden t = 0 og blir så raskt svakere igjen.

Ved å benytte oscillatorens stigeoperatorer, er da matriseelementet for en overgang til en eksitert tilstand n > 0 gitt som

Her er nå og slik at matriseelementet er forskjellig fra null kun for n = 1. Dette er eksempel på en «utvalgsregel». Eksitasjon til høyere nivå vil kunne fremkomme i høyere ordens perturbasjonsteori.[6]

Overgangsamplituden til første energinivå er nå gitt som

Ved et mye senere tidspunkt t → ∞ kan dette finnes fra det bestemte Gauss-integralet

Da blir overgangssannsynligheten

og er størst når støtet varer mye kortere enn perioden T = 2π /ω til oscillatoren.

Harmonisk perturbasjon

[rediger | rediger kilde]

For å studere hvordan et kvantemekanisk system vekselvirker med elektromagnetisk stråling, består denne av bølger. Hvis man betrakter bidraget fra en frekvens ω, vil det tilsvare en tidsvariabel perturbasjon

hvor er en konstant operator. Den kan generelt være avhengig av andre operatorer som beskriver systemet under betraktning.

Amplituden for at systemet skal gå fra en opprinnelig tilstand til en annen tilstand ved et mye senere tidspunkt, blir da til første orden i denne perturbasjonen

der nå Argumentet til Diracs deltafunksjon som opptrer her, avhenger av om man betrakter det første eller det andre leddet i uttrykket for perturbasjonen. I det ene tilfellet gir den at Det betyr at systemet i sluttilstanden er tilført energien som er naturlig å kalle et kvant. I det andre tilfellet avgir systemet den samme energien som tilsvarer emisjon av et slikt kvant.[7]

Sannsynligheten for at denne overgangen skal finne sted, er nå gitt ved det absolutte kvadratet

Det inneholder produktet av to singulære deltafunksjoner som kan reguleres ved å skrive den ene som grensen av et regulært integral,

hvor man må la T → ∞. Den gjenværende deltafunksjonen fører nå til at integralet gaanske enkelt gir T/2 - (-T/2) = T. Overgangssannsynligheten øker derfor proporsjonalt med tiden prosessen varer. Antall overganger per tidsenhet er dermed konstant og bestemt av

Det engelse ordet for denne «overgangshastigheten» er transition rate.

Tidsavhengig perturbasjonsteori ble først utviklet og benyttet av Paul Dirac i 1927. Han kunne på den måten beregne Einsteins strålingskoeffisienter og dermed gi en kvantemekanisk forklaring på absorpsjon og emisjon av lys fra atomer.[8]

Eksempel: To-nivå system

[rediger | rediger kilde]

I det enkleste tilfellet kan man tenke seg et system med bare to tilstander. Det befinner seg opprinnelig i den laveste tilstanden med energi E1 < E2 og man utsetter det for en harmonisk perturbasjon med A21 = ħa ved tiden t = 0. Da vil amplituden for å finne det i den øvre tilstanden ved et senere tidspunkt, være i laveste orden gitt ved

der man kan sette c1 = 1. En enkel integrasjon gir dermed

når man setter Sannsynligheten for at systemet vil befinne seg i den eksisterte tilstanden kan nå skrives som

og vil oscillere med en frekvens Når frekvensen ω til den påtrykte perturbasjonen nærmer seg systemfrekvensen Δ, blir denne sannsynligheten større enn én. Første ordens perturbasjonsteori kan derfor ikke benyttes og må utvides med en mer nøyaktig beregning. For dette spesielle systemet kan man også finne en eksakt løsning med samme form som denne, men med en frekvens Ω som inneholder amplituden a  til perturbasjonen.[2]

Fermis gyldne regel

[rediger | rediger kilde]

Deltafunksjonen i den generelle overgangssannsynligheten er et uttrykk for at energien er bevart i prosessen. Den er strengt tatt kun gyldig for en diskret slutttilstand mens det ofte vil finnes et kontinuum av tilstander med nesten sammenfallende energier . Den totale overgangssannsynlighet er da gitt ved en sum over disse tilstandene. Med de samme antagelsene kan den da erstattes med integralet

hvor fremdeles gjelder. Her er funksjonen tettheten av tilstander med energi E. Det betyr at deres antall i et lite intervall kan skrives som

Resultatet for overgangshastigheten kalles for «Fermis gyldne formel» da den kan anvendes på svært mange prosesser. Den har samme form uansett hvor mange partikler eller atomer som finnes ved begynnelsen eller som opptrer i slutttilstanden. Også hvis disse beveger seg med relativistiske hastigheter, vil den se likedan ut. Den har samme form ved bruk av høyere ordens perturbasjonsteori som kun vil forandre matriseelementet . Men dette og tetthetsfunksjonen vil forandre seg fra prosess til prosess.[5]

Potensialspredning

[rediger | rediger kilde]

Et viktig eksempel på bruk av Fermis gyldne regel er spredning av en ikke-relativistisk partikkel på et statisk potensial Partikkelen som kommer inn mot potensialet, har impuls og er derfor beskrevet ved bølgefunksjonen

Spredningstverrsnittet for prosessen angir sannsynligheten for at den kommer ut i en slutttilstand med impuls som tilsvarer bølgefunksjonen

Til første orden er da matriselementet av det perturberende potensialet gitt ved integral

hvor gir forandringen av impulsen til partiklen. Matriseelementet er derfor ikke noe annet enn det Fourier-transformerte potensialet, det vil si

Spredningstverrsnitt

[rediger | rediger kilde]

Istedenfor å regne ut tettheten av slutttilstander i Fermis regel er det enklere å benytte at denne tettheten i impulsrommet er

hvor den infinitesemale romvinkelen angir retningen til partikkelen i slutttilstanden. Da dens energi da er slik at kan man integrere over denne energien. På den måten forsvinner deltafunksjonen etter å ha gitt energibevarelse som karakteriserer elastisk spredning uten energitap.[3]

Den differensielle overgangssannsynligheten for at den innkommende partikkelen skal bli spredt i retningen er dermed

Fluksen eller strømtettheten av innkommende partikler finnes fra sannsynlighetsstrømmen beregnet fra bølgefunksjonen Når impulsen velges å være langs z-aksen, blir denne ganske enkelt hvor for elastisk spredning. Dermed kan det differensielle spredningstverrsnittet finnes fra

Det avhenger av retningen til den spredte partikkelen gjennom impulsforandringen Sammenhengen følger fra

der spredningsvinkelen er definert som vinkelen mellom innkommende impuls og utgående impuls Begge impulsene er like store og gitt som hvor er hastigheten til partiklene langt unna potensialet hvor de antas å bevege seg fritt.[6]

For Coulomb-spredning har potensialet formen som etter Fourier-transformasjon gir

Det gir i et differensielt spredningstverrsnitt som kan skrives på formen

hvor Plancks konstant har fallt ut. En første ordens beregning i kvantemekanisk perturbasjonsteori gir derfor det samme spredningstverrsnittet som for klassisk Rutherford-spredning.

Referanser

[rediger | rediger kilde]
  1. ^ a b D.J. Griffiths, Quantum Mechanics, Pearson Education International, Essex (2005). ISBN 1-292-02408-9.
  2. ^ a b c d B.H. Bransden and C.J. Joachain, Quantum Mechanics, Prentice Hall, New York (2000). ISBN 978-0-582-35691-7.
  3. ^ a b c R. Shankar, Principles of Quantum Mechanics, Plenum Press, New York (1982). ISBN 0-306-40397-8.
  4. ^ a b A. Messiah, Quuantum Mechanics, Volume II, North-Holland Publishing Company, Amsterdam (1962).
  5. ^ a b c J.J. Sakurai, Modern Quantum Mechanics, Benjamin Cummings Publishing Company, Menlo Park CA (1985). ISBN 0-8053-7501-5.
  6. ^ a b c E. Merzbacher, Quantum Mechanics, John Wiley & Sons, New York (1961).
  7. ^ a b J.J. Sakurai, Advanced Quantum Mechaniics, Addison -Wesley Publishing Company, Reading MA (1967).
  8. ^ W. Heitler, The Quantum Theory of Radiation, Oxford Universiity Press, England (1960).

Eksterne lenker

[rediger | rediger kilde]