Reflektometrische Dünnschichtmessung

Die Reflektometrische Dünnschichtmessung ist ein Messverfahren, das auf den Prinzipien der Dünnschichtinterferenz beruht, und dazu dient, die Schichtdicke von dielektrischen Dünnschichten zu bestimmen. Außerdem kann aus den Daten der Brechungsindex des Materials gewonnen werden. Diese Informationen sind von zunehmender Bedeutung, da Dünnschichten heutzutage vielfach eingesetzt werden, z. B. bei Antireflexionsbeschichtungen, der Beschichtung von medizinischen Implantaten oder in der Solarindustrie.

Theorie

Als Dünnschicht wird eine Materialschicht bezeichnet, dessen zweidimensionale Ausdehnung um einige Größenordnungen höher als die Schichtdicke ist. Im Normalfall werden Dünnschichten mit einer Dicke im Nano- bis Mikrometerbereich verwendet. Einfallendes Licht wird an der Grenze Material–Umgebungsmedium (zumeist Luft) entweder reflektiert oder transmittiert; das transmittierte Licht wiederum trifft auf die Grenzfläche Dünnschicht–Substrat und wird hier wiederum reflektiert. Wie viel Licht hierbei reflektiert bzw. transmittiert wird, lässt sich mit Hilfe der Fresnel-Gleichungen bestimmen. Zwischen den an den beiden Grenzflächen reflektierten Anteilen kommt es in der Folge zu konstruktiver bzw. destruktiver Interferenz. Diese ist dabei abhängig von der Schichtdicke $d$, dem Brechungsindex $n$ des Materials sowie dem Einfallswinkel $\alpha$ des Lichts. Außerdem kann es zu einer Phasenverschiebung um 180° bzw. π kommen.

Strahlengang bei der Dünnschichtinterferenz

Der optische Weg (OPD) des einfallenden Lichts lässt sich anhand der Abbildung folgendermaßen berechnen:

$\mathrm {OPD} =n_{2}({\overline {\mathrm {AB} }}+{\overline {\mathrm {BC} }})-n_{1}({\overline {\mathrm {AD} }})$

sowie

${\overline {\mathrm {AB} }}={\overline {\mathrm {BC} }}={\frac {d}{\cos(\beta )}}$

${\overline {\mathrm {AD} }}=2d\tan(\beta )\sin(\alpha )$

Nach Snellius folgt

$n_{1}\sin(\alpha )=n_{2}\sin(\beta )$

und damit

$\mathrm {OPD} =n_{2}\left({\frac {2d}{\cos(\beta )}}\right)-2d\tan(\beta )n_{2}\sin(\beta )$

$\mathrm {OPD} =2n_{2}d\left({\frac {1-\sin ^{2}(\beta )}{\cos(\beta )}}\right)$

$\mathrm {OPD} =2n_{2}d\cos(\beta )$

Bei der reflektometrischen Schichtdickenmessung wird das Licht einer Weißlichtquelle senkrecht auf die zu messende Schicht gelenkt und mit Hilfe eines Spektrometers die Reflexion in Abhängigkeit von der Wellenlänge gemessen. Aus der obigen Gleichung ergibt sich mit $\cos(0^{\circ })=1$, dass nur für bestimmte Wellenlängen konstruktive Interferenz auftritt und diese Wellenlängen direkt abhängig sind von der Schichtdicke und dem Brechungsindex.

$\lambda ={\frac {2n_{2}d}{m}}$

So ergibt sich ein charakteristisches Oszillationsmuster, dessen Frequenz direkt von der Schichtdicke abhängt. Dieses Muster in der gemessenen Reflexion sei nun hergeleitet. Betrachtet man eine elektromagnetische Welle die auf eine Grenzfläche trifft, so wird sie teilweise reflektiert und transmittiert. Die transmittierte Welle trifft nun im Falle eines Dünnfilms wiederum auf eine Grenzfläche an der ebenfalls Transmission und Reflexion stattfinden. Dieser Vorgang wiederholt sich unendlich oft, wobei jedes Mal eine elektromagnetische Welle vom Dünnfilm reflektiert wird. Um dies mathematisch zu erfassen sei die folgende Abbildung betrachtet, in der die ersten vier Reflexionsvorgänge illustriert sind. Die Einzelreflexionen können durch die folgenden vier Gleichungen beschrieben werden, wobei $E_{i,0}$ die Amplitude des einfallenden elektrischen Feldvektors darstellt. Der Wellenvektor der reflektierten Wellen ist durch ${\vec {k_{r}}}$, der Ortsvektor durch ${\vec {r}}$ und deren Kreisfrequenz ist durch $\omega$ gegeben.

