Beugungsintegral

Das Beugungsintegral ermöglicht es, in der Optik die Beugung von Licht durch eine beliebig geformte Blende zu berechnen. Speziell wird dabei die an einem Punkt des Beobachtungsschirms auftreffende Intensität des Lichtes berechnet, ausgehend von einer einfallenden Elementarwelle und der Blendenfunktion, welche die Lichtdurchlässigkeit der Blende beschreibt.

Zwei Grenzfälle des Beugungsintegrals sind die Näherungen für das Fernfeld (Fraunhofer-Beugung) und für das Nahfeld (Fresnel-Beugung). Siehe dazu die entsprechenden Teilabschnitte.

Versuchsaufbau zur Beugung von Licht an einer Blende

Die nebenstehende Skizze zeigt die experimentelle Anordnung, bestehend aus einer Lichtquelle , einer Blende , an der das einfallende Licht gebeugt wird, und einem Beobachtungsschirm, auf dem die auftreffende Lichtintensität an untersucht wird. Die Form und die Eigenschaften der Blende bestimmen dabei, wie die Intensitätsverteilung auf dem Beobachtungsschirm aussieht.

Hat die Blende z.B. die Form eines Doppelspalts, so ergibt sich als Intensitätsverteilung das bekannte Interferenzmuster. Weitere Anwendungen des Beugungsintegrals sind z.B. Beugungsscheibchen und Klotoide.

Das Kirchhoffsche Beugungsintegral

Skizze zur Fraunhofer-/Fresnel-Näherung des Beugungsintegrals

Das Kirchhoffsche Beugungsintegral, auch Fresnel-Kirchhoffsches Beugungsintegral genannt, lautet

$\psi _{P}={\frac {a_{Q}\,k_{0}}{2\pi \,\mathrm {i} }}\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,{\frac {e^{\mathrm {i} \,k_{0}(d+d_{1})}}{d\cdot d_{1}}}\left[{\frac {\cos \theta +\cos {\theta _{1}}}{2}}\right].$

Dabei bezeichnen

$a_{Q}$ die Amplitude der Quelle,
$k_{0}=2\pi /\lambda$ den Betrag des Wellenvektors,
$\lambda$ die Wellenlänge des Lichtes,
$\mathrm {d} S$ ein infinitesimales Flächenelement der Blende,
$f_{S}$ die Blendenfunktion,
$(\cos \theta +\cos \theta _{1})/2$ den Neigungsfaktor und schließlich
$\psi _{P}$ die Amplitude im Punkt auf dem Beobachtungsschirm.

Da die Abstände $d_{1}$ und in den meisten Anwendungen hinreichend senkrecht zur Blende sind, kann der Neigungsfaktor in diesen Fällen gleich Eins gesetzt werden. Dabei sind $\theta _{1}$ bzw. $\theta$ die Winkel zwischen den mit $d_{1}$ bzw. gekennzeichneten Linien und einem Lot auf die Blendenebene im Schnittpunkt der Linien.

Die Intensität am Punkt ergibt sich als Betragsquadrat von $\psi _{P}$

$I(P)=|\psi _{P}|^{2}={\frac {a_{Q}^{2}\,k_{0}^{2}}{4\pi ^{2}}}\left|\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,{e^{\mathrm {i} \,k_{0}(d+d_{1})} \over d\cdot d_{1}}\left[{\frac {\cos \theta +\cos {\theta _{1}}}{2}}\right]\right|^{2}.$

Fraunhofer- und Fresnel-Beugung

Prinzip der Fresnelbeugung erläutert anhand einer Schlitzblende

Prinzip der Fraunhoferbeugung erläutert anhand einer Schlitzblende

Prinzip der Fraunhoferbeugung erläutert anhand eines Linsensystems und einer Schlitzblende

Für die Lichtwege und $d_{1}$ gelten die geometrischen Zusammenhänge (siehe Skizze)

$d={\sqrt {L^{2}+|{\vec {s}}+(-{\vec {p}})|^{2}}}$ und

$d_{1}={\sqrt {L_{1}^{2}+|{\vec {s}}|^{2}}}$ .

Unter den Annahmen $L_{1}\gg |{\vec {s}}|={\sqrt {x^{2}+y^{2}}}$ und $L\gg |{\vec {p}}|={\sqrt {x'^{2}+y'^{2}}}$ können die Wurzeln durch eine Taylor-Entwicklung angenähert werden.

Diese Näherung entspricht gerade dem Fall, dass $\theta \approx \theta _{1}\approx 0$ , d.h., für diese Betrachtungen kann der Neigungsfaktor näherungsweise gleich 1 gesetzt werden. Damit lautet das Beugungsintegral

$\psi _{P}={\frac {a_{Q}\,k_{0}}{2\pi \,\mathrm {i} }}\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,{e^{\mathrm {i} \,k_{0}(d+d_{1})} \over d\cdot d_{1}}.$

