Was ist die Transfermatrix-Methode?

1. Einleitung

Die Transfermatrixmethode (TMM) ist eine schnelle und rechnerisch effiziente Technik, die zur Modellierung der Ausbreitung elektromagnetischer Wellen in mehrschichtigen, dünnschichtigen optischen Strukturen verwendet wird. Sie findet breite Anwendung in der Optik, Photonik und bei der Simulation von Halbleiterbauelementen, um die wellenlängenabhängige Reflexion, Transmission, Absorption, Interferenz und die Verteilung des elektrischen Feldes in optischen Schichtsystemen zu berechnen, die aus geschichteten Materialien mit unterschiedlichen Brechungsindizes bestehen. Zu den Anwendungsbereichen gehören Dünnschicht-Solarzellen, OLED-Mikrohohlräume, dielektrische Spiegel, Antireflexbeschichtungen, verteilte Bragg-Reflektoren (DBRs), optische Filter, Wellleiterbeschichtungen und andere planare photonische Bauelemente.

Die TMM-Methode basiert direkt auf den Maxwell-Gleichungen und bietet eine vereinfachte, eindimensionale Lösung der elektromagnetischen Wellengleichung für geschichtete Medien. Anstatt das vollständige Vektor-elektromagnetische Feld in einem mehrdimensionalen Simulationsbereich zu berechnen, stellt die Methode das Licht innerhalb jeder Schicht als sich vorwärts und rückwärts ausbreitende Wellen dar, die sich durch ein komplexes Brechungsindexprofil bewegen.

Da die Transfermatrixmethode das optische Problem auf eine Reihe von kompakten Matrixoperationen reduziert, ist sie im Vergleich zu vollwertigen numerischen Methoden, wie zum Beispiel, extrem schnell. FDTDObwohl das TMM-Verfahren (Transfer Matrix Method) relativ einfach zu berechnen ist, erfasst es die kohärenten Interferenzphänomene dünner Schichten präzise, einschließlich der Bildung stehender Wellen, der Fabry-Perot-Resonanzen, des Verhaltens optischer Hohlräume sowie der Feldverstärkung in dielektrischen Schichtsystemen und mehrschichtigen optischen Beschichtungen. Dies macht die Methode zu einem der am weitesten verbreiteten Ansätze zur Modellierung der kohärenten Optik in geschichteten Dünnschichtbauelementen.

2. Physikalische Grundlage der Transfermatrix-Methode

Die Transfermatrixmethode stellt das optische Feld innerhalb jeder Schicht als die Superposition einer sich vorwärts und einer sich rückwärts ausbreitenden elektromagnetischen Welle dar:

\[ E(z)=E^{+}e^{-jkz}+E^{-}e^{jkz} \]

wobei \(E^{+}\) und \(E^{-}\) die Amplituden der sich vorwärts bzw. rückwärts ausbreitenden Wellen sind, \(z\) die Ausbreitungsrichtung durch die mehrschichtige Struktur angibt und

\[ k=\frac{2\pi \tilde{n}}{\lambda}, \qquad \tilde{n}=n+j\kappa \]

Hierbei handelt es sich um den komplexen optischen Wellenvektor, wobei \(n\) der Brechungsindex ist, \(\kappa\) die optische Absorption beschreibt und \(\lambda\) die Wellenlänge im Vakuum darstellt. Wie in ?? dargestellt, sind die sich vorwärts und rückwärts ausbreitenden Wellen an jeder Materialgrenzfläche durch partielle Reflexion und Transmission miteinander gekoppelt.

Diese Wellendarstellung wird durch die Lösung der eindimensionalen Helmholtz-Wellengleichung erhalten, die aus den Maxwell-Gleichungen in einem homogenen Medium abgeleitet wird:

\[ \frac{d^2E(z)}{dz^2}+k^2E(z)=0 \]

An jeder Grenzfläche werden Maxwells elektromagnetische Randbedingungen angewendet, um die Kontinuität der tangentialen elektrischen und magnetischen Felder sicherzustellen. Diese Bedingungen führen zu den Fresnel-Reflexions- und -Transmissionskoeffizienten:

\[ r=\frac{n_1-n_2}{n_1+n_2}, \qquad t=\frac{2n_1}{n_1+n_2} \]

wobei \(n_1\) und \(n_2\) die komplexen Brechungsindizes der benachbarten Schichten sind. Diese Koeffizienten bestimmen, wie die sich vorwärts und rückwärts ausbreitenden Wellen an jeder Materialgrenzfläche teilweise reflektiert und durchgelassen werden.

3. Erstellung der Übergangsmatrix

