Bragg-Spiegel
# Bragg-Spiegel

Ein Bragg-Spiegel (abgekürzt DBR von distributed Bragg reflector) bezeichnet einen effizienten Reflektor, der in Lichtleitern oder in optischen Resonatoren eingesetzt wird. Er besteht aus alternierenden, dünnen Schichten unterschiedlicher Brechungsindizes. Meist bestehen die Schichten aus Dielektrika. Darum verwendet man bei solch einem Reflektor auch den Begriff dielektrischer Spiegel. An jeder Grenzschicht wird ein Teil der elektromagnetischen Welle des Lichtes gemäß den fresnelschen Formeln reflektiert. Wenn die Wellenlänge nahe dem Vierfachen der optischen Weglänge der Schichten liegt, dann interferieren die reflektierten Strahlen konstruktiv und es entsteht ein hochqualitativer Reflektor. Der Bereich, in dem die Reflexion sehr hoch ist, heißt Stoppband. Licht, dessen Wellenlänge innerhalb des Stoppbands liegt, kann sich in der Struktur nicht ausbreiten.

## Charakteristika

Bragg-Spiegel bestehen aus alternierenden, dielektrischen dünnen Schichten mit niedrigem und hohem Brechungsindex. Die maximale Reflektivität für eine Wellenlänge wird erreicht, wenn alle Schichten eine optische Dicke von genau einem Viertel der Wellenlänge aufweisen. In der Skizze rechts sind 4 Schichten eines Bragg-Spiegels illustriert. Trifft Licht senkrecht auf den Bragg-Spiegel, so kommt es an den Grenzflächen von niedrigem zu hohem Brechungsindex ($n_{L}$,$n_{H}$) zu einer Phasenverschiebung des elektrischen Feldvektors des Lichts von einer halben Wellenlänge ${\frac {\lambda }{2}}$. Bei den Übergängen von hohem zu niedrigem Brechungsindex ist dies jedoch nicht der Fall. Um konstruktive Interferenz des Bragg-Spiegels zu erzeugen, muss die gesamte Phasendifferenz jedoch ein Vielfaches der Wellenlänge des einfallenden Lichtes sein. Um konstruktive Interferenz an allen Grenzschichten zu erreichen, muss also die optische Weglänge eines jeden Dünnfilmes ${\frac {\lambda }{2}}$ sein. Die Bedingung für maximale Reflektivität des Bragg-Spiegels kann nun wie folgt ermittelt werden. Für konstruktive Interferenz muss die gesamte Phasenverschiebung ein ganzzahliges Vielfaches $m$ der Wellenlänge des einfallenden Lichtes sein.

    $2n_{L}d_{L}+{\frac {\lambda }{2}}=2n_{H}d_{H}+{\frac {\lambda }{2}}=m\lambda$

Es ergibt sich somit folgende Bedingung für die optische Weglänge der Dünnschichten des Bragg-Spiegels.

    $n_{H}d_{H}=n_{L}d_{L}={\frac {(2m-1)\lambda }{4}}$

Ein Bragg-Spiegel zeigt daher konstruktive Interferenz bei mehreren Wellenlängen. Es gibt also mehrere Wellenlängenbereiche konstruktiver Interferenz bei denen ein Maximum der Reflexion auftritt. Die Wellenlänge, für die maximale Reflektivität bei $m=1$ erfüllt ist, wird als $\lambda _{0}$ bezeichnet und liegt in der Mitte des sogenannten Stoppbands eines Bragg-Spiegels. Im Bild rechts ist das berechnete Reflexionsspektrum eines Bragg-Spiegels im Bereich des Stoppbands und $\lambda _{0}$ abgebildet. Die Reflektivität für diese Wellenlänge zu:

    $R=\left[{\frac {n_{o}(n_{2})^{2N}-n_{s}(n_{1})^{2N}}{n_{o}(n_{2})^{2N}+n_{s}(n_{1})^{2N}}}\right]^{2},$

wobei $n_{o}$ der Brechungsindex des Umgebungsmediums ist, $n_{1}$ und $n_{2}$ die Brechungsindizes der beiden Materialien und $n_{s}$ der Brechungsindex des Substrats. $N$ ist die Anzahl der Schichtpaare. Unter der Voraussetzung, dass beide Materialien unterschiedliche Brechungsindizes haben, ergibt $\lim \limits _{N\to \infty }R=1$. Es ist also möglich eine beliebig hohe Reflektivität zu erreichen, wenn nur genug Schichtpaare verwendet werden.

