Fermis Goldene Regel
Fermis Goldene Regel, benannt nach dem Physiker Enrico Fermi, bezeichnet eine viel benutzte Gleichung aus der quantenmechanischen Störungstheorie. Die Gleichung ergibt die theoretische Voraussage für die Übergangsrate (Übergangswahrscheinlichkeit pro Zeit), mit der ein Anfangszustand unter dem Einfluss einer Störung in einen anderen Zustand übergeht. Wenn nicht zusätzlich noch Übergänge in weitere Zustände möglich sind, gibt der Kehrwert der Übergangsrate die mittlere Lebensdauer des Anfangszustands an. Anschaulich gesagt ist das die Zeit, die der Quantensprung in den neuen Zustand im Mittel noch auf sich warten lassen wird. Mit der Goldenen Regel kann man sowohl spontane Umwandlungen behandeln (z. B. den radioaktiven Zerfall, die Emission von Lichtquanten, den Zerfall von instabilen Elementarteilchen) als auch die Absorption (z. B. von Lichtquanten), aber auch den Wirkungsquerschnitt beliebiger Reaktionen zwischen zwei Teilchen.
Formel
Wird ein Anfangszustand $ |i\rangle $ einer Störung $ V $ ausgesetzt, durch die er in einen Endzustand $ |f\rangle $ in einem Energiekontinuum übergehen kann, so ist in 1. störungstheoretischer Näherung die Übergangsrate dafür (d. h. die Übergangswahrscheinlichkeit pro Zeiteinheit) durch
- $ \lambda _{i\to f}={\frac {2\pi }{\hbar }}\rho (E_{f})\left|V_{fi}\right|^{2} $
gegeben. Dabei ist $ \hbar $ das reduzierte plancksche Wirkungsquantum, $ \rho (E_{f}) $ die Zustandsdichte der beobachteten Endzustände bei der Energie $ E_{f} $, und $ V_{fi} $ das zu diesem Übergang gehörende Matrixelement des Störoperators $ V $. Die Zustände $ |i\rangle $ und $ |f\rangle $ sind Eigenzustände zu einem ungestörten Hamilton-Operator, zu dem die Störung $ V $ zusätzlich betrachtet werden soll. Die Übergangsrate $ \lambda _{i\to f} $ hat die Dimension 1/Zeit. Für spontane Zerfälle (Beispiel: Radioaktivität) ist $ \lambda _{i\to f} $ die Zerfallskonstante im Exponentialgesetz. $ \Gamma _{i\to f}=\hbar \lambda _{i\to f} $ hat die Dimension Energie und gibt die Energieunschärfe oder Halbwertsbreite des Anfangszustands an. Sind Umwandlungen in verschiedener Form möglich, ergeben sich die totale Zerfallskonstante $ \lambda $ bzw. die totale Halbwertsbreite $ \Gamma $ aus der Summe der einzelnen partiellen Werte für jede Art des Übergangs. $ \tau =1/\lambda $ ist die mittlere Lebensdauer des Systems im Anfangszustand.
Geschichte
Der störungstheoretische Formalismus der „Goldenen Regel“ wurde 1927 von Paul Dirac entwickelt, um die Absorption und Emission von Photonen erstmals quantenmechanisch zu behandeln[1]. Wenig später wurde er von Gregor Wentzel in einer Arbeit zur Berechnung der Übergangswahrscheinlichkeit für den (strahlungslosen) Auger-Meitner-Effekt in Atomen noch einmal entwickelt, der ebenfalls ein Übergang von einem diskreten Zustand im Atom in den kontinuierlichen Bereich des Spektrums ist. [2] Nach Fermi ist diese "Regel"[Anmerkung 1] benannt, da er sie 1950 in einem Kernphysik-Lehrbuch als „Golden Rule No. 2“[3] aufführte. In der Literatur finden sich aber manchmal auch die Bezeichnungen Wentzel-Fermi Golden Rule und Fermi-Wentzel Golden Rule.
Als Golden Rule Nr.1 wird bei Fermi[4] der Einsatz des Terms 2. Ordnung Störungstheorie für solche Übergänge bezeichnet, die nach der 1. Ordnung verboten wären.
Anwendungen
Auf Grund ihrer allgemeinen Gültigkeit können für Fermis Goldene Regel vielfältige Anwendungen gefunden werden, z. B. in der Atomphysik, Kernphysik und Festkörperphysik bei der Absorption und Emission von Photonen, Phononen oder Magnonen.
Herleitungsskizze
Als Grundannahme wird ein zeitlich konstantes System mit exakt lösbarem Hamiltonoperator $ \,H_{0} $ durch einen Störoperator $ \,V $ erweitert.
- $ \,H=H_{0}+V $
Fermis Goldene Regel gilt für beliebige, konstante oder zeitabhängige Störoperatoren. Sie können die Wechselwirkung mit einem äußeren (konstanten oder zeitabhängigen) Feld darstellen oder eine zusätzliche Art der Wechselwirkung zwischen den Teilchen des Systems, die in $ \,H_{0} $ nicht berücksichtigt war (z.B. die Möglichkeit der Erzeugung eines Photons). Hier wird die Herleitung für einen zeitlich konstanten Störoperator gezeigt.
Für den vervollständigten Hamiltonoperator muss die zeitabhängige Schrödingergleichung
- $ \mathrm {i} \hbar \partial _{t}\left|\Psi (t)\right\rangle =H\left|\Psi (t)\right\rangle $
gelöst werden. Am Anfang (t = 0) soll das System sich in einem Eigenzustand zu $ \,H_{0} $ befinden. Man entwickelt die gesuchte Funktion $ \left|\Psi (t)\right\rangle $ nach den Eigenfunktionen des ungestörten Hamiltonoperators $ \left|\Psi _{n}\right\rangle $ (mit Energie-Eigenwerten $ E_{n} $) und hat damit zeitabhängige Koeffizienten:
- $ \left|\Psi (t)\right\rangle =\sum _{n=0}^{\infty }a_{n}(t)\,\left|\Psi _{n}\right\rangle \,e^{-{\frac {\mathrm {i} E_{n}t}{\hbar }}} $.
Die Anfangsbedingung ist $ a_{i}(0)=1 $, alle anderen $ a_{n}(0)=0 $.
Nach Einsetzen von Hamiltonoperator und Wellenfunktion in die Schrödingergleichung ergibt sich durch Koeffizientenvergleich:
- $ \mathrm {i} \hbar \,\partial _{t}a_{f}(t)=\sum _{n=0}^{\infty }V_{fn}\,a_{n}(t)\,e^{i\omega _{fn}t} $,
wobei $ \omega _{fn} $ eine Kurzschreibweise für $ {\tfrac {E_{f}-E_{n}}{\hbar }} $ darstellt.
Diese Gleichung beschreibt, wie sich die Koeffizienten $ a_{f}(t) $ zeitlich ändern. Zur näherungsweisen Lösung wird angenommen, dass die Koeffizienten sich gegenüber ihren Anfangswerten stetig ändern, sodass für kleine Zeiten von der Summe nur das Glied n = i zu betrachten ist. (Das ist der Sinn der 1. störungstheoretischen Näherung.) Da $ V_{fi} $ zeitlich konstant ist, kann die Gleichung integriert werden:
- $ \mathrm {i} \hbar a_{f}(t)=V_{fi}\int _{0}^{t}dt'e^{i\omega _{fi}t'} $
Als Ergebnis erhält man
- $ \mathrm {i} \hbar a_{f}(t)=V_{fi}\,{\frac {1}{i\omega _{fi}}}\,(e^{i\omega _{fi}t}-1)=2V_{fi}\;e^{\frac {i\omega _{fi}t}{2}}\,\left({\frac {\sin {\tfrac {\omega _{fi}t}{2}}}{\omega _{fi}}}\right) $.
Es sei daran erinnert, dass Aufgrund der Herleitung diese Formel nur gültig bleiben kann, solange $ \left|a_{f}(t)\right|\ll 1 $ gilt. Die Wahrscheinlichkeit, zur Zeit t das System im Zustand f zu finden, ist das Betragsquadrat
- $ \left|a_{f}(t)\right|^{2}={\frac {4\left|V_{fi}\right|^{2}}{\hbar ^{2}}}\left({\frac {\sin ^{2}{\tfrac {\omega _{fi}t}{2}}}{\omega _{fi}^{2}}}\right) $
Wir betrachten einen Übergang ins Kontinuum, so dass der Endzustand f zahlreiche Nachbarzustände ähnlicher Struktur, aber mit kontinuierlich variierender Energie besitzt, die auch als Endzustände möglich sind. Das Matrixelement kann deshalb für alle gleich angenommen werden, die jeweilige Übergangswahrscheinlichkeit ist wegen des eingeklammerten Faktors (in nachfolgender Graphik $ g $) in der letzten Formel aber unterschiedlich. In Abhängigkeit von $ \omega _{fi} $ betrachtet, ist dieser Faktor $ g $ eine Funktion mit einem spitzen Maximum der Höhe $ {\tfrac {t^{2}}{4}} $ bei $ \omega _{fi}=0 $. Die benachbarten Nullstellen liegen bei $ \omega _{fi}=\pm {\tfrac {2\pi }{t}} $, das Maximum dazwischen kann man gut durch ein Dreieck mit der Grundlinie $ {\tfrac {4\pi }{t}} $ annähern. Der Mittelwert des Klammerfaktors in diesem Intervall ist daher näherungsweise die Hälfte des Maximums, $ {\tfrac {t^{2}}{8}} $.
Das spitze Maximum dieser Funktion bei $ \omega _{fi}=0 $ zeigt, dass Übergänge bevorzugt in Zustände gleicher Energie $ E_{f}=E_{i} $ führen, aber Abweichungen in einem Bereich der Breite $ \Delta E={\tfrac {4\pi \hbar }{t}}={\tfrac {2h}{t}} $ möglich sind. Mit wachsender Zeit wird dieser Schwankungsbereich kleiner. Dies ist eine der Formen der Unschärferelation für Energie und Zeit und begründet die bei allen Spektrallinien zu findende natürliche Linienbreite. Auch wird mit wachsendem t die Spitze höher, und die gesamte Fläche unter dem Maximum (genähert durch halbe Höhe mal Breite) wächst proportional zu t an. Ist W die Summe der Übergangswahrscheinlichkeiten in alle Zustände im Bereich des Maximums, wächst also W proportional zu t, und $ {\tfrac {W}{t}} $ ist die gesuchte konstante Übergangsrate. Aber es sei wiederholt, dass die ganze Herleitung nur solange gültig bleiben kann, wie $ W\ll 1 $ bleibt.
Zur Berechnung von W wird nun die mittlere Übergangswahrscheinlichkeit im Maximum einfach mit der Zahl der Endzustände im Intervall $ \Delta E $ multipliziert. Diese Zahl der ergibt sich aus der Dichte der Zustände auf der Energieachse
- $ \rho (E_{f})={\frac {dn}{dE_{f}}} $
multipliziert mit der Breite $ \Delta E={\tfrac {4\hbar \pi }{t}} $.
Für die Summe W aller einzelnen Übergangswahrscheinlichkeiten $ \left|a_{f}(t)\right|^{2} $ erhält man damit
- $ W={\frac {4\left|V_{fi}\right|^{2}}{\hbar ^{2}}}\left({\frac {t^{2}}{8}}\right)\left({\frac {4\hbar \pi }{t}}\right)\rho (E_{f})\ =\ {\frac {2\pi }{\hbar }}\;\rho (E_{f})\;\left|V_{fi}\right|^{2}\;t $.
Division durch die Zeit ergibt somit Fermis Goldene Regel für die Übergangsrate:
- $ \lambda _{i\to f}={\frac {2\pi }{\hbar }}\rho (E_{f})\left|V_{fi}\right|^{2} $
Literatur
Auf Grund ihrer Wichtigkeit für die quantenmechanische Störungstheorie wird Fermis Goldene Regel in den meisten einführenden Büchern zur Quantenmechanik behandelt.
- J.J. Sakurai: Modern Quantum Mechanics. Addison-Wesley, 1994, ISBN 0-201-53929-2.
- W. Greiner: Quantenmechanik Teil I - Eine Einführung. 4. Auflage. Harry Deutsch, Thun 1989, ISBN 978-3-81711-064-3.
Weblinks
- Transition Probabilities and Fermi's Golden Rule
- Ausführliche Herleitung (PDF, engl.)
- Freier Matlab-Code zur Visualisierung der Regel
Anmerkungen
- ↑ Warum sie als "Regel" bezeichnet wird anstatt "Gleichung" oder "Formel" wie sonst eher üblich, ist nicht klar. Jedenfalls kann dies Wort an die Anfangszeit der neuen Quantenmechanik erinnern, als noch herumprobiert werden musste, welche Rezepte oder eben Regeln zum Erfolg führen.
Einzelnachweise
- ↑ Paul Dirac The Quantum Theory of the Emission and Absorption of Radiation, Proceedings Royal Society A, Band 114, 1927, S. 243
- ↑ Gregor Wentzel: Über strahlungslose Quantensprünge. In: Zeitschrift für Physik. 43, Nr. 8, 1927, S. 524–530, doi:10.1007/BF01397631.
- ↑ Enrico Fermi: Nuclear Physics. University of Chicago Press, 1950, S. 142.
- ↑ Fermi, loc. cit. S. 136