Wie in ?? dargestellt, wird das optische Feld innerhalb jeder Schicht durch eine sich vorwärts ausbreitende Welle \(E^+\) und eine sich rückwärts ausbreitende Welle \(E^-\) dargestellt:

\[ \mathbf{E}= \begin{bmatrix} E^{+} \\ E^{-} \end{bmatrix} \]

Während sich die elektromagnetische Welle durch eine Schicht der Dicke \(d\) ausbreitet, akkumuliert sich das optische Feld, was zu einer Phasenverschiebung und Dämpfung führt. Diese Ausbreitung durch die Schicht lässt sich wie folgt darstellen:

\[ \begin{bmatrix} E^{+}(d) \\ E^{-}(d) \end{bmatrix} = \begin{bmatrix} e^{-jkd} & 0 \\ 0 & e^{jkd} \end{bmatrix} \begin{bmatrix} E^{+}(0) \\ E^{-}(0) \end{bmatrix} \]

An der Grenze zwischen zwei Materialien wird ein Teil des optischen Feldes durchgelassen und ein Teil reflektiert. Diese Beziehung an der Grenzfläche lässt sich auch in Matrixform darstellen:

\[ \begin{bmatrix} E^{+}_2 \\ E^{-}_2 \end{bmatrix} = \frac{1}{t} \begin{bmatrix} 1 & r \\ r & 1 \end{bmatrix} \begin{bmatrix} E^{+}_1 \\ E^{-}_1 \end{bmatrix} \]

Die vollständige, mehrschichtige optische Struktur wird daher als eine Folge von Propagations- und Grenzflächenmatrizen dargestellt:

\[ \mathbf{E}_{out} = M_N I_N \cdots M_2 I_2 M_1 I_1 \mathbf{E}_{in} \]

wobei \(M\) die Ausbreitung durch eine Schicht und \(I\) die Übertragung und Reflexion an einer Grenzfläche darstellt. Die optische Reaktion des gesamten Geräts kann daher mithilfe einer Reihe einfacher Matrixmultiplikationen berechnet werden.

Dieser Ansatz ermöglicht es, die Reflexion, Transmission, Absorption, Interferenz und die Bildung stehender Wellen in komplexen, mehrschichtigen optischen Strukturen äußerst effizient zu berechnen, wobei gleichzeitig eine direkte Verbindung zu den Maxwell-Gleichungen erhalten bleibt. Da die Transfermatrixmethode die optische Wellengleichung direkt löst, erfasst sie auf natürliche Weise die Interferenzeffekte dünner Schichten. Wenn Licht von mehreren Grenzflächen reflektiert wird, interferieren die reflektierten Wellen konstruktiv oder destruktiv, je nach der optischen Phase, die sich in jeder Schicht aufbaut. \( \phi = \frac{2\pi n d}{\lambda} \)wo \(n\) ist der Brechungsindex. \(d\) die Schichtdicke sowie \(\lambda\) die Wellenlänge. Infolgedessen können selbst geringfügige Änderungen der Dicke im Nanometerbereich die optischen Eigenschaften eines Bauelements erheblich beeinflussen. Diese Interferenzphysik bestimmt die Funktionsweise von Antireflexbeschichtungen, Fabry-Pérot-Resonatoren, verteilten Bragg-Reflektoren (DBRs), optischen Filtern, OLED-Mikroresonatoren und Dünnschicht-Solarzellen.

4. Was berechnet die Transfermatrixmethode?

Sobald die globale Transfermatrix der mehrschichtigen Struktur erstellt wurde, kann das elektromagnetische Feld an jeder beliebigen Position innerhalb des Bauelements berechnet werden. Dies ermöglicht die Bestimmung der räumlichen optischen Feldverteilung, der Reflexion, der Transmission und der Absorption innerhalb der gesamten Struktur.

Da die Transfermatrixmethode die optische Wellengleichung direkt löst, ergeben sich bei der Simulation auf natürliche Weise kohärente Interferenzerscheinungen. Wenn sich Licht mehrfach zwischen den Grenzflächen reflektiert, entstehen durch konstruktive und destruktive Interferenz stehende Wellenmuster im Inneren des Bauelements.

Die in ?? dargestellte Verteilung der Photonenflussdichte veranschaulicht, wie sich die Intensität des optischen Feldes innerhalb einer mehrschichtigen Struktur verändert. Bereiche mit hoher optischer Feldintensität entsprechen einer erhöhten Photonenflussdichte, die durch optische Interferenz und Kavitätseffekte verursacht wird.

Die lokale Verteilung des optischen Feldes bestimmt direkt, wo die Photonen innerhalb des Bauelements absorbiert werden. Durch die Kombination des elektromagnetischen Feldes mit dem Materialabsorptionskoeffizienten kann das optische Erzeugungsprofil berechnet werden:

\[ G(z)=\alpha(z)\Phi(z) \]

wobei \(G(z)\) die lokale Ladungsträgererzeugungsrate, \(\alpha(z)\) der optische Absorptionskoeffizient und \(\Phi(z)\) der lokale Photonenfluss ist. Ein Beispiel für ein mit der Transfermatrixmethode berechnetes Erzeugungsprofil ist in ?? dargestellt.

Diese Berechnungen werden in großem Umfang in Dünnschicht-Solarzellen, OLEDs, optischen Beschichtungen, dielektrischen Spiegeln, optischen Sensoren und resonanten Mikrohöhlenbauelementen eingesetzt. Da die Transfermatrixmethode wenig Rechenaufwand erfordert, lassen sich in der Regel vollständige optische Spektren und interne Feldverteilungen innerhalb von Sekundenbruchteilen berechnen.

5. Vorteile und Einschränkungen der Transfermatrix-Methode

Die Transfermatrixmethode (TMM) wird häufig verwendet, da sie eine hohe physikalische Genauigkeit mit einem extrem geringen Rechenaufwand verbindet. Bei geschichteten optischen Strukturen kann die TMM die gleiche physikalische Grundlage für die Interferenz dünner Schichten simulieren wie deutlich aufwendigere, vollständige Wellensimulationen, wobei sie jedoch um Größenordnungen weniger Speicherplatz und Rechenzeit benötigt.

Zu den wichtigsten Vorteilen von TMM gehören:

• Sehr schnelle Simulationen: ideal für die Durchführung von Parameterstudien und zur Optimierung von Bauelementen.
• Präzise Berechnung der Interferenz dünner Schichten: Erfasst auf zuverlässige Weise kohärente Interferenz- und stehende Welleneffekte.
• Direkte optische Ausgangssignale: Das elektrische Feld, die Photonenanzahl und die Absorptionsprofile lassen sich einfach ermitteln.
• Hervorragend geeignet für planare Multischichten: besonders nützlich für OLEDs, Solarzellen, optische Beschichtungen und DBRs.
• Einfache Anbindung an elektrische Modelle: optische Erzeugungsprofile können direkt in Drift-Diffusions-Simulationen eingegeben werden.

TMM weist jedoch auch wichtige Einschränkungen auf. Die Methode geht grundsätzlich davon aus, dass die optische Struktur näherungsweise als eindimensional behandelt werden kann, wobei sich die Materialeigenschaften hauptsächlich entlang der Wachstumsrichtung ändern. Infolgedessen ist es schwierig, beliebige dreidimensionale Streugeometrien präzise zu modellieren.

• Im Wesentlichen eindimensional: auf geschichtete Strukturen beschränkt.
• Eingeschränkte Streuphysik: Diffuse Streuung und zufällige Rauheit werden nicht umfassend berücksichtigt.
• Es dürfen keine willkürlichen Geometrien verwendet werden: Gitterstrukturen, Nanopartikel, strukturierte Oberflächen und photonische Kristalle erfordern in der Regel die Verwendung von vollwertigen Wellenlösungsverfahren.
• Kohärenzannahmen: Das Standard-TMM-Modell geht von einer kohärenten Ausbreitung innerhalb jeder Schicht aus.

Bei komplexen, dreidimensionalen optischen Problemen kommen numerische Methoden wie zum Beispiel FDTDIn den meisten Fällen sind Finite-Elemente-Methoden (FEM) oder eine präzise Analyse gekoppelter Wellen (RCWA) besser geeignet. Daher eignet sich die TMM am besten für planare Mehrschichtstrukturen, bei denen Interferenzeffekte eine wichtige Rolle spielen, die Struktur aber im Wesentlichen eindimensional bleibt.

8. Zusammenfassung

Die Transfermatrixmethode ist eine der wichtigsten numerischen Methoden in der Optik dünner Schichten. Sie bietet eine schnelle und physikalisch fundierte Möglichkeit, die Reflexion, Transmission, Interferenz, Absorption und optische Erzeugung innerhalb von mehrschichtigen Strukturen zu modellieren.

Obwohl die Methode hauptsächlich auf schichtförmige Geometrien beschränkt ist, bleibt sie eine der recheneffizientesten Methoden, die für die Simulation planarer optoelektronischer Bauelemente zur Verfügung stehen.

In praktischen Simulationsabläufen für elektronische Bauelemente wird die TMM-Methode häufig zusammen mit elektrischen Transportmodellen eingesetzt, um die gekoppelte Wechselwirkung zwischen Optik und Halbleiterphysik zu untersuchen.