Beugungsintegral

Das Beugungsintegral ermöglicht es, in der Optik die Beugung von Licht durch eine beliebig geformte Blende zu berechnen. Speziell wird dabei meistens 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. Es kann aber auch die Phase der gebeugten Wellenfront berechnet werden.

Zwei Spezialfä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 $Q$, einer Blende $S$, an der das einfallende Licht gebeugt wird, und einem Beobachtungsschirm, auf dem die auftreffende Lichtintensität an $P$ 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 (auch als Inklinationsfaktor bezeichnet) und schließlich
• $\psi _{P}$ die Amplitude im Punkt $P$ auf dem Beobachtungsschirm.

Da die Abstände $d_{1}$ und $d$ 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. $d$ gekennzeichneten Linien und einem Lot auf die Blendenebene im Schnittpunkt der Linien.

Der Neigungsfaktor korrigiert eine Inkonsistenz des huygensschen Prinzips. Nach diesem Prinzip kann jeder Punkt einer Wellenfront als Ausgangspunkt von Elementarwellen aufgefasst werden, deren Überlagerung dann die Wellenfront an einem anderen Ort ergibt. Allerdings würde dies auch bedeuten, dass die Elementarwellen sich in alle Richtungen ausbreiten, also auch entgegen der Ausbreitungsrichtung, zurück in Richtung auf die Lichtquelle. In Richtung auf die Lichtquelle beträgt der Winkel $\theta =\pi$, und $\theta _{1}=0$. Dann ergibt sich für den Neigungsfaktor 0 und die Intensität in Richtung auf die Quelle ist ebenfalls 0 – so wie es experimentell beobachtet wird.

Im Spezialfall einer ebenen Welle, die senkrecht auf die Blende fällt, ist der Abstand $d_{1}$ konstant. Dann kann ${\frac {e^{ik_{0}d_{1}}}{d_{1}}}$ als konstanter Faktor ohne Beschränkung der Allgemeinheit gleich Eins gesetzt werden und es gilt:

$\psi _{PEB}={\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}}\left[{\frac {\cos \theta +\cos {\theta _{1}}}{2}}\right].$

Die Intensität am Punkt $P$ 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}.$

Neben der Intensität kann aber auch die Phase der gebeugten Wellenfront berechnet werden:

$\phi =\arctan {\frac {\Im (\psi _{P})}{\Re (\psi _{P})}}$

wobei $\Re$ Realteil und $\Im$ Imaginärteil bedeutet. Diese Phasenberechnung wird z. B. in der digitalen Holografie angewendet.

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

Die folgenden Ableitungen findet man in ähnlicher Form in verschiedenen Lehrbüchern, siehe z. B.

Für die Lichtwege $d$ 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 $d$ 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 $L$ 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 $x$ und $y$ 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 $x$ und $y$ 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 $L'$ 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}}}.$

Die Fresnel-Näherung gilt, wenn folgende Bedingung erfüllt ist:

$d\gg {\sqrt[{3}]{\frac {((x-\,x')^{2}+(y-\,y')^{2})^{2}}{8\lambda }}}$

Dazu ein praktisches Beispiel: Für eine Blendengröße von 1 cm ($(x-\,x')=1\,\mathrm {cm}$; $(y-\,y')=1\,\mathrm {cm}$) und $\lambda =500\,\mathrm {nm}$ ergibt sich ein Mindestabstand von 22 cm. Das Attribut Nahfeld bedeutet also keineswegs, dass man die Fresnel-Näherung schon direkt hinter der Blende anwenden darf, sondern dass ein Mindestabstand eingehalten werden muss. Für noch kürzere Abstande muss das Beugungsintegral ohne Näherung verwendet werden.

Heuristische Herleitung

Bezeichnungen

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

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

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

$\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 $S$ 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 $S$.

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

Die Sekundärwelle $\psi _{S}$ erzeugt im Punkt $P$ 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 $P$ 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}}.$