Die Frequenz-Breite $\Delta f_{0}$ des Stoppbands berechnet sich wie folgt:

    ${\frac {\Delta f}{f_{0}}}={\frac {4}{\pi }}\arcsin \left\vert {\frac {n_{1}-n_{2}}{n_{1}+n_{2}}}\right\vert$.

## Physikalische Wirkungsweise

### Reflexionsgrad

In diesem Abschnitt wird die Berechnung der Reflektivität $R$ eines Bragg-Spiegels erläutert. Der Leser kann sich dabei am Bild rechts orientieren, wo ein Dünnfilm auf einem Substrat skizziert ist. An der Grenzfläche zwischen Luft und Dünnfilm kommt es zur Reflexion einer elektromagnetischen Welle. Die einfallende Welle und reflektierte Wellenlänge haben die Wellenvektoren ${\vec {k_{0}}}$ und ${\vec {k'_{0}}}$. Weiters hat die an der ersten Grenzfläche transmittierte Welle den Wellenvektor ${\vec {k_{1}}}$ und die an der zweiten Grenzfläche reflektierte Welle wird durch ${\vec {k'_{1}}}$ beschrieben. Letztendlich hat die Welle, die sich ins Substrat ausbreitet, den Wellenvektor ${\vec {k_{2}}}$. Die Vektoren der magnetischen Feldstärke ${\vec {H}}$ sowie der elektrischen Feldstärke ${\vec {E}}$ sind in analoger Weise beschriftet. Die Vektoren der magnetischen Feldstärke zeigen in die Bildebene (markiert durch ein Kreuz) bei reflektierten Wellen und aus der Bildebene für einlaufende Wellen (markiert durch einen Punkt). Am ersten Phasenübergang muss die folgende Bedingung für die entsprechenden Amplituden der elektrischen Feldstärke-Vektoren gelten. Dies folgt aus den fresnelschen Gleichungen, die besagen, dass die Tangentialkomponenten der elektrischen Feldvektoren an einer Grenzfläche stetig sein müssen.

    $E_{0}+E_{0}'=E_{1}+E_{1}'$

Gleiches gilt für die Amplituden der magnetischen Feldstärke. Da jedoch bei der Reflexion die Orientierung der Vektoren umgekehrt wird, werden die reflektierten Amplituden abgezogen.

    $H_{0}-H_{0}'=H_{1}-H_{1}'$

Die Amplituden der magnetischen Feldstärke können durch die entsprechenden Amplituden der elektrischen Feldstärke ausgedrückt werden. Man bedient sich dabei der Beziehung $H={\frac {n}{\mu c}}E$. Der Brechungsindex wird durch $n$ repräsentiert, während $c_{0}$ die Lichtgeschwindigkeit im Vakuum darstellt und die magnetische Permeabilität mit $\mu$ gekennzeichnet ist. Unter der Annahme, dass die relative magnetische Permeabilität für alle beteiligten Materialien annähernd 1 ist, erhält man folgende Gleichungen, da man $c_{0}$ und $\mu$ kürzen kann.

    $n_{0}(E_{0}-E_{0}')=n_{1}(E_{1}-E_{1}')$

Somit sind die Stetigkeitsbedingungen an der ersten Grenzfläche durch die Amplituden der elektrischen Feldstärken mit zwei Gleichungen ausgedrückt. Selbiges kann auch für die zweite Grenzfläche gemacht werden. Dazu müssen die Amplituden $E_{1}$ und $E'_{1}$ jedoch mit einem zusätzlichen Phasenterm angeschrieben werden. Beide Amplituden ohne Phasenterm entsprechen den Bedingungen an der ersten Grenzfläche. Da der Wellenvektor der transmittierten Welle in positive x-Richtung zeigt, wird die Amplitude an der zweiten Grenzfläche durch eine positive Phasenverschiebung $k_{1}d_{1}$ beschrieben. Da der Wellenvektor ${\vec {k'_{1}}}$ in negative x-Richtung zeigt, kann die entsprechende Phasenverschiebung mit $-k_{1}d_{1}$ angeschrieben werden. In beiden Fällen kennzeichnet $d_{1}$ die Schichtdicke des Dünnfilmes. Mit Hilfe dieser Überlegungen gelangt man zu zwei Gleichungen für die zweite Grenzfläche.

    $E_{1}e^{\mathrm {i} k_{1}d_{1}}+E_{1}'e^{-\mathrm {i} k_{1}d_{1}}=E_{2}$
    $n_{1}E_{1}e^{\mathrm {i} k_{1}d_{1}}-n_{1}E_{1}'e^{-\mathrm {i} k_{1}d_{1}}=n_{2}E_{2}$

