Accuracy assessment of scalar wave propagation methods for diffractive optics design: from thin elements to thick binary grating

Dieser Beitrag bewertet systematisch die Genauigkeit der Näherung dünner Elemente sowie der Strahl- und Wellenausbreitungsmethoden anhand einer rigorosen Referenz für binäre diffraktive Gitter und erstellt Genauigkeitskarten, um die Auswahl geeigneter Vorwärtsmodelle in inversen Entwurfsprozessen auf Basis der räumlichen Frequenz und der Gitterdicke zu leiten.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, eine komplexe, unebene Oberfläche (wie ein mikroskopisches Labyrinth) zu entwerfen, die Licht auf sehr spezifische Weise so bricht, dass ein Hologramm oder eine spezielle Linse entsteht. Um dies zu tun, benötigen Sie ein Computerprogramm, das exakt vorhersagt, wie sich das Licht durch dieses Labyrinth bewegen wird.

Dieser Artikel ist wie ein „Führerschein-Test" für drei verschiedene Computerprogramme, die versuchen, die Reise des Lichts vorherzusagen. Die Autoren wollten wissen: Welches Programm liefert die genaueste Karte, und unter welchen Bedingungen versagt jedes einzelne?

Hier ist die Aufschlüsselung ihrer Erkenntnisse unter Verwendung einfacher Analogien:

Die drei „Navigatoren" (Die Methoden)

Die Forscher testeten drei verschiedene Möglichkeiten, die Bewegung von Licht durch diese mikroskopischen Strukturen zu simulieren:

1. Der „Sofort-Teleporter" (TEA – Thin-Element Approximation / Dünnelement-Näherung):

   • Funktionsweise: Diese Methode tut so, als wäre die unebene Struktur so dünn, dass sie nicht existiert. Sie berechnet einfach den Weg des Lichts, als würde das Licht sofort durch die Oberfläche „teleportieren", wobei es seine Richtung basierend auf der Form ändert, aber die Zeit ignoriert, die es benötigt, um durch das Material zu reisen.
   • Die Analogie: Es ist wie der Versuch vorherzusagen, wie ein Auto durch einen Tunnel fährt, indem man nur auf die Schilder am Eingang und Ausgang schaut und den Tunnel selbst ignoriert.
   • Das Ergebnis: Es ist superschnell und einfach, funktioniert aber nur, wenn der Tunnel sehr kurz ist. Wenn der Tunnel länger wird (dicker), wird die Vorhersage völlig falsch, weil sie die Kurven und Wendungen im Inneren vergisst.

2. Der „Geradlinige Wanderer" (BPM – Beam Propagation Method / Strahlausbreitungsmethode):

   • Funktionsweise: Diese Methode zerlegt den Tunnel in viele dünne Scheiben und berechnet das Licht schrittweise. Sie geht jedoch davon aus, dass sich das Licht hauptsächlich gerade nach vorne bewegt und nur kleine, sanfte Kurven macht.
   • Die Analogie: Stellen Sie sich vor, Sie gehen durch einen Wald. Diese Methode geht davon aus, dass Sie geradeaus laufen und nur gelegentlich ein wenig nach links oder rechts treten. Wenn der Pfad erfordert, dass Sie eine scharfe 90-Grad-Wende machen, verirrt sich dieser Wanderer, da er nicht darauf programmiert ist, große Winkel zu bewältigen.
   • Das Ergebnis: Es ist besser als die „Teleporter"-Methode für dickere Tunnel, aber wenn das Licht scharfe Kurven (große Winkel) machen muss oder der Tunnel sehr lang ist, summieren sich die kleinen Fehler in seiner „geradlinigen" Annahme auf, und die Karte wird unscharf.

3. Der „Wahre Navigator" (WPM – Wave Propagation Method / Wellenausbreitungsmethode):

   • Funktionsweise: Dies ist die anspruchsvollste der drei. Wie die zweite Methode geht sie schrittweise durch den Tunnel Scheibe für Scheibe, verwendet aber eine komplexere mathematische Formel, die jeden Winkel der Kurve zulässt, nicht nur kleine.
   • Die Analogie: Dieser Wanderer kennt die exakten Regeln der Physik. Er kann geradeaus gehen, scharf abbiegen oder sogar perfekt zickzack laufen. Er geht nicht davon aus, dass der Pfad einfach ist; er berechnet die exakte Kurve jedes Schritts.
   • Das Ergebnis: Es ist am genauesten, besonders für lange Tunnel oder Pfade mit scharfen Kurven. Es bleibt viel länger dem „wahren" Pfad treu als die anderen beiden.

Der „Goldstandard" (Die Referenz)

Um zu wissen, wer das Rennen gewonnen hat, verwendeten die Forscher eine supergenaue, leistungsstarke Methode namens FMM (Fourier Modal Method / Fourier-Modalmethode).

