And Yet Another FEM-Based Mode Solver for Dielectric Waveguides

Dit artikel presenteert een open-source, full-vector FEM-modesolver voor diëlektrische golfgeleiders, gebaseerd op een gemengde Nedelec-Lagrange-discretisatie die nauwkeurige modellering van hybride modi mogelijk maakt en spurious oplossingen onderdrukt, met validatie die een afwijking van minder dan 0,05% ten opzichte van COMSOL Multiphysics aantoont.

De kern

Stel je voor dat je een heel klein, onzichtbaar lichtpad wilt bouwen voor een computerchip. Dit pad, een golfgids, moet licht van de ene kant naar de andere sturen zonder dat het verdwaalt of verdampt. Om dit te doen, moeten ingenieurs precies weten hoe het licht zich gedraagt binnen dit pad.

Deze paper is als het ware een nieuwe, open-source bouwpakket voor een computerprogramma dat precies uitrekent hoe dat licht zich gedraagt. De auteur, Ergun Simsek, heeft een slimme manier bedacht om dit te doen met wiskunde, en hij heeft het gratis beschikbaar gemaakt voor iedereen.

Hier is de uitleg in simpele taal, met een paar leuke vergelijkingen:

1. Het Probleem: Licht is lastig

Licht in een chip is niet zomaar een rechte lijn. Het is een complex trillend veld dat in alle richtingen kan bewegen (horizontaal, verticaal en diep).

• De oude manier: Vroeger probeerden mensen dit te simuleren alsof het licht alleen maar in één richting trilde (zoals een slinger die alleen voor en achter gaat). Dit werkt goed voor simpele paden, maar faalt bij de moderne, complexe chips waar licht in alle richtingen beweegt.
• Het gevaar: Als je de wiskunde niet perfect doet, krijg je "spookantwoorden". Dit zijn oplossingen die er wiskundig mooi uitzien, maar die in de echte wereld onmogelijk zijn. Het is alsof je een tekening maakt van een vliegtuig dat vliegt, maar dat in de realiteit direct neerstort.

2. De Oplossing: De "Mix" van twee werelden

De auteur gebruikt een techniek genaamd FEM (Finite Element Method). Denk hierbij aan het opdelen van een complexe vorm in duizenden kleine, simpele driehoekjes (net als een mozaïek of een pixelated foto).

Het echte genie van deze paper zit in de mix:

• Voor het licht dat over de zijkanten "kruipt" (de transversale delen), gebruikt hij rand-elementen (edge elements). Stel je voor dat je een touw langs de randen van de driehoekjes legt. Dit zorgt ervoor dat het licht soepel over de randen blijft glijden zonder te "steken".
• Voor het licht dat "in" de driehoekjes zit (de lengterichting), gebruikt hij knooppunten-elementen (nodal elements). Dit zijn als het ware de hoekpunten van de driehoekjes waar je de hoogte van het licht meet.

De analogie:
Stel je voor dat je een dansgroep hebt.

• De rand-elementen zorgen ervoor dat de dansers hand in hand blijven terwijl ze rond een hoek draaien (geen loslaten, geen gaten).
• De knooppunten-elementen zorgen ervoor dat ze op de juiste hoogte springen.
  Door deze twee technieken te combineren, krijg je een dans die perfect is en geen enkele "spookdanser" (foute berekening) produceert.

3. Waarom is dit speciaal?

• Het is gratis en toegankelijk: Veel van deze gereedschappen zijn dure, gesloten software (zoals COMSOL, waartegen deze nieuwe tool het opneemt). Deze nieuwe tool is geschreven in Python en MATLAB, talen die veel studenten en onderzoekers al kennen. Je kunt het zelfs in je browser draaien (via Google Colab), dus je hoeft geen dure computer te kopen.
• Het werkt net zo goed als de dure tools: De auteur heeft zijn tool getest tegen de beroemde, dure software COMSOL. Het resultaat? Ze geven bijna exact hetzelfde antwoord (verschil is minder dan 0,05%). Het is alsof je een zelfgebouwde auto hebt die net zo snel rijdt als een Ferrari, maar dan voor een fractie van de prijs.
• Het is snel en betrouwbaar: De paper laat zien dat je kunt kiezen tussen snelheid en precisie. Wil je een heel snel antwoord? Dan gebruik je een grover raster (minder driehoekjes). Wil je het perfecte antwoord? Dan gebruik je meer driehoekjes, maar duurt het iets langer.