Ferner kann wegen der Näherung im Nenner $d\cdot d_{1}\approx L_{1}\,L$ gesetzt werden. Der Exponent enthält die für die Interferenz wesentliche Phaseninformation und darf nicht auf diese Weise vereinfacht werden. Daraus folgt

$\psi _{P}={\frac {a_{Q}\,k_{0}}{2\pi \,\mathrm {i} }}{\frac {1}{L_{1}\,L}}\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,e^{\mathrm {i} \,k_{0}(d+d_{1})}.$

Die Näherung für die Ausdrücke und $d_{1}$ , explizit ausgeführt bis zur 2. Ordnung, ergibt

$d={\sqrt {L^{2}+|{\vec {s}}-{\vec {p}}|^{2}}}\approx L\left(1+{|{\vec {s}}-{\vec {p}}|^{2} \over 2L^{2}}+\ldots \right)=L\left(1+{s^{2}+p^{2}-2{\vec {s}}\cdot {\vec {p}} \over 2L^{2}}+\ldots \right)$

sowie

$d_{1}={\sqrt {L_{1}^{2}+s^{2}}}\approx L_{1}\left(1+{s^{2} \over 2L_{1}^{2}}+\ldots \right).$

Ausgedrückt durch die Koordinaten $(x,\,y)$ und $(x',\,y')$ ergibt das

$d\approx L\left(1+{\frac {x^{2}+y^{2}+x'^{2}+y'^{2}-2(x\,x'+y\,y')}{2L^{2}}}+\ldots \right)$

und

$d_{1}\approx L_{1}\left(1+{\frac {x^{2}+y^{2}}{2L_{1}^{2}}}+\ldots \right).$

Fraunhofer-Näherung

Die Fraunhofer-Näherung entspricht einer Fernfeld-Näherung, das heißt, dass nicht nur die Blendenöffnung als klein, sondern auch die Entfernung des Beobachtungsschirms als groß angenommen werden. Als Beugungsintegral ergibt sich dabei im Wesentlichen gerade die Fourier-Transformierte der Blendenfunktion. Deshalb spricht man im Rahmen der Fraunhofer-Beugung auch von der Fourier-Optik.

Entsprechend diesen Annahmen werden nur Terme berücksichtigt, die linear in und sind, das heißt

$d\approx L\left(1+{x'^{2}+y'^{2}-2(x\,x'+y\,y') \over 2L^{2}}\right)=L\left(1+{|{\vec {p}}|^{2}-2(x\,x'+y\,y') \over 2L^{2}}\right),$

$d_{1}\approx L_{1}.$

In diesem Fall vereinfacht sich das Beugungsintegral zu

$\psi _{P}\approx {a_{Q}\,k_{0} \over 2\pi \,\mathrm {i} }{e^{\mathrm {i} \,k_{0}(L_{1}+L+{p^{2} \over 2L})} \over L_{1}\,L}\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{-\mathrm {i} \,k_{0}{(x\,x'+y\,y') \over L}}.$

Definiert man einen neuen Wellenvektor ${\vec {K}}={k_{0} \over L}{\vec {p}}$ , so ergibt sich für das Integral

$\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{-\mathrm {i} \,k_{0}{(x\,x'+y\,y') \over L}}=\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{-\mathrm {i} \,{k_{0} \over L}{\vec {p}}\cdot {\vec {s}}}=\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{-\mathrm {i} \,{\vec {K}}\cdot {\vec {s}}}$ .

Dies ist gerade die Fourier-Transformierte der Blendenfunktion $f_{S}\,(x,y)$ .

Fresnel-Näherung

Die Fresnel-Näherung entspricht einer Nahfeld-Näherung. In ihr werden auch quadratische Terme im Exponenten berücksichtigt. Das in die Form einer Fourier-Transformierten gebrachte Beugungsintegral ist dann durch einen zusätzlichen Term im Allgemeinen nicht mehr analytisch, sondern nur numerisch lösbar.

Unter Berücksichtigung quadratischer Terme in und ergibt sich

$d\approx L\left(1+{x^{2}+y^{2}+x'^{2}+y'^{2}-2(x\,x'+y\,y') \over 2L^{2}}+\ldots \right),$

$d_{1}\approx L_{1}\left(1+{x^{2}+y^{2} \over 2L_{1}^{2}}+\ldots \right).$

In diesem Fall lautet das Beugungsintegral

$\psi _{P}\approx {a_{Q}\,k_{0} \over 2\pi \,\mathrm {i} }{e^{\mathrm {i} \,k_{0}(L+L_{1}+{p^{2} \over 2L})} \over L\,L_{1}}\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{\mathrm {i} \,k_{0}({x^{2}+y^{2}-2(x\,x'+y\,y') \over 2L}+{x^{2}+y^{2} \over 2L_{1}})}.$