• Die Analogie: Denken Sie an FMM als eine Hochgeschwindigkeits-Drohne, die über den Wald fliegt, Millionen von Fotos macht, um eine perfekte 3D-Karte zu erstellen, die genau zeigt, wo jedes Blatt und jeder Ast ist. Es erfordert viel Rechenleistung und Zeit, daher würde man es nicht für jede einzelne Schätzung verwenden, aber es ist die „Wahrheit", an der die anderen drei gemessen werden.

Das Experiment: Zufällige Labyrinthe

Die Forscher testeten nicht nur ein Labyrinth. Sie generierten 1.210 zufällige mikroskopische Labyrinthe mit zwei variierenden Merkmalen:

1. Dicke: Wie tief der Tunnel ist (von 1 Schicht bis 11 Schichten dick).
2. Komplexität: Wie uneben und scharf die Kurven sind (von sanften Hügeln bis zu gezackten, scharfen Spitzen).

Sie ließen alle drei „Navigatoren" durch diese Labyrinthe laufen und verglichen ihre Karten mit der „Goldstandard"-Drohnentkarte.

Das Urteil: Wann welche Methode verwenden?

Der Artikel erstellte „Genauigkeitskarten" (wie Wetterkarten, die zeigen, wo es sicher zu fahren ist), die Ihnen sagen, welche Methode Sie wählen sollen:

• Verwenden Sie den „Sofort-Teleporter" (TEA) nur, wenn: Die Struktur extrem dünn ist (weniger als die Breite einer einzelnen Lichtwelle). Wenn sie auch nur ein wenig dicker wird, hören Sie auf, sie zu verwenden; sie wird Ihnen ein schlechtes Design liefern.
• Verwenden Sie den „Geradlinigen Wanderer" (BPM), wenn: Die Struktur dünn ist ODER wenn die Struktur dick ist, das Licht aber nur sehr sanfte Kurven machen muss. Es ist ein gutes Mittelmaß.
• Verwenden Sie den „Wahren Navigator" (WPM), wenn: Sie dicke Strukturen entwerfen, die moderate scharfe Kurven erfordern. Dies ist der Sweet Spot, in dem die anderen beiden Methoden zu versagen beginnen, WPM aber immer noch das Richtige liefert.

Der Haken

Die Forscher testeten diese Methoden an „Binären Gittern", die wie Labyrinthe mit sehr scharfen, gezackten Wänden sind (wie eine Treppe). Sie stellten fest, dass dies ein „Hard-Mode"-Test ist. Wenn Sie glattere, sanftere Strukturen entwerfen (wie einen sanften Hügel), würden alle drei Methoden wahrscheinlich noch besser abschneiden als die hier gezeigten Ergebnisse.

Kurz gesagt: Wenn Sie komplexe, dicke optische Geräte entwerfen wollen, verlassen Sie sich nicht auf die einfache „Teleporter"-Methode. Wenn die Struktur dick ist und das Licht abbiegen muss, ist der „Wahre Navigator" (WPM) der einzige, der Sie nicht verirren lässt.

Technische Zusammenfassung

Technische Zusammenfassung: Genauigkeitsbewertung skalare Wellenausbreitungsmethoden für das Design diffraktiver Optik

Problemstellung
Das Design diffraktiver optischer Elemente (DOEs) stützt sich zunehmend auf gradientenbasiertes inverses Design, wobei die Genauigkeit des Vorwärtsausbreitungsmodells die Leistung des resultierenden Designs fundamental begrenzt. Während rigorose Methoden wie die Fourier-Moden-Methode (FMM) und Finite-Difference Time-Domain (FDTD) eine hohe Genauigkeit bieten, sind sie für iterative Optimierungsloops rechnerisch prohibitiv. Umgekehrt ist die Dünnelement-Näherung (TEA) rechnerisch trivial, versagt jedoch bei dicken Strukturen oder großen Beugungswinkeln aufgrund ihrer Vernachlässigung innerer Ausbreitungseffekte und Kopplung an Grenzflächen. Intermediäre skalare Methoden, insbesondere die Strahlausbreitungsmethode (BPM) und die Wellenausbreitungsmethode (WPM), bieten einen Mittelweg, doch eine systematische, quantitative Bewertung ihrer Genauigkeitsbereiche in Abhängigkeit von Gitterparametern hat in der Literatur gefehlt.

Methodik
Die Autoren führten eine systematische Genauigkeitsbewertung unter Verwendung der rigorosen FMM (implementiert via FMMAX) als Referenzstandard durch. Die Studie nutzte folgendes experimentelles Design:

• Teststrukturen: Ensembles zufälliger binärer Gitter wurden auf einer $7.5 \, \mu\text{m}$ periodischen Zelle generiert. Die Gitter wurden durch den räumlichen Frequenzcutoff ($f_c$) und die physikalische Dicke ($h$) parametrisiert.
  • $f_c$ wurde von 0 bis $1.33 \, \mu\text{m}^{-1}$ abgetastet, was maximalen Beugungswinkeln ($\theta_{\text{max}}$) von bis zu $\approx 45^\circ$ entspricht.
  • Die Dicke ($h$) wurde von $1\lambda$ bis $11\lambda$ abgetastet (wobei $\lambda = 532 \, \text{nm}$).
  • Insgesamt wurden 1.210 Gitter (11 $\times$ 11 Parameterpaare $\times$ 10 zufällige Realisierungen) ausgewertet.
• Ausbreitungsmodelle: Drei skalare Methoden wurden im GPU-beschleunigten, differenzierbaren Framework FluxOptics.jl implementiert:
  • TEA: Modelliert als unendlich dünne Phasenmaske, wobei Ausbreitungsschritte in einem effektiven Medium approximiert werden.
  • BPM: Verwendet eine paraxiale Näherung, wobei das Gittervolumen in Schritte von $\Delta z = 25 \, \text{nm}$ unterteilt wird.
  • WPM: Setzt einen exakten nicht-paraxialen Propagator ($k_z = \sqrt{(nk_0)^2 - k_\perp^2}$) ein, mit einer festen Anzahl von 30 Schichten unabhängig von der Dicke und unter Verwendung glatter Übergangsschichten für Indexdiskontinuitäten.
• Metrik: Die Genauigkeit wurde mittels der normierten Feldüberlappung ($\eta$) zwischen dem vom skalaren Modell übertragenen Feld und dem FMM-Referenzfeld quantifiziert, gemittelt über die 10 Realisierungen für jeden Parametersatz. Kreuzpolarisierte Komponenten wurden ausgeschlossen, um konsistent mit der skalaren Näherung zu bleiben.

Hauptergebnisse
Die Studie generierte Genauigkeitskarten im $(f_c, h)$-Parameterraum und enthüllte distincte Gültigkeitsbereiche für jede Methode:

• TEA: Die Genauigkeit verschlechtert sich mit zunehmender Dicke rapide. Der Überlapp fällt bei moderaten Winkeln unter 90% bei $h \approx 3\lambda$ und sinkt unter 50% für $h > 6\lambda$, was ihre Beschränkung auf sehr dünne Strukturen bestätigt.
• BPM: Performt signifikant besser als TEA, leidet jedoch unter akkumulierten paraxialen Phasenfehlern. Während sie den gültigen Bereich auf moderate Dicken erweitert, stabilisiert sich ihre Genauigkeit bei einem Überlapp von 0,4–0,5 für große Dicken innerhalb des getesteten Bereichs.
• WPM: Behält über den gesamten Parameterraum die höchste Genauigkeit bei. Für kleine Winkel ($\theta_{\text{max}} \lesssim 12^\circ$) behält WPM $\eta > 0,90$ selbst für Dicken bis zu $11\lambda$ bei. Der Vorteil von WPM gegenüber BPM ist am ausgeprägtesten für dicke Strukturen ($h > 3\lambda$) bei moderaten Winkeln ($\theta_{\text{max}} \lesssim 30^\circ$), wobei die exakte Dispersionsrelation die fortschreitende Dephasierung verhindert, die BPM begrenzt.

Bedeutung und Behauptungen
Die Arbeit beansprucht, die erste systematische quantitative Bewertung dieser skalaren Methoden speziell für binäre diffraktive Gitter zu liefern, einen konservativen Benchmark aufgrund ihrer scharfen Indexdiskontinuitäten. Der primäre Beitrag ist die Generierung von Genauigkeitskarten, die als praktische Leitlinien für die Auswahl von Vorwärtsmodellen in inversen Design-Pipelines für diffraktive Optik dienen:

1. TEA ist nur für dünne Strukturen ($h \lesssim \lambda$) angemessen.
2. BPM ist eine wettbewerbsfähige Alternative für dünne Strukturen bei großen Winkeln oder moderate Dicken bei kleinen Winkeln.
3. WPM wird als Methode der Wahl für das inverse Design dicker diffraktiver Strukturen bei moderaten Winkeln identifiziert und bietet einen signifikanten Genauigkeitsvorteil gegenüber BPM ohne die prohibitiven Kosten rigoroser Methoden.

Die Autoren weisen darauf hin, dass binäre Gitter einen „konservativen Benchmark" darstellen (Maximierung des spektralen Inhalts an Indexdiskontinuitäten), sodass die berichteten Genauigkeitskarten wahrscheinlich Untergrenzen für die Leistung dieser Methoden bei glatten, kontinuierlich reliefierten Strukturen darstellen. Die Studie erkennt an, dass alle drei Methoden Einschränkungen bezüglich unidirektionaler Ausbreitung und der Vernachlässigung vektorieller Effekte teilen, wobei Erweiterungen auf Split-Step-nicht-paraxiale Methoden zukünftigen Arbeiten vorbehalten bleiben.