Bethe-Formel

Die Bethe-Formel (auch Bethe-Gleichung, Bethe-Bloch-Formel, Bethe-Bloch-Gleichung oder Bremsformel) gibt den Energieverlust pro Weglängeneinheit an, den schnelle geladene schwere Teilchen (z. B. Protonen, Alphateilchen, Ionen) beim Durchgang durch Materie durch inelastische Stöße mit den Elektronen erleiden; die übertragene Energie bewirkt im Material Anregung oder Ionisation. Dieser Energieverlust, auch als elektronische Abbremsung oder ungenau als Ionisationsverlust bezeichnet, hängt ab von Geschwindigkeit und Ladung der Projektilteilchen und vom Targetmaterial.

Die klassische nicht-relativistische Formel hat schon 1913 Niels Bohr aufgestellt, die quantenmechanisch nicht-relativistische Formel wurde dann 1930, die unten gezeigte quantenmechanisch-relativistische Version 1932 von Hans Bethe aufgestellt.

Die Bethe-Bloch-Formel gilt nicht für einfallende Elektronen. Zum einen ist für diese der Energieverlust wegen ihrer Ununterscheidbarkeit mit den Hüllenelektronen des Materials anders. Zum anderen kommt bei Elektronen aufgrund ihrer geringen Masse ein bedeutender Energieverlust durch Bremsstrahlung hinzu. Der Energieverlust von Elektronen kann stattdessen mit Hilfe der Berger-Seltzer-Formel beschrieben werden.

Weitere Mechanismen, die zum Gesamt-Energieverlust schneller geladener schwerer Teilchen in Materie beitragen können, sind die nukleare Abbremsung (elastische Coulomb-Stöße mit den Atomkernen, siehe Bremsvermögen) und die Bremsstrahlung.

Die Formel

Bremsvermögen von Aluminium für Protonen als Funktion der Energie des Protons, mit Bethe-Formel ohne (rot) bzw. mit Korrekturen (blau)

Bewegen sich schnelle geladene Teilchen durch Materie, führen sie inelastische Stöße mit den Hüllenelektronen des Materials aus. Dies führt zur Anregung oder zur Ionisation der Atome. Dadurch erleidet das durchquerende Teilchen einen Energieverlust, der durch die folgende Formel näherungsweise angegeben wird. Ihre relativistische Form lautet:

$-{\frac {\mathrm {d} E}{\mathrm {d} x}}={\frac {4\pi nz^{2}}{m_{\rm {e}}c^{2}\beta ^{2}}}\cdot \left({\frac {e^{2}}{4\pi \varepsilon _{0}}}\right)^{2}\cdot \left[\ln \left({\frac {2m_{\rm {e}}c^{2}\beta ^{2}}{I\cdot (1-\beta ^{2})}}\right)-\beta ^{2}\right]\qquad \qquad$ (1)

wobei

$\beta$	= $v/c$
$v$	= momentane Geschwindigkeit des Teilchens
$c$	= Lichtgeschwindigkeit
$E$	= Energie des Teilchens
$x$	= Weglänge
$z$	= Ladungszahl des Teilchens ($z\cdot e$ = Ladung des Teilchens)
$\varepsilon _{0}$	= Elektrische Feldkonstante
$e$	= Elementarladung
$n$	= Elektronendichte des Materials
$m_{\rm {e}}$	= Masse des Elektrons
$I$	= mittleres Anregungspotential des Materials (s. u.)

Die Elektronendichte $n$ lässt sich dabei mit $n={\frac {Z\cdot \rho }{A\cdot u}}$ berechnen; $\rho$ ist die Dichte des Materials, $Z$ und $A$ Ordnungs- bzw. Massenzahl des Materials und $u$ die atomare Masseneinheit.

Im Bild rechts bedeuten die kleinen Kreise Messergebnisse von verschiedenen Arbeitsgruppen; die rote Kurve stellt die Bethe-Formel dar. Offenbar ist die Übereinstimmung von Bethes Theorie mit den Experimenten oberhalb von 0,5 MeV sehr gut, besonders wenn die Korrekturen (s. u.) hinzugefügt werden (blaue Kurve).

Für kleine Energien, d. h. kleine Teilchengeschwindigkeiten $(\beta \ll 1)$, reduziert sich die Bethe-Formel auf

$-{\frac {\mathrm {d} E}{\mathrm {d} x}}={\frac {4\pi nz^{2}}{m_{e}v^{2}}}\cdot \left({\frac {e^{2}}{4\pi \varepsilon _{0}}}\right)^{2}\cdot \ln \left({\frac {2m_{e}v^{2}}{I}}\right).$

Bei kleinen Energien ist die Bethe-Formel nur dann gültig, wenn diese noch hoch genug sind, dass das durchquerende Teilchen keine Hüllenelektronen mit sich führt. Anderenfalls wird seine effektive Ladung dadurch reduziert, und das Bremsvermögen ist kleiner. Es gibt für kleine Energien eine verfeinerte Theorie der elektronischen Abbremsung von Jens Lindhard, Morten Scharff und Hans E. Schiøtt (LSS-Theorie). Näherungsweise wird auch die Barkas-Formel für die effektive Ladungszahl benutzt:

$z_{\mathrm {eff} }=z\left(1-\exp \left(-125\beta z^{-{\frac {2}{3}}}\right)\right)$

Allgemein fällt der Energieverlust mit steigender Energie zunächst etwa mit $1/v^{2}$ ab und erreicht ein Minimum bei etwa Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle E = 3m_Tc^2} , wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle m_T} die Masse des Teilchens ist (also z. B. für Protonen etwa bei 3 GeV, was im Bild nicht mehr sichtbar ist). Da für viele in der Teilchenphysik relevante Strahlungsteilchen und Absorbermaterialien der Energieverlust in der Nähe des Minimums ungefähr den gleichen Wert hat, werden Teilchen mit einer Energie in der Nähe des Minimums häufig zusammengefasst und als MIPs (

Minimum Ionizing Particles

, dt. minimal ionisierende Teilchen) bezeichnet. Als Faustformel für den spezifischen Energieverlust der MIPs gilt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle -\frac{1}{\rho}\frac{\mathrm d E}{\mathrm d x}\approx 2 \frac{\mathrm{MeV}\mathrm{cm}^{2}}{\mathrm g\;}} .

Bei noch höherer Energie steigt der Energieverlust wieder an. Bei sehr hohen Energien müssen auch Teilchenreaktionen berücksichtigt werden, die zu Sekundärteilchen führen. Der Energieverlust kann daher in materialabhängiger Weise noch stärker ansteigen.

In der Strahlenbiologie nennt man die Energieabgabe ionisierender Teilchen gemäß der Bethe-Bloch-Gleichung den Linearen Energietransfer (Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle LET_{\infty}} ) und verwendet die Einheit Kiloelektronenvolt pro Mikrometer (keV/µm).

Das mittlere Anregungspotential

Im Gültigkeitsbereich der Bethe-Formel (1) wird das durchdrungene Material neben der Teilchendichte Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle n} nur durch einzige Konstante, das mittlere Anregungspotential Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle I} , beschrieben.

Felix Bloch hat 1933 gezeigt, dass das mittlere Anregungspotential der Atome im Mittel etwa

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle I = (10\,\mathrm{eV}) \cdot Z \qquad \qquad } (2)

beträgt, wo Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Z} die Ordnungszahl der Atome des Materials bedeutet. Setzt man diese Größe in Formel (1) oben ein, so führt das zu einer Gleichung, die oft als Bethe-Bloch-Gleichung bezeichnet wird. Es gibt aber genauere Tabellen von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle I} als Funktion von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Z} . Mit ihnen erhält man bessere Resultate als mit Formel (2).

Das mittlere Anregungspotential Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle I} von Elementen, dividiert durch die Ordnungszahl Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Z} , aufgetragen über der Ordnungszahl

Im Bild ist das mittlere Anregungspotential der verschiedenen Elemente gezeigt, das die Information über das jeweilige Atom enthält. Die Daten stammen aus dem genannten ICRU-Report. Den Spitzen und Tälern in der Darstellung („Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Z_2} -Oszillationen“, wobei Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle Z_2} die Ordnungszahl des Materials bedeutet) entsprechen niedrigere bzw. höhere Werte des Bremsvermögens; diese Oszillationen beruhen auf der Schalenstruktur der Atome. Wie das Bild zeigt, gilt Formel (2) nur näherungsweise.

Korrekturen

Die Bethe-Formel wurde von Bethe mit Hilfe der quantenmechanischen Störungstheorie abgeleitet, das Ergebnis ist daher dem Quadrat der Ladung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} proportional. Eine bessere Beschreibung erhält man, wenn man auch Abweichungen berücksichtigt, die höheren Potenzen von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z} entsprechen, und zwar den Barkas-Andersen-Effekt (proportional Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z^3} nach Walter H. Barkas und Hans Henrik Andersen) und die Bloch-Korrektur (proportional Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle z^4} ). Auch muss die Bewegung der Hüllenelektronen im Atom des Materials berücksichtigt werden („Schalenkorrektur“).

Diese Korrekturen sind beispielsweise in den Programmen PSTAR und ASTAR des National Institute of Standards and Technology (NIST), die das Bremsvermögen für Protonen bzw. Alphateilchen berechnen, eingebaut. Die Korrekturen sind groß bei niedrigen Energien und werden immer kleiner, je größer die Energie wird.

Zusätzlich kommt bei sehr hohen Energien noch Fermis Dichtekorrektur hinzu.

Literatur

• P. Sigmund: Particle Penetration and Radiation Effects, General Aspects and Stopping of Swift Point Charges (= Springer Series in Solid State Sciences. Vol. 151). Springer, Berlin/Heidelberg 2006, ISBN 978-3-540-72622-7.

Weblinks

• Durchgang geladener Teilchen durch Materie, inklusive Plot (engl.; PDF-Datei; 512 kB)
• Bremsvermögen für Protonen und Alphateilchen (engl.)
• Stopping Power Daten und Kurven (engl.)