Man schreibt nun beide Gleichungen mit Imaginär- und Realteil an:

    $(E_{1}+E_{1}')\cos(k_{1}d_{1})+\mathrm {i} (E_{1}-E_{1}')\sin(k_{1}d_{1})=E_{2}$
    $n_{1}(E_{1}-E_{1}')\cos(k_{1}d_{1})+\mathrm {i} n_{1}(E_{1}+E_{1}')\sin(k_{1}d_{1})=n_{2}E_{2}$

Mit den eingangs erläuterten Stetigkeitsbedingungen an der ersten Grenzfläche kann man die beiden letzten Gleichungen in Matrix-Form anschreiben.

    $M{\binom {E_{0}+E_{0}'}{n_{0}(E_{0}-E_{0}')}}={\binom {E_{2}}{n_{2}E_{2}}}$

und

    $M={\begin{pmatrix}\cos(k_{1}d_{1})&{\frac {\mathrm {i} }{n_{1}}}\sin(k_{1}d_{1})\\\mathrm {i} n_{1}\sin(k_{1}d_{1})&\cos(k_{1}d_{1})\end{pmatrix}}$

Dividiert man schlussendlich noch durch $E_{0}$ und invertiert die Matrix $M$ zu $M'$ so gelangt man zu folgenden Ausdrücken. Die erste Gleichung enthält den Reflexionskoeffizienten $r={\frac {E'_{0}}{E_{0}}}$, sowie den Transmissionskoeffizienten $t={\frac {E_{2}}{E_{0}}}$.

    ${\binom {1}{n_{0}}}+{\binom {1}{-n_{0}}}r=M'{\binom {1}{n_{2}}}t$
    $M'={\begin{pmatrix}\cos(k_{1}d_{1})&-{\frac {\mathrm {i} }{n_{1}}}\sin(k_{1}d_{1})\\-\mathrm {i} n_{1}\sin(k_{1}d_{1})&\cos(k_{1}d_{1})\end{pmatrix}}$

Man kann somit sowohl den Transmissions- als auch Reflexionskoeffizienten für einen Dünnfilm auf einem Substrat berechnen. Um die Reflektivität eines Bragg-Spiegels berechnen zu können muss man jedoch mehrere Lagen berücksichtigen. Dies kann man tun, indem man in die obige Gleichung eine Matrix für jeden Dünnfilm einsetzt. Anschaulich gesprochen bedeutet dies, dass man ausgehend von der n-ten Schicht die Bedingungen an der ersten Grenzfläche ermitteln kann. Dies ist in der Abbildung rechts verdeutlicht. Man nennt dieses Verfahren Matrix-Transfer-Methode.

    ${\binom {1}{n_{0}}}+{\binom {1}{-n_{0}}}r=M_{1}'M_{2}'\ldots M_{N-1}'M_{N}'{\binom {1}{n_{s}}}t$

Letztendlich gelangt man mit der obigen Formel zu der maximalen Reflektivität des Bragg-Spiegels $R=r^{2}$ bei der Wellenlänge $\lambda _{0}$. Man bedient sich dabei der Tatsache, dass bei der Wellenlänge $\lambda _{0}$ alle $kd$ Terme den Wert von ${\frac {\pi }{2}}$ haben.

    $R=\left[{\frac {n_{o}(n_{2})^{2N}-n_{s}(n_{1})^{2N}}{n_{o}(n_{2})^{2N}+n_{s}(n_{1})^{2N}}}\right]^{2}$

### Bandbreite

Im vorigen Abschnitt wurde erläutert, dass der Reflexionskoeffizient an der obersten Schicht des Bragg-Spiegels mit Hilfe der Matrix-Transfer-Methode berechnet werden kann. Durch Umformung der Matrix-Transfer Gleichung des vorigen Kapitels gelangt man zu folgendem Ausdruck, der die Amplituden der magnetischen und elektrischen Feldstärke bei $z=0$ mit den Amplituden bei $z=d_{1}+d_{2}=d$ miteinander verknüpft. Hier sind die Dicken der ersten beiden Lagen des Bragg-Spiegels zu $d$ zusammengefasst. Gemäß den Ausführungen im letzten Abschnitt ist die Transfer-Matrix durch $M_{d}=M'_{1}M'_{2}$ gegeben.

    ${\binom {E(0)}{H(0)}}=M_{d}{\binom {E(d)}{H(d)}}$