4. Wat betekent dit voor de wereld?

Dit is een hulpmiddel voor onderzoekers en studenten die aan geïntegreerde fotonica werken (lichtgebaseerde computers).

• Het helpt bij het ontwerpen van snellere internetverbindingen.
• Het helpt bij het bouwen van sensoren voor medische toepassingen.
• Het maakt onderwijs makkelijker, omdat studenten nu zelf kunnen experimenteren met lichtgolven zonder een dure licentie te hoeven kopen.

Kortom:
Deze paper introduceert een nieuwe, gratis "rekenmachine" voor licht. In plaats van te raden hoe licht zich gedraagt in een chip, kunnen ingenieurs nu met een open, betrouwbare tool precies zien wat er gebeurt, net zoals je met een GPS precies kunt zien welke route je moet nemen, zonder vast te lopen in verkeersopstoppingen (de spookantwoorden).

Technische samenvatting

Titel: Een andere FEM-gebaseerde modusoplosser voor diëlektrische golfgeleiders

Auteur: Ergun Simsek (University of Maryland Baltimore County)
Datum: 15 april 2026

---

1. Probleemstelling

De modale analyse van diëlektrische optische golfgeleiders is essentieel voor het ontwerp van geïntegreerde fotonische apparaten. Hoewel de Finite Element Method (FEM) een krachtige techniek is vanwege zijn flexibiliteit in geometrie en materiaaldistributie, bestaan er aanzienlijke uitdagingen bij de numerieke implementatie:

• Scalar vs. Vector: Traditionele scalar FEM-formuleringen (gebaseerd op alleen $E_z$ of $H_z$) zijn efficiënt maar onvoldoende voor sterk geleidende, anisotrope of hoog-contrast structuren waar hybride modi optreden.
• Spurious Modes (Ongewenste oplossingen): Volledige vectoriële FEM-formuleringen, die alle veldcomponenten oplossen, lijden vaak aan het ontstaan van niet-fysische, "spurious" modi. Deze vervuilen het numerieke spectrum en maken de resultaten onbetrouwbaar.
• Toegankelijkheid en Reproduceerbaarheid: Bestaande hoogwaardige oplossers zijn vaak gesloten (commercieel) of moeilijk toegankelijk voor educatieve doeleinden en cloud-platforms. Er is behoefte aan een open-source, reproduceerbare implementatie die zowel in MATLAB als Python werkt.

2. Methodologie

Het paper introduceert een volledig vectoriële FEM-modusoplosser voor golfgeleiders in het frequentiedomein, gebaseerd op een gemengde Nédélec–Lagrange-discretisatie van de Maxwell-vergelijkingen.

• Wiskundige Formulering:

  • De methode start met de Maxwell-curl-vergelijkingen voor een bronvrij, lineair en isotroop medium.
  • Het elektrische veld wordt ontbonden in een transversaal deel ($E_t$) en een longitudinaal deel ($E_z$), met een propagatiefactor $e^{-j\beta z}$.
  • In plaats van de sterke vorm direct op te lossen, wordt een zwakke vorm (weak form) afgeleid door de vergelijkingen te vermenigvuldigen met testfuncties en te integreren over de dwarsdoorsnede.
  • Dit leidt tot een gekoppeld systeem van vergelijkingen dat wordt omgezet in een veralgemeend eigenwaardeprobleem: $A(E, v) = \lambda B(E, v)$, waarbij $\lambda = \beta^2/k_0^2 = n_{eff}^2$.

• Discretisatie en Elementen:

  • Transversale componenten ($E_t$): Gebruikmakend van Nédélec randelementen (edge elements). Deze zorgen voor continuïteit van de tangentiële veldcomponenten over elementgrenzen en voldoen inherent aan de divergentievrije voorwaarde, wat cruciaal is voor het onderdrukken van spurious modi.
  • Longitudinale component ($E_z$): Gebruikmakend van Lagrange nodale elementen (P1).
  • Deze hybride formulering combineert de voordelen van beide elementtypes: nauwkeurige modellering van hybride modi en effectieve onderdrukking van niet-fysische oplossingen.