Einführung von mit ${1 \over L'}={1 \over L}+{1 \over L_{1}}$ und ${\vec {K}}={k_{0} \over L}{\vec {p}}$ ergibt dann das Beugungsintegral in Nahfeld-Näherung

$\psi _{P}\approx {a_{Q}\,k_{0} \over 2\pi \,\mathrm {i} }{e^{\mathrm {i} \,k_{0}(L+L_{1}+{p^{2} \over 2L})} \over L\,L_{1}}\iint _{\mathrm {Blende} }\mathrm {d} x\,\mathrm {d} y\,f_{S}(x,y)\,e^{i\,k_{0}{x^{2}+y^{2} \over 2L'}}e^{-\mathrm {i} \,{\vec {K}}\cdot {\vec {s}}}.$

Heuristische Herleitung

Bezeichnungen

Aus der Quelle mit Amplitude $a_{Q}$ bei ${\vec {r}}_{Q}$ tritt die Kugelwelle $\psi _{Q}$ , deren Amplitude reziprok mit der Entfernung ( $1/|{\vec {r}}|$ ) abnimmt. Wellenvektor mal Abstand $|{\vec {r}}|$ gibt die Phasenverschiebung der Welle am Ort ${\vec {r}}$ , Kreisfrequenz $\omega$ mal Zeit die Phasenverschiebung zur Zeit . Die Welle ist beschrieben durch die Phase am Ort ${\vec {r}}$ zur Zeit :

$\psi _{Q}({\vec {r}},t)=a_{Q}\,{e^{\mathrm {i} (k\,|{\vec {r}}|-\omega \,t)} \over |{\vec {r}}|}.$

Am Punkt bei ${\vec {r}}_{S}$ trifft die Welle im Abstand $d_{1}$ auf die Blende. Es sei $\psi _{1}$ die Feldverteilung der Welle am Punkt .

$\psi _{1}(t)=\psi _{Q}(d_{1},t)=a_{Q}\,{e^{\mathrm {i} (k\,d_{1}-\omega \,t)} \over d_{1}}$

Nach dem Huygensschen Prinzip ist der Punkt Ausgangspunkt einer Elementarwelle, der Sekundärwelle $\psi _{S}$ .

$\psi _{S}({\vec {r}},t)=a_{S}\,{e^{i(k\,|{\vec {r}}|-\omega \,t)} \over |{\vec {r}}|}.$

Die Amplitude von $\psi _{S}$ ist proportional zur Quellen-Amplitude $a_{Q}$ und zur Blendenfunktion $f_{S}$ . Die Blendenfunktion gibt die Durchlässigkeit der Blende an. Im einfachsten Fall ist $f_{S}=1$ , wenn die Blende geöffnet ist, und $f_{S}=0$ , wenn die Blende geschlossen ist. > $\mathrm {d} S$ ist das infinitesimale Flächenelement der Blendenöffnung am Punkt .

$a_{S}(t)\sim f_{S}\,\mathrm {d} S\,\psi _{1}(t)$

Die Sekundärwelle $\psi _{S}$ erzeugt im Punkt bei ${\vec {r}}_{P}$ auf dem Schirm die Wellenintensität $\mathrm {d} \psi _{P}(t)$ . Sie ist infinitesimal, da nur der Beitrag von $\mathrm {d} S$ und nicht aller anderen Punkte auf der Blende betrachtet wird.

$\mathrm {d} \psi _{P}(t)=\psi _{S}({\vec {d}};t)=a_{S}\,{e^{\mathrm {i} (k\,d-\omega \,t)} \over d}$

Die Zeitabhängigkeit $e^{-\mathrm {i} \,\omega \,t}$ kann vernachlässigt werden, da sie später beim Rechnen mit Intensitäten ohnehin durch die zeitliche Mittelung verschwindet.
Durch Einsetzen erhält man:

$\mathrm {d} \psi _{P}\sim f_{S}\,\mathrm {d} S\,\psi _{1}(t)\,{\frac {e^{\mathrm {i} \,k\,d}}{d}}=f_{S}\,\mathrm {d} S\,a_{Q}\,{\frac {e^{\mathrm {i} \,k\,d_{1}}}{d_{1}}}\,{\frac {e^{\mathrm {i} \,k\,d}}{d}}=f_{S}\,\mathrm {d} S\,a_{Q}\,{\frac {e^{\mathrm {i} \,k\,(d_{1}+d)}}{d_{1}\,d}}$

Von jedem Punkt auf der Blende geht eine Sekundärwelle aus. Die Intensität $\psi _{P}$ im Beobachtungspunkt wird durch die Überlagerung aller Einzelbeiträge erzeugt:

$\psi _{P}\sim a_{Q}\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,{\frac {e^{\mathrm {i} \,k(d+d_{1})}}{d\cdot d_{1}}}.$

Diese Gleichung erinnert bereits stark an das oben gegebene Beugungsintegral. Mit dem Proportionalitätsfaktor $k/(2\pi \,\mathrm {i} )$ ergibt sich (Neigungswinkel vernachlässigt):

$\psi _{P}=a_{Q}\,{\frac {k}{2\pi \,\mathrm {i} }}\,\int _{\text{Blende}}\mathrm {d} S\,f_{S}\,{e^{\mathrm {i} \,k(d+d_{1})} \over d\cdot d_{1}}.$

Basierend auf einem Artikel in:

Wikipedia.de