Um weiter fortzufahren, bedient man sich des Bloch-Theorems. Die Amplituden einer elektromagnetischen Welle in einem Medium mit periodisch variierendem Brechungsindex mit Periodenlänge $d$ sind bei $z=d$ durch deren Werte bei $z=0$ gegeben, multipliziert mit einem Phasenfaktor $e^{iKd}$. Der Wellenvektor der sich im Bragg-Spiegel ausbreitenden elektromagnetischen Welle ist durch $K$ gegeben. Man gelangt schlussendlich zu folgendem Ausdruck:

    ${\binom {E(0)}{H(0)}}=e^{-\mathrm {i} Kd}{\binom {E(d)}{H(d)}}$

Kombiniert man die ersten beiden Gleichungen dieses Abschnitts, so gelangt man zu folgender Eigenwertgleichung:

    $M_{d}{\binom {E(d)}{H(d)}}=e^{-\mathrm {i} Kd}{\binom {E(d)}{H(d)}}$

Ziel ist es nun, einen Ausdruck für den Wellenvektor $K$ zu finden. Wenn dieser imaginäre Werte annimmt, so fallen die Komponenten der elektromagnetischen Welle exponentiell mit der Schichtdicke ab. Das heißt, dass sich die elektromagnetische Welle nicht im Spiegel ausbreiten kann. Der Bereich in dem $K$ imaginär ist, entspricht jenen Wellenlängenbereichen, in denen Bragg-Spiegel ihre Reflexionsmaxima zeigen. Um nun also einen Ausdruck für $K$ in Abhängigkeit von der Wellenlänge zu finden, bedient man sich erst der Tatsache, dass die Determinante von $M_{d}$ 1 ist. Die Determinante einer Matrix ist zugleich aber auch durch das Produkt ihrer Eigenwerte gegeben. Ein Eigenwert ist bereits durch $e^{-\mathrm {i} Kd}$ gegeben, was bedingt, dass der verbleibende $e^{iKd}$ gegeben sein muss. Schlussendlich ist die Spur der Matrix $M_{d}$ durch die Summe der Eigenwerte definiert. Unter Verwendung der Definition des Kosinus mit imaginären Exponentialfunktionen erhält man letztendlich folgende Gleichung:

    $e^{-iKd}+e^{iKd}=2\cos(Kd)=2\cos(k_{1}d_{1})\cos(k_{2}d_{2})-{\Bigl{\frac {n_{2}}{n_{1}}}+{\frac {n_{1}}{n_{2}}}{\Bigr )}\sin(k_{1}d_{1})\sin(k_{2}d_{2})$

Durch weiteres Umstellen und Ausdrücken der Wellenzahlen durch die Wellenlänge $\lambda$ erhält man weiters:

    $K=\arccos \left(\cos(k_{1}d_{1})\cos(k_{2}d_{2})-{\frac {1}{2}}\left({\frac {n_{2}}{n_{1}}}+{\frac {n_{1}}{n_{2}}}\right)\sin(k_{1}d_{1})\sin(k_{2}d_{2})\right)$
    $K=\arccos \left(\cos \left({\frac {2\pi n_{1}d_{1}}{\lambda }}\right)\cos \left({\frac {2\pi n_{2}d_{2}}{\lambda }}\right)-{\frac {1}{2}}\left({\frac {n_{2}}{n_{1}}}+{\frac {n_{1}}{n_{2}}}\right)\sin \left({\frac {2\pi n_{1}d_{1}}{\lambda }}\right)\sin \left({\frac {2\pi n_{2}d_{2}}{\lambda }}\right)\right)$

Die letzte Gleichung beschreibt nun den Wellenvektor $K$ in Abhängigkeit von der Wellenlänge des senkrecht einfallenden Lichtes. Ein Plot davon ist in der Abbildung rechts zu sehen, wobei $K$ auf der x-Achse in Einheiten von ${\frac {\pi }{d}}$ aufgetragen ist. Bei ${\frac {\pi }{d}}$ nimmt $K$ komplexe Werte an, was den Peaks im Reflexionsspektrum des Bragg-Spiegels entspricht. Passend dazu sind im Bild rechts auch die Reflexionsspektren von dementsprechenden Bragg-Spiegeln mit mehreren Dünnfilm-Lagen aufgetragen. Je mehr Lagen der Spiegel hat, umso besser entspricht er einem idealen Bragg-Spiegel und das Stoppband stimmt besser mit dem Bereich von komplexen Wellenvektoren überein. Weiters sieht man wie im vorigen Abschnitt hergeleitet, dass die Reflektivität mit der Anzahl der Dünnfilmschichten zunimmt.