• Implementatie Details:

  • Mesh: De dwarsdoorsnede wordt getrianguleerd met lineaire driehoeken. De mesh wordt gegenereerd door punten te zaaien langs materiaalgrenzen (kern, cladding, BOX) en een Delaunay-triangulatie toe te passen.
  • Numerieke Integratie: Er wordt gebruikgemaakt van een symmetrische Gaussische integratieregel met drie punten (op de middelpunten van de randen van de referentietriehoek), wat exact is voor polynomen tot graad 2.
  • Oplosser: Het resulterende grote, dunne (sparse) eigenwaardeprobleem wordt opgelost met de shift-invert spectrale transformatie en de Implicitly Restarted Arnoldi Method (via scipy.sparse.linalg.eigs in Python en vergelijkbaar in MATLAB).
  • Randvoorwaarden: Voor open grenzen wordt $E_z = 0$ opgelegd op de buitenste knopen (absorberende voorwaarde). Voor geleidende wanden (PEC) worden de tangentiële rand-DOF's op nul gezet.

3. Belangrijkste Bijdragen

1. Open Source Implementatie: De oplossing is beschikbaar gesteld in zowel MATLAB als Python, met specifieke aandacht voor reproduceerbaarheid en compatibiliteit met cloud-platforms zoals Google Colab.
2. Hybride Discretisatie: Een robuuste implementatie van de Nédélec-Lagrange-methode die spurious modi effectief onderdrukt zonder de complexiteit van kwadratische eigenwaardeproblemen die bij sommige andere transverse formuleringen optreden.
3. Validatie en Benchmarking: Uitgebreide numerieke validatie tegen de industriestandaard COMSOL Multiphysics, met een focus op nauwkeurigheid en convergentie.
4. Educatieve Toegankelijkheid: Het biedt een transparante, code-beschreven methode voor studenten en onderzoekers om de onderliggende FEM-principes voor golfgeleiders te begrijpen, zonder afhankelijk te zijn van gesloten software.

4. Resultaten

De solver werd getest op twee representatieve golfgeleiderstructuren:

• Case Study 1 (Si3N4 golfgeleider):

  • Een rechthoekige Si3N4-golfgeleider in SiO2.
  • De resultaten toonden een uitstekende overeenkomst met COMSOL. De relatieve fouten voor de effectieve index ($n_{eff}$) lagen onder de 0,004% voor alle geanalyseerde modi.
  • Convergentiestudies toonden aan dat de fout systematisch afneemt bij verfijning van het mesh (toename van punten per golflengte, PPW).

• Case Study 2 (Hoog-index-contrast):

  • Een golfgeleider met een hoge brekingsindex ($n_{core} = 3.5$), wat leidt tot sterke veldconfinement en meer hybride modi.
  • De fouten bleven onder de 0,05% vergeleken met COMSOL.
  • Rekenkosten: De MATLAB-implementatie was iets sneller dan de Python-versie, maar beide waren iets trager dan COMSOL voor zeer fijne meshes (vanwege de geoptimaliseerde commerciële solvers van COMSOL). Echter, de nauwkeurigheid was vergelijkbaar.

• Algemene Observaties:

  • De methode convergeert zoals verwacht voor FEM-methoden.
  • Er is een duidelijke trade-off tussen rekenkosten en precisie.
  • De solver kan de polarisatie (TE/TM fractie) en modale overlappingen nauwkeurig berekenen.

5. Betekenis en Conclusie

Dit werk levert een betrouwbare, flexibele en open-source tool voor onderzoek en onderwijs in geïntegreerde fotonica. Hoewel de rekenprestaties nog niet volledig gelijk zijn aan die van geoptimaliseerde commerciële software, biedt de solver een uitstekende balans tussen nauwkeurigheid en toegankelijkheid.

De belangrijkste waarde ligt in de reproduceerbaarheid en de open aard van de code, waardoor onderzoekers en studenten de complexiteit van vectoriële FEM-modellering kunnen doorgronden en aanpassen aan specifieke golfgidsconfiguraties (zoals trapeziumvormige doorsneden of anisotrope materialen in de toekomst). De succesvolle validatie tegen COMSOL bevestigt dat de gemengde Nédélec-Lagrange-aanpak een solide basis vormt voor de analyse van moderne, hoog-contrast diëlektrische golfgeleiders.