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.
Die nebenstehende Skizze zeigt die experimentelle Anordnung, bestehend aus einer Lichtquelle <math>Q</math>, einer Blende <math>S</math>, an der das einfallende Licht gebeugt wird, und einem Beobachtungsschirm, auf dem die auftreffende Lichtintensität an <math>P</math> 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
Das kirchhoffsche Beugungsintegral, auch fresnel-kirchhoffsches Beugungsintegral genannt, lautet<ref>Helmut Rauch: Neutronenoptik. In: Ludwig Bergmann, Clemens Schaefer: Lehrbuch der Experimentalphysik. Bd. 3: Optik. 9. Auflage. de Gruyter, Berlin/New York 1993, ISBN 3-11-012973-6, Kapitel 11.4.</ref>
- <math>\psi_P = \frac{a_Q\,k_0}{2\pi\,\mathrm i} \int_\text{Blende}\mathrm dS\,f_S\,\frac{e^{\mathrm i\,k_0(d+d_1)}}{d\cdot d_1}\left[\frac{\cos\theta+\cos{\theta_1}}{2}\right].</math>
Dabei bezeichnen
- <math>a_Q</math> die Amplitude der Quelle,
- <math>k_0=2\pi/\lambda</math> den Betrag des Wellenvektors,
- <math>\lambda</math> die Wellenlänge des Lichtes,
- <math>\mathrm dS</math> ein infinitesimales Flächenelement der Blende,
- <math>f_S</math> die Blendenfunktion,
- <math>(\cos\theta+\cos\theta_1)/ 2</math> den Neigungsfaktor (auch als Inklinationsfaktor bezeichnet) und schließlich
- <math>\psi_P</math> die Amplitude im Punkt <math>P</math> auf dem Beobachtungsschirm.
Da die Abstände <math>d_1</math> und <math>d</math> in den meisten Anwendungen hinreichend senkrecht zur Blende sind, kann der Neigungsfaktor in diesen Fällen gleich Eins gesetzt werden. Dabei sind <math>\theta_1</math> bzw. <math>\theta</math> die Winkel zwischen den mit <math>d_1</math> bzw. <math>d</math> 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 <math>\theta = \pi</math>, und <math>\theta_1 = 0</math>. 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 <math>d_1</math> konstant. Dann kann <math>\frac{e^{ik_0d_1}}{d_1}</math> als konstanter Faktor gleich Eins gesetzt werden und es gilt:
- <math>\psi_{PEB} = \frac{a_Q\,k_0}{2\pi\,\mathrm i} \int_\text{Blende}\mathrm dS\,f_S\,\frac{e^{\mathrm i\,k_0d}}{d}\left[\frac{\cos\theta+\cos{\theta_1}}{2}\right].</math>
Die Intensität am Punkt <math>P</math> ergibt sich als Betragsquadrat von <math>\psi_P</math>
- <math>I(P)=| \psi_P |^2 = \frac{a_Q^2\,k_0^2}{4\pi^2 }\left|\int_\text{Blende}\mathrm dS\,f_S\,\frac{e^{\mathrm i\,k_0(d+d_1)}}{d\cdot d_1}\left[{\frac{\cos\theta+\cos{\theta_1}}{2}}\right]\right|^2.</math>
Neben der Intensität kann aber auch die Phase der gebeugten Wellenfront berechnet werden:
<math>\phi = \arctan\frac{\Im(\psi_P)}{\Re(\psi_P)}</math>
wobei <math>\Re</math> Realteil und <math>\Im</math> Imaginärteil bedeutet. Diese Phasenberechnung wird z. B. in der digitalen Holografie angewendet.
Fraunhofer- und Fresnel-Beugung
Die folgenden Ableitungen findet man in ähnlicher Form in verschiedenen Lehrbüchern, siehe z. B. Klein, Furtak.<ref>Miles V. Klein, Thomas E. Furtak: Optics. Second Edition, John Wiley and Sons, 1986, ISBN 0-471-87297-0, chapter 6 and 7.</ref>
Für die Lichtwege <math>d</math> und <math>d_1</math> gelten die geometrischen Zusammenhänge (siehe Skizze)
- <math>d=\sqrt{L^2+|\vec{s}+(-\vec{p})|^2}</math> und
- <math>d_1=\sqrt{L_1^2+|\vec{s}|^2}</math>.
Unter den Annahmen <math>L_1\gg|\vec{s}|=\sqrt{x^2+y^2}</math> und <math>L\gg|\vec{p}|=\sqrt{x'^2+y'^2}</math> können die Wurzeln durch eine Taylor-Entwicklung angenähert werden.
Diese Näherung entspricht gerade dem Fall, dass <math>\theta\approx\theta_1\approx 0</math>, d. h., für diese Betrachtungen kann der Neigungsfaktor näherungsweise gleich 1 gesetzt werden. Damit lautet das Beugungsintegral
- <math>\psi_P = \frac{a_Q\,k_0}{2\pi\,\mathrm i}\int_\text{Blende}\mathrm dS\,f_S\,\frac{e^{\mathrm i\,k_0(d+d_1)}}{d\cdot d_1}.</math>
Ferner kann wegen der Näherung im Nenner <math>d\cdot d_1\approx L_1\,L</math> gesetzt werden. Der Exponent enthält die für die Interferenz wesentliche Phaseninformation und darf nicht auf diese Weise vereinfacht werden. Daraus folgt
- <math>\psi_P = \frac{a_Q\,k_0}{2\pi\,\mathrm i}\frac{1}{L_1\,L} \int_\text{Blende}\mathrm dS\,f_S\,e^{\mathrm i\,k_0(d+d_1)}.</math>
Die Näherung für die Ausdrücke <math>d</math> und <math>d_1</math>, explizit ausgeführt bis zur 2. Ordnung, ergibt
- <math>d = \sqrt{L^2+|\vec{s}-\vec{p}|^2}\approx L\left( 1+\frac{|\vec{s}-\vec{p}|^2}{ 2L^2}+ \ldots \right)= L\left( 1+\frac{s^2+p^2-2\vec{s}\cdot\vec{p}}{2L^2}+ \ldots \right)</math>
sowie
- <math>d_1=\sqrt{L_1^2+s^2}\approx L_1\left( 1+\frac{s^2}{2 L_1^2}+ \ldots \right).</math>
Ausgedrückt durch die Koordinaten <math>(x,\,y)</math> und <math>(x',\,y')</math> ergibt das
- <math>d \approx L\left( 1 + \frac{x^2+y^2+x'^2+y'^2-2(x\,x'+y\,y')}{2L^2}+ \ldots \right)</math>
und
- <math>d_1 \approx L_1\left( 1 + \frac{x^2+y^2}{2 L_1^2} + \ldots \right).</math>
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 <math>L</math> 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 <math>x</math> und <math>y</math> sind, das heißt
- <math>d\approx L\left( 1+\frac{x'^2+y'^2-2(x\,x'+y\,y')}{ 2L^2} \right)= L\left( 1+\frac{|\vec{p}|^2-2(x\,x'+y\,y')}{2L^2} \right),</math>
- <math>d_1\approx L_1.</math>
In diesem Fall vereinfacht sich das Beugungsintegral zu
- <math>\psi_P\approx\frac{a_Q\,k_0}{2\pi\,\mathrm i}\frac{e^{\mathrm i\,k_0\left(L_1+L+\frac{p^2}{ 2L}\right)}}{L_1\,L}\iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{-\mathrm i\,k_0\frac{(x\,x'+y\,y')}{ L}}.</math>
Definiert man einen neuen Wellenvektor <math>\vec{K}= \frac{k_0}{L}\vec{p}</math>, so ergibt sich für das Integral
- <math>\iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{-\mathrm i\,k_0\frac{(x\,x'+y\,y')}{L}} = \iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{-\mathrm i\,\frac{k_0}{L}\vec{p}\cdot\vec{s}}= \iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{-\mathrm i\,\vec{K}\cdot\vec{s}}</math>.
Dies ist gerade die Fourier-Transformierte der Blendenfunktion <math>f_S</math>.
Fresnel-Näherung
{{#invoke:Vorlage:Anker|f |errCat=Wikipedia:Vorlagenfehler/Vorlage:Anker |errHide=1}} 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 <math>x</math> und <math>y</math> ergibt sich
- <math>d\approx L\left( 1+\frac{x^2+y^2+x'^2+y'^2-2(x\,x'+y\,y')}{ 2L^2}+ \ldots \right),</math>
- <math>d_1\approx L_1\left( 1+\frac{x^2+y^2}{ 2 L_1^2}+ \ldots \right).</math>
In diesem Fall lautet das Beugungsintegral
- <math>\psi_P\approx\frac{a_Q\,k_0}{ 2\pi\,\mathrm i}\frac{ e^{\mathrm i\,k_0\left(L+L_1+\frac{p^2}{ 2L}\right)}
}{ L\,L_1} \iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{\mathrm i\,k_0\left(\frac{x^2+y^2-2(x\,x'+y\,y')}{ 2L}+\frac{x^2+y^2}{ 2 L_1}\right)}.</math>
Einführung von <math>L'</math> mit <math>\frac{1}{ L'}=\frac{1}{L}+\frac{1}{ L_1}</math> und <math>\vec{K}= \frac{k_0}{L}\vec{p}</math> ergibt dann das Beugungsintegral in Nahfeld-Näherung
- <math>\psi_P\approx\frac{a_Q\,k_0}{2\pi\,\mathrm i}\frac{ e^{\mathrm i\,k_0\left(L+L_1+\frac{p^2}{ 2L}\right)}
}{ L\,L_1} \iint_\mathrm{Blende}\mathrm dx\,\mathrm dy\,f_S(x,y)\,e^{i\,k_0\frac{x^2+y^2}{ 2L'}}e^{-\mathrm i\,\vec{K}\cdot\vec{s}}.</math>
Die Fresnel-Näherung gilt, wenn folgende Bedingung erfüllt ist:<ref>Miles V. Klein, Thomas E. Furtak: Optics. Second Edition, John Wiley and Sons, 1986, ISBN 0-471-87297-0, S. 413</ref>
- <math>d \gg \frac{((x-x')^2+(y-y')^2)^{2/3}}{(8\lambda)^{1/3}}</math>
Dazu ein praktisches Beispiel: Für eine Blendengröße von 1 cm (<math>x-x' = 1\,\text{cm}</math>, <math>y-y' = 1\,\text{cm}</math>) und <math>\lambda = 500\,\text{nm}</math> 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 Abstände muss das Beugungsintegral ohne Näherung verwendet werden.
Heuristische Herleitung
Aus der Quelle <math>Q</math> mit Amplitude <math>a_Q</math> bei <math>\vec{r}_Q </math> tritt die Kugelwelle <math>\psi_Q</math>, deren Amplitude reziprok mit der Entfernung (<math>1/|\vec{r}|</math>) abnimmt. Wellenvektor <math>k</math> mal Abstand <math>|\vec{r}|</math> gibt die Phasenverschiebung der Welle am Ort <math>\vec{r}</math>, Kreisfrequenz <math>\omega</math> mal Zeit <math>t</math> die Phasenverschiebung zur Zeit <math>t</math>. Die Welle ist beschrieben durch die Phase am Ort <math>\vec{r}</math> zur Zeit <math>t</math>:
- <math>\psi_Q(\vec{r},t)=a_Q\, \frac{e^{\mathrm i(k\,|\vec{r}|-\omega\,t)}}{ |\vec{r}|}.</math>
Am Punkt <math>S</math> bei <math>\vec{r}_S </math> trifft die Welle im Abstand <math>d_1</math> auf die Blende. Es sei <math>\psi_1</math> die Feldverteilung der Welle am Punkt <math>S</math>:
- <math>\psi_1(t)=\psi_Q(d_1,t)=a_Q\,\frac{e^{\mathrm i(k\,d_1-\omega\,t)}}{ d_1}.</math>
Nach dem huygensschen Prinzip ist der Punkt <math>S</math> Ausgangspunkt einer Elementarwelle, der Sekundärwelle <math>\psi_S </math>:
- <math>\psi_S(\vec{r},t)=a_S\, \frac{e^{i(k\,|\vec{r}|-\omega\,t)}}{ |\vec{r}|}.</math>
Die Amplitude von <math>\psi_S </math> ist proportional zur Quellen-Amplitude <math>a_Q</math> und zur Blendenfunktion <math>f_S</math>. Die Blendenfunktion gibt die Durchlässigkeit der Blende an. Im einfachsten Fall ist <math>f_S=1</math>, wenn die Blende geöffnet ist, und <math>f_S=0</math>, wenn die Blende geschlossen ist. <math>\mathrm dS</math> ist das infinitesimale Flächenelement der Blendenöffnung am Punkt <math>S</math>:
- <math>a_S(t)\sim f_S\,\mathrm dS\,\psi_1(t).</math>
Die Sekundärwelle <math>\psi_S </math> erzeugt im Punkt <math>P</math> bei <math>\vec{r}_P</math> auf dem Schirm die Wellenintensität <math>\mathrm d\psi_P(t)</math>. Sie ist infinitesimal, da nur der Beitrag von <math>\mathrm dS</math> und nicht aller anderen Punkte auf der Blende betrachtet wird:
- <math>\mathrm d\psi_P(t)=\psi_S(\vec{d};t)=a_S\,\frac{e^{\mathrm i(k\,d-\omega\,t)} }{ d}.</math>
Die Zeitabhängigkeit <math>e^{-\mathrm i\,\omega\,t}</math> kann vernachlässigt werden, da sie später beim Rechnen mit Intensitäten ohnehin durch die zeitliche Mittelung verschwindet. Durch Einsetzen erhält man:
- <math>\mathrm d\psi_P \sim f_S\,\mathrm dS\,\psi_1(t)\,\frac{e^{\mathrm i\,k\,d}}{d} = f_S\,\mathrm dS\,a_Q\,\frac{e^{\mathrm i\,k\,d_1}}{d_1}\,\frac{e^{\mathrm i\,k\,d}}{d} = f_S\,\mathrm dS\,a_Q\,\frac{e^{\mathrm i\,k\,(d_1+d)}}{d_1\,d}.</math>
Von jedem Punkt auf der Blende geht eine Sekundärwelle aus. Die Intensität <math>\psi_P</math> im Beobachtungspunkt <math>P</math> wird durch die Überlagerung aller Einzelbeiträge erzeugt:
- <math>\psi_P\sim a_Q\int_\text{Blende}\mathrm dS\,f_S\, \frac{e^{\mathrm i\,k(d+d_1)}}{d \cdot d_1}.</math>
Diese Gleichung erinnert bereits stark an das oben gegebene Beugungsintegral. Mit dem Proportionalitätsfaktor <math>k/(2\pi\,\mathrm i)</math> ergibt sich (Neigungswinkel vernachlässigt):<ref>Eine exakte Herleitung unter Berücksichtigung des Neigungsfaktors findet man z. B. in Miles V. Klein, Thomas E. Furtak: Optics. Second Edition, John Wiley and Sons, 1986, ISBN 0-471-87297-0, Appendix A, S. 649–652.</ref>
- <math>\psi_P= a_Q\,\frac{k}{2\pi\,\mathrm i}\,\int_\text{Blende}\mathrm dS\,f_S\,\frac{e^{\mathrm i\,k(d+d_1)}}{ d\cdot d_1}.</math>
Einzelnachweise
<references />