Für einen idealen Bragg-Spiegel mit unendlich vielen Dünnfilm-Lagen lässt sich auch ein analytischer Ausdruck für die Breite des Stoppbands finden. Man setzt dazu das Argument des Arkuskosinus-Terms -1 oder 1, da dieser ab diesen Werten nicht mehr im reellen Bereich definiert ist. Weiters führt man die Phasendifferenz $\delta _{e}$ ein. Somit kann man die Sinus- und Kosinus-Terme zusammenfassen, da deren Argumente in einem Bragg-Spiegel die gleichen Werte haben.

    $-1=\cos ^{2}(\delta _{e})-{\frac {1}{2}}{\Bigl{\frac {n_{2}}{n_{1}}}+{\frac {n_{1}}{n_{2}}}{\Bigr )}\sin ^{2}(\delta _{e})$

Durch Umstellen gelangt man zu:

    $\cos ^{2}(\delta _{e})={\Bigl{\frac {n_{1}-n_{2}}{n_{1}+n_{2}}}{\Bigl )}^{2}$

Weiters kann man $\delta _{e}$ mit Hilfe der Hilfsvariable $\Delta g$ schreiben:

    $\delta _{e}={\frac {\pi \lambda _{0}}{2\lambda _{e}}}={\frac {\pi }{2}}(1\pm \Delta g)$

Mit Hilfe der letzten Gleichung gelangt man zu:

    $\cos ^{2}(\delta _{e})=\sin ^{2}(\pm {\frac {\pi \Delta g}{2}})$

Schlussendlich erhält man einen Ausdruck für die Breite des Stoppbands. Anstatt der Wellenlänge wie bisher wird das Stoppband durch Frequenzen ausgedrückt. Die zentrale Frequenz des Stoppbands ist $f_{0}$ während $\Delta f$ die Breite des Stoppbands charakterisiert.

    ${\frac {\Delta f}{f_{0}}}={\frac {4}{\pi }}\arcsin \left\vert {\frac {n_{1}-n_{2}}{n_{1}+n_{2}}}\right\vert$

## Herstellung

Die alternierenden Lagen von Dielektrika eines Bragg-Spiegels können mit unterschiedlichen Beschichtungsverfahren hergestellt werden, zum einen durch physikalische Gasphasenabscheidung wie Sputtern oder Aufdampfen zum anderen chemische Gasphasenabscheidung oder mit Hilfe der Sol-Gel-Methode. Eine weitere Methode um Bragg-Spiegel herzustellen ist das elektrochemische Porösizieren von Silizium-Wafern. Dabei kann man die Porosität maßgenau einstellen. Variiert man die Porosität zwischen hoher und niedriger Porosität, so erhält man eine Abfolge von Schichten mit niedrigem und hohem Brechungsindex. Das elektrochemische Porösizieren erlaubt im Gegensatz zu den vorher genannten Methoden auch die einfache Realisierung von Spiegeln mit stetig variierenden (z. B. sinusförmigen) Brechungsindex-Profilen. Solche Spiegel werden Rugate-Filter genannt.

## Anwendung

Bragg-Spiegel werden bei vielen Halbleiterlasern wie Oberflächenemittern (VCSEL), optisch gepumpten Halbleiterlasern (VECSEL), Laserdioden, DFB- und DBR-Lasern eingesetzt. Bei vielen Lasern werden Bragg-Spiegel als Spiegel verwendet, da die Wellenlänge meist genau festgelegt ist. Somit kann man mit Bragg-Spiegeln deutlich höhere Reflektivitäten erreichen als mit metallischen Spiegeln. Außerdem lassen sich Bragg-Spiegel als dichroitische Spiegel verwenden, die eine Farbe fast vollständig reflektieren und andere Farben annähernd vollständig transmittieren. Durch die Verwendung von λ/2- anstelle von λ/4-Schichten ergibt sich ein Interferenzfilter und bei Verwendung von dielektrischen Materialien ein dielektrisches Filter.

Bragg-Spiegel lassen sich auch gut in Glasfasern integrieren, wobei man von Faser-Bragg-Gittern spricht. Hier gelten die gleichen Gesetzmäßigkeiten wie auch bei anderen Bragg-Spiegeln.

Neben den bisher beschriebenen Anwendungsfeldern ist die potentielle Anwendung von porösen Bragg-Spiegeln ein aktueller Forschungsgegenstand. Mögliche Anwendungsgebiete gibt es im Bereich der analytischen Chemie sowie in der Gassensorik.