Strahlengang bei der Dünnschichtinterferenz

${E}_{r,1}=r_{01}E_{i,0}e^{{\vec {k_{r}}}{\vec {r}}-\omega t}$ 1. Reflexion

${E}_{r,2}=t_{01}r_{12}t_{10}E_{i,0}e^{{\vec {k_{r}}}{\vec {r}}-\omega t}e^{i\delta }$ 2. Reflexion

${E}_{r,3}=t_{01}r_{12}r_{10}r_{12}t_{10}E_{i,0}e^{{\vec {k_{r}}}{\vec {r}}-\omega t}e^{i2\delta }$ 3. Reflexion

${E}_{r,4}=t_{01}r_{12}(r_{10}r_{12})^{2}t_{10}E_{i,0}e^{{\vec {k_{r}}}{\vec {r}}-\omega t}e^{i3\delta }$ 4. Reflexion

An jedem Interface wird die ursprüngliche Amplitude der einfallenden Welle gemäß der Fresnelschen Gleichungen vermindert. Dem wird durch die Reflexions- und Transmissionskoeffizienten, $r$ und $t$ Rechnung getragen. Die Subskripte bezeichnen die jeweilige Richtung. So steht der Koeffizient $r_{01}$ für die Reflexion der einfallenden Welle an der 1. Grenzfläche. Man erhält somit Gleichungen für unendlich viele reflektierte Wellen $E_{r,n}$. Betrachtet man die vorgestellten 4 Gleichungen für die ersten vier Reflexionen so sieht man ein Muster. Die zweite reflektierte Welle resultiert aus einer Welle die erst an der ersten Grenzfläche transmittiert, an der zweiten reflektiert und schließlich an der ersten Grenzfläche wiederum transmittiert wird. Die resultierende Welle erfährt dabei auch eine Phasenverschiebung $\delta$ wie im Exponentialterm zu sehen ist. Die Phasenverschiebung ergibt sich aus der Wellenzahl und der optischen Weglänge. Wenn es weiters zu einem Phasensprung 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 \Delta_\text{phase}} an einer Grenzfläche kommt, so muss auch dieser berücksichtigt werden. So ein Phasensprung tritt auf, wenn es zum Übergang von niedrigem zu hohem Brechungsindex kommt und hat stets einen Wert von $\pi$. Für weitere Details hierzu siehe Fresnel-Gleichungen. Zusammenfassend erhält man für die Phasenverschiebung:

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 \delta = \frac{4\pi n_2 d \cos(\beta)}{\lambda} + \Delta_\text{phase}}

Ab der dritten Welle kommt es zusätzlich zu zwei Reflexionen an der ersten und zweiten Grenzfläche, repräsentiert durch 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 r_{10}} und 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 r_{12}} . Weiters wird die Phasenverschiebung vervielfacht. Um nur die Verminderung der Reflexion zu betrachten definiert man einen totalen Reflexionskoeffizienten 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 r_\mathrm{tot}} .

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 r_\mathrm{tot}=r_{01}+t_{01}r_{12}t_{10}e^{i \delta}+t_{01}r_{12}r_{10}r_{12}t_{10}e^{i2\delta}+t_{01}r_{12}(r_{10}r_{12})^2t_{10}e^{i3\delta} + \dots }

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 r_\mathrm{tot}=r_{01}+t_{01}r_{12}t_{10}e^{i \delta}\sum_{n=0}^\infty (r_{12}r_{10})^ne^{i\delta n}}

Da die Reflexionskoeffizienten in ihrem Betrag stets kleiner als eins sind, handelt es sich bei der Summe in der letzten Gleichung um eine geometrische Reihe. Somit kann der totale Reflexionskoeffizient wie folgt vereinfacht werden:

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 r_\mathrm{tot} = r_{01}+\frac{t_{01}r_{12}t_{10}e^{i \delta}}{1-r_{12}r_{10}e^{i\delta}}}

Somit erhält man einen Ausdruck für die gesamte reflektierte Welle 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_r} an der ersten Grenzfläche. Der Ortsvektor wird hierbei als Nullvektor angesetzt. Weiters kann man aufgrund der Zeit-Reversibilität von Wellen (Stokes Relation) die beiden Transmissionskoeffizienten an der ersten Grenzfläche mit dem dementsprechenden Reflexionskoeffizienten als 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 t_{01}t_{10}=1-{r_{01}}^2} anschreiben.

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_r = r_\mathrm{tot}E_0e^{i\omega t}=\frac{r_{01}+r_{01}e^{i\delta}}{1+r_{12}r_{01}e^{i\delta}}E_0e^{i\omega t}}

Reflektivität von Dünnfilmen mit unterschiedlicher Schichtdicke

