Practical Vectorial Mode Solver for Dielectric Waveguides Based on Finite Differences

Titel: Praktische Vectoriële Modusoplosser voor Dielektrische Golfgeleiders Gebaseerd op Eindige Verschillen

1. Het Probleem

De analyse van elektromagnetische golfvoortplanting in dielektrische golfgeleiders is fundamenteel voor de moderne fotonica en telecommunicatie. Hoewel er diverse numerieke methoden bestaan (zoals de eindige-elementenmethode en eindige-verschillenmethode), blijven er uitdagingen bestaan bij het oplossen van de volledige vector-golfvergelijking:

• Spurious oplossingen: Numerieke methoden kunnen vaak kunstmatige, fysisch onjuiste oplossingen genereren.
• Randvoorwaarden: Het correct afdwingen van continuïteit en randvoorwaarden op interfaces waar de permittiviteit ($\varepsilon$) discontinu is, is moeilijk.
• Complexiteit: Bestaande hoog-accurate oplossers zijn vaak complex en rekenintensief, wat ze minder geschikt maakt voor snelle exploratieve analyse of integratie met datagedreven benaderingen.

De specifieke uitdaging waar deze paper op focust, is het vinden van een numerieke oplossing voor de geleide modi in longitudinally uniforme golfgeleiders die optisch lineair en niet-magnetisch zijn, zonder de complexiteit van geavanceerde hoog-orde schema's, maar wel met voldoende nauwkeurigheid voor praktische toepassingen.

2. Methodologie

De auteur presenteert een numerieke solver gebaseerd op de eindige-verschillenmethode (Finite Difference Method - FDM) die specifiek is ontworpen voor de elektrische veldformulering van de vector-golfvergelijking.

• Fundamentele Vergelijkingen:
  De methode start met de vector-golfvergelijking voor het elektrische veld $\mathbf{E}$:
  $$\nabla^2 \mathbf{E} + \nabla \left( \frac{1}{\varepsilon} \nabla \varepsilon \cdot \mathbf{E} \right) + k_0^2 \varepsilon_r \mathbf{E} = 0$$
  In plaats van alleen de transversale componenten ($E_x, E_y$) te gebruiken en $E_z$ achteraf te berekenen (wat fouten kan introduceren bij permittiviteitsdiscontinuïteiten), worden alle drie componenten ($E_x, E_y, E_z$) simultaan opgelost.

• Koppeling en Eigenwaardeprobleem:
  De vergelijkingen worden gekoppeld tot een systeem dat leidt tot een genormaliseerd eigenwaardeprobleem. Omdat de voortplantingsconstante $\beta$ zowel lineair als kwadratisch in de vergelijkingen voorkomt, wordt het probleem getransformeerd naar een lineair veralgemeend eigenwaardeprobleem door een hulpvariabele $F = \beta E$ in te voeren. Dit resulteert in een vergrote matrixsysteem dat opgelost kan worden met standaard lineaire algebra-routines.

• Discretisatie:

  • Er wordt een rooster gebruikt waarbij de veldcomponenten op de randen van de roostercellen worden geëvalueerd en de relatieve permittiviteit ($\varepsilon_r$) in het midden.
  • Voor ruimtelijke afgeleiden wordt een vierde-orde centraal verschil-schema gebruikt (in plaats van het standaard eerste-orde schema) om de nauwkeurigheid te verhogen:
    $$\frac{\partial A}{\partial \rho} \approx \frac{-A(\rho+2h) + 8A(\rho+h) - 8A(\rho-h) + A(\rho-2h)}{12h}$$
  • De magnetische veldcomponenten ($H_x, H_y, H_z$) worden achteraf berekend uit de gevonden elektrische velden via Maxwell's vergelijkingen, wat nauwkeuriger is dan het gebruik van de afgeleide van $\varepsilon$ in de $E_z$-formulering.

• Oplossing:
  Het probleem wordt opgelost met de eigs-functie in MATLAB (geoptimaliseerd voor grote, schaarse matrices). Alleen fysisch zinvolle modi (met een effectieve brekingsindex $n_{eff}$ kleiner dan de index van de omringende ring) worden behouden.

3. Belangrijkste Bijdragen

• Zelfconsistent Formulier: Door alle drie de elektrische veldcomponenten in één formulering op te nemen, worden randvoorwaarden op permittiviteitsinterfaces beter afgedwongen en worden numerieke artefacten (zoals spurious modes) geminimaliseerd.
• Praktische Implementatie: De methode biedt een evenwicht tussen eenvoud en nauwkeurigheid. Het is minder complex dan hoog-geordende full-vector solvers, maar nauwkeuriger dan benaderende methoden die de longitudinale component verwaarlozen.
• Toepasbaarheid: De solver is ontworpen als een praktisch, betrouwbaar instrument voor exploratieve analyse en kan worden geïntegreerd in datagedreven workflows, zonder de noodzaak van machine-precisie (wat vaak overkill is voor ontwerpfasen).

4. Resultaten

De solver is gevalideerd aan de hand van twee voorbeeldgolfgleiders, vergeleken met de commerciële solver Tidy3D (die als "ground truth" fungeerde met een zeer fijn rooster van 80 punten per golflengte):

1. Rectangulaire Si3N4-golfgeleider:

   • Een golfgeleider van $1.6 \, \mu m \times 0.7 \, \mu m$ omgeven door SiO2.
   • De berekende effectieve brekingsindex ($n_{eff}$) voor de eerste vier modi toonde een afwijking van minder dan 0,23% ten opzichte van Tidy3D.
   • De veldverdelingen ($|E_x|, |H_y|$, etc.) kwamen overeen met de verwachte fysische gedragingen (concentratie in de kern).

2. Trapeziumvormige golfgeleider (Lithium Niobate-achtig):

   • Een golfgeleider met variërende breedte (van $1.2 \, \mu m$ naar $1.0 \, \mu m$) op een dunne film.
   • Ook hier was de afwijking in $n_{eff}$ voor de eerste vier modi kleiner dan 1%.
   • De solver slaagde erin de complexe veldverdelingen bij de schuine randen correct te modelleren.

5. Betekenis en Conclusie

De studie demonstreert dat een goed geïmplementeerde eindige-verschillenmethode, die de volledige vectoriële aard van het elektrische veld respecteert, een krachtig en efficiënt alternatief is voor complexe hogere-orde methoden.

De belangrijkste waarde van dit werk ligt in het bieden van een interpreteerbare en flexibele tool voor de modale analyse van dielektrische golfgeleiders. Hoewel het niet de uiterste precisie van machine-geoptimaliseerde solvers nastreeft, biedt het voldoende nauwkeurigheid voor de meeste ontwerptaken en biedt het een transparante basis voor verdere uitbreidingen of integratie met machine learning. De methode lost het probleem van spurious oplossingen en randvoorwaarden op een elegante manier op door de volledige koppeling van de veldcomponenten binnen een gestructureerd rooster.