Die Intensität einer elektromagnetischen Welle ist in Vakuum durch 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 = 0{,}5\epsilon_0cEE^*} gegeben. Hier ist 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 \epsilon_0} die Dielektrizitätskonstante von Vakuum, 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 c} die Lichtgeschwindigkeit und 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} und 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^*} die elektrische und komplex konjugierte elektrische Feldstärke. Wenn man nun die Phase Null als Vakuum (oder näherungsweise Luft) annimmt so ergibt sich für die reflektierte Intensität.

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_r = 0{,}5\epsilon_0c{E_0}^2\frac{{r_{01}}^2+{r_{12}}^2+2r_{01}r_{12}\cos(\delta)}{1+(r_{01}r_{12})^2+2r_{01}r_{12}\cos(\delta)}}

Der Vorfaktor 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 0{,}5\epsilon_0c{E_0}^2} kann als die Intensität der einfallenden Strahlung identifiziert werden. Dividiert man also die letzte Gleichung durch diesen Ausdruck, so erhält man den Reflexionsgrad 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 R} ; also den prozentuellen Anteil des Lichts der vom Dünnfilmsystem wieder reflektiert wird:

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 R=\frac{{r_{01}}^2+{r_{12}}^2+2r_{01}r_{12}\cos(\delta)}{1+(r_{01}r_{12})^2+2r_{01}r_{12}\cos(\delta)}}

Im Bild rechts ist die ermittelte Reflexion für Dünnfilmsystem mit den Brechungsindizes 1, 1,5 und 2 in aufsteigender Reihenfolge ausgerechnet. Die jeweiligen Reflexionskoeffizienten wurden mit Hilfe der Fresnelschen Formeln ermittelt, wobei wie in der reflektometrischen Dünnschichtmessung angewendet, von einem senkrecht einfallendem Lichtstrahl ausgegangen wurde. Somit ist 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 \beta} Null und die Polarisation spielt für die Berechnung keine Rolle mehr. Man erkennt, dass sich die Anzahl der Minima und Maxima erhöht, wenn man die Schichtdicke des Dünnfilmes (Phase 2) von 500 nm auf 1500 nm erhöht. Weiters ist für den Fall eines 500 nm dicken Dünnfilms zusätzlich die Intensität geplottet, die nur von den ersten beiden reflektierten Wellen verursacht wird. Die Position der Minima und Maxima bleibt unverändert, wenn man alle Reflexionsvorgänge berücksichtigt. Jeder Reflexionsvorgang verstärkt aber den Effekt positiver oder destruktiver Interferenz.

Auswertung der Messung

Auswertung der gemessenen Reflexion

In der Beschreibung des Messprinzips wurde gezeigt, dass man bei der Messung der Reflexion ein Signal mit Minima und Maxima abhängig von der Wellenlänge erhält. Die Position der Maxima der Reflektivität ist durch die optische Weglänge 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 OPD} sowie durch den Brechungsindex des Dünnfilmes gegeben. Plottet man die inverse Wellenlänge (Wellenzahl) als Funktion der Interferenzordnung 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} so erhält man Punkte, durch die man eine Gerade legen kann. Dies ist aus der im vorigen Kapitel abgeleiteten Gleichung für Interferenzmaxima ersichtlich.

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 \lambda = \frac{2n_2 d }{m} }

Für konstruktive Interferenz nimmt 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} den Wert ganzer Zahlen, für destruktive Interferenz nimmt 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} Vielfache von 0.5 an. Für die Auswertung werden beide Fälle zur Interferenzordnung zusammengefasst. Invertiert man die Bedingung für Interferenzminima und Maxima so erhält man:

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}{\lambda} = \frac{m}{2n_2 d }}

Genau dieser Zusammenhang ist in der Abbildung rechts veranschaulicht. Als Steigung erhält man den Wert 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}{2n_2d}} . Kennt man den Brechungsindex des Dünnfilms 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_2} , so kann man dessen Dicke 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 d} berechnen bzw. bei Kenntnis der Schichtdicke kann man auf den Brechungsindex schließen. Hierbei handelt es sich um ein idealisiertes Beispiel. In der Realität hängt der Brechungsindex eines Materiales mehr oder weniger stark von der Wellenlänge ab. Diesem Effekt muss gegebenenfalls Rechnung getragen werden, indem man die Messdaten mit komplizierteren Zusammenhängen fittet.

Literatur

• Eugene Hecht: Optics. Addison-Wesley, 2002, ISBN 0-321-18878-0.
• Olaf Stenzel: Optical Coatings: Material Aspects in Theory and Practice. Springer Series in Surface Sciences, 2014, ISBN 978-3-642